R0-II1S3  998 
UNCLASSIFIED 


DEEP  SEN  SHIP  HOOK  VOLUNE  2<U>  GENEML  ELECTRIC  CO  1/2 

SVRRCUSE  NV  ELECTRONIC  SVSTEHS  DIV  R  L  MEBSTER  OCT  78 
CHES/NflVFflC-FPO-1-78  <  1«  )  -2 

F/G  13/2  NL 


DEEP  SEA  SHIP  MOOR 
USER'S  MANUAL 
SEADYN/DSSM 

VOLUME  2 


FPO-1-78  (16) 


AUGUST  1978 


Thla  docucrant  hoi  bMO  op( 
tor  put''  ond  Botoi 

'tlatributi'ju  in  unUmltwL 


DTIC 

^ELBCTE 

^^JAN3  0e86 


Ocean  Engineering 

and 

Construction 

86  1  28  068 


REPORT  DOCUMENTATION  PAGE 

la.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 

lb.  RESTRICTIVE  MARKINGS 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3.  DISTRIBUTION  AVAILABILITY  OF  REP. 
Approved  for  public  release: 
distribution  is  unlimited 

2b.  DECLASSIFICATION/DOWNGRADING  SCHEDULE 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER 
FPO-l-78(16) 

5.  MONITORING  ORGANIZATION  REPORT  # 

V  v' *. 

y.-r.' 


6a.  NAME  OF  PERFORM.  ORG.  6b.  OFFICE  SYM 
Electronic  Systems  Div. 

General  Electric  Company 


•i 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Ocean  Engineering 

Si  Construction 

Project  Office 

CHESNAVFACENGCOM 


6c.  ADDRESS  (City.  State,  and  Zip  Code) 

Syracuse,  N.Y. 

7b.  ADDRESS  (City,  State,  and  Zip  ) 
BLDG.  212.  Washington  Navy  Yard 
Washington.  D.C.  20374-2121 

■ 

8a,  NAME  OF  FUNDING  ORG.  8b.  OFFICE  SYM 

9.  PROCUREMENT  INSTRUMENT 

INDENT  # 

8c.  ADDRESS  (City,  State  &  Zip) 

10.  SOURCE  OF  FUNDING  NUMBERS 

s'  . 

1  •  - 

PROGRAM  PROJECT  TASK 

ELEMENT  #  #  » 

WORK  UNIT 
ACCESS  # 

i.  •  • 

11.  TITLE  (Including  Security  Classification) 

Deep  Sea  Ship  Moor  User's  Manual  DESMOOR 

••.Vo 
V.  .'O 

12,  PERSONAL  AUTHOR(S) 

R.L.  Webster 

13a.  TYPE  OF  REPORT  13b.  TIME  COVERED 

FROM  TO 

14,  DATE  OF  REP.  (YYMMDD) 
76-08 

15.  PAGES 

5  vol. 

**■  •’ .  •’ 

16.  SUPPLEMENTARY  NOTATION 

'**.**'  * 

17. 


COSAT I  CODES 


FIELD 


GROUP 


SUB-GROUP 


Mooring  systems,  DESMOOR 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  &  identify  by  block  number) 

This  report  documents  the  DESMOOR  computer  program  which  selects  trial  designs  '..--I;-; 
of  mooring  systems  for  surface  ships  in  deep  water.  Either  slack  (negatively  'v-!! 
buoyant  catenary)  or  taut  (neutrally  buoyant)  moors  can  be  treated.  Given 
design  coefficients  for  the  mooring  line  material,  the  number  and _ (Con* t) 

20.  DISTRIBUTION/AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION  ' 

UNCLASS/UNLIMITED  X  SAME  AS  RPT.  DTIC _ _ _ 

22a.  NAK^  OF  RESPONSIBLE  INDIVIDUAL  22b.  TELEPHONE  22c.  OFFICE  SYMBOL  ; 

Jacgueline  B.  Riley _ 202-433-3881 _  • 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  ; 


DD  FORM  1473.  84MAR 


BLOCK  19  (Con’t) 

headings  of  the  mooring  lines,  and  specifications  of  design  constraints,  the 
program  selects  the  minimum  require  leg  diameter.  In  the  case  of  slack  moors 
the  anchor  locations  are  also  determined.  A  search  is  then  made  of  defined 
component  inventories  to  select  buoys,  anchors,  chains,  and  hawsers  as 
required  by  the  design  configuration.  The  program  utilized  the  equations  of 
static  equilibrium  at  a  reference  point  on  the  ship  and  the  ship  is  presumed 
to  be  a  rigid  body  with  zero  heave,  pitch  and  roll.  The  external  loadings  on 
the  ship  are  obtained  from  user-supplied  wind  and  current  load  tables  or 
approximate  analytical  functions.  Similarity  scaling  is  used  when  the  load 
tables  from  one  ship  (or  model)  are  used  to  approximate  the  loads  on  another. 
The  DESMOOR  program  was  developed  as  a  preliminary  design  tool  to  be  used  in 
conjunction  with  the  deep  sea  ship  mooring  static  and  dynamic  analysis 
features  of  the  SEADYN  computer  program. 


USER'S  MANUAL 


SEADYN/DSSM 

GENERAL  PURPOSE  CABLE  AND  DEEP  SEA 
SHIP  MOOR  ANALYSIS  COMPUTER  PROGRAM 

by 

R.  L.  WEBSTER 


NOVEMBER  1976 
REVISED  JULY  1978 
REVISED  SEPT  1978 
REVISED  OCT  1978 


Accession  F 

KTIS  GRA&I 
DTIC  TAB 
Uiiaiinounoed 
I  Justlficnti 

or 

% 

□ 

on. 

By 

Distribution/ 

Aval 

Diet 

M 

labim 

Avail' 

Spec 

'.y  Codes 

and/or 

lal 

PREPARED  FOR 


OCEAN  ENGINEERING  AND  CONSTRUCTION  PROJECT  OFFICE 
CHESAPEAKE  DIVISION 
NAVAL  FACILITIES  ENGINEERING  COMMAND 
WASHINGTON,  DC 


by 

ELECTRONICS  SYSTEM  DIVISION 
GENERAL  ELECTRIC  COMPANY 
SYRACUSE,  NY 


Revised  JULY  1978 

Revised  SEPT  1978 
Revised  OCT  1978  ' 


.  JWJyw^wwm«ww.w^^  'a.uwjjjji  n 

DSSM-2 


ERRATA 


The  July  1978  revision  to  this  document  has  changes  on  the  pages  listed 
below: 


V 

22 

24 

41 

50 

64 

69 

71 

77 

82 

83 

84 
89 

107 

130 

136 

142 

147 

157 

166 

168 

177 

178 
179. 


The  September  1978  revision  to  this  document  has  changes  on  the  pages 


listed  below: 

25 

50 

54 

63 

64 

69 

70 
84 
98 


The  October  1978  revision  to  this  document  has  changes 


This  page  added 
JULY  1978. 

Revised  SEPT  1978 

Revised  OCT  1978 


on  the  pages  listed  below: 


ABSTRACT 


The  basic  theory  and  Input  Instructions  for  a  general  purpose  nonlinear 
analysis  program  for  underwater  cable  and  truss  structures  Is  described.  A 
combination  of  the  finite  element  and  lumped  parameter  methods  Is  used.  Options 
are  provided  for  the  treatment  of  mooring  systems  for  surface  ships.  The 
SEADYN/DSSM  computer  program  treats  nonllnearltles  arising  from  large  displacements 
large  strains,  velocity-squared  drag  loading,  position  dependent  loading,  hyper¬ 
elastic  materials  and  constraints  on  the  surface  and  bottom  of  the  fluid  field. 

The  structure  may  be  partly  In  water  and  partly  In  air.  Analysis  options  provide 
for  iterative/incremental  solutions  for  static  and  transient  dynamic  responses. 

In  addition,  natural  frequencies  and  modes  may  be  determined  and  response 
spectrum  calculations  can  be  made  to  estimate  the  response  of  mooring  systems  to 
surface  waves  including  the  effects  of  steady  state  approximations  of  the  wave 
Induced  drift  forces.  A  unique  option  to  check  the  adequacy  of  lines,  buoys, 
and  anchors  is  also  provided. 


ADMINISTRATIVE  INFORMATION 


The  SEADYN/DSSM  computer  program  has  been  developed  over  the  past  three 
years  under  various  activities  supported  by  the  Electronics  Systems  Division 
of  the  General  Electric  Company,  the  United  States  Navy  and  the  private  efforts 
of  the  author  as  part  of  his  Doctoral  program  at  Cornell  University.  The  recent 
improvements  in  the  program,  the  generalization  to  treat  ship's  moorings  and 
the  production  of  this  manual  were  supported  by  the  Naval  Facilities  Engineering 
Command  through  contract  number  N62477-76-C-0002.  This  report  is  a  combined 
user's  manual  which  meets  the  requirements  of  Items  (10)  and  (13)  of  the 
Milestones  and  Deliverable  Schedule  dated  17  June  1975  as  revised  4  August  1975. 
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1.0 


INTRODUCTION 


This  report  provides  the  basic  Information  required  to  use  the  SEADYN/ 
DSSM  digital  computer  program.  The  analytical  capability  described  herein  Is 
considered  by  the  author  to  be  the  most  advanced  available  for  approximate 
analysis  of  cable  structures  and  mooring  systems  using  discrete  element  modeling 
techniques.  Advanced  state-of-the-art  methods  have  been  employed  In  all  of  the 
major  aspects  of  the  program.  Capability  Is  provided  for  the  static  and  dynamic 
analysis  of  cable  and  mooring  systems  In  which  all  of  the  major  nonl Inearl  ties 
have  been  given  consideration.  Emphasis  has  been  placed  on  versatility  and 
generality. 

The  SEADYN  program  has  developed  through  various  stages  since  the  summer 
of  1974.  The  starting  point  was  a  computer  program  known  as  NLIN[1]  which  was 
developed  by  the  Bechtel  Corporation  In  1973  under  contract  to  the  Electronic 
Systems  Division  of  the  General  Electric  Company,  Syracuse,  NY.  The  NLIN 
program  began  with  an  unpublished  cable  finite  element  program  written  by 
Or.  John  Leonard  of  the  Illinois  Institute  of  Technology  In  Chicago,  Ill.  The 
NLIN  effort  added  fluid  loading  effects  but  retained  the  linearized  Incremental 
solution  procedure.  This  author  was  technical  advisor  on  that  effort. 

Noting  that  various  deficiencies  existed  In  the  NLIN  formulation,  an 
effort  was  launched  to  study  the  general  underwater  structural  response  of  cable 
systems.  This  effort  eventually  led  to  the  doctoral  thesis  contained  In 
Reference  2  and  the  SEADYN  computer  program.  The  theoretical  work  and  program 
development  was  supported  by  the  author's  own  activities  and  by  various  work 
assignments  at  GE.  By  the  fall  of  1975  the  SEADYN  program  had  the  capability 
of  dealing  with  general  cable/truss  systems  using  Incremental,  Iterative  and 
generalized  Iterative/Incremental  solution  methods.  Both  the  static  and  time 
domain  dynamic  analysis  options  were  available.  The  discrete  equations  were 
derived  using  finite  element  theory,  and  the  program  had  a  single  element  In  Its 
library  (the  3-D  nonlinear  truss  element  otherwise  known  as  the  1-D  simplex 
element  In  3-D  space  [3]).  In  addition,  lumped  bodies  representing  buoys  and 
anchors  could  be  located  at  each  node  of  the  system.  The  program  was  written 
with  fixed  dimensions  which  limited  it  to  problems  with  no  more  than  3D  nodes 
and  IDO  elements. 


In  the  development  of  SEADYN  It  was  assumed  that  the  cable  system  could 
be  completely  submerged  1n  a  fluid  such  as  sea  water  and  that  imposed  loads  would 
Include  hydrodynamic  effects  from  the  relative  velocity  between  an  arbitrary 
flow  field  and  the  cable  motion.  Point  loads  could  also  be  Imposed  at  any  node. 
Constraints  were  provided  which  could  prevent  nodes  from  moving  out  of  the  region 
between  the  fluid  surface  and  the  bottom. 

In  the  fall  of  1975,  with  funding  from  the  Naval  Facilities  Engineering 
Command,  an  effort  was  begun  which  led  to  significant  generalization  of  the 
SEADYN  program.  The  main  thrust  of  this  effort  was  to  develop  an  effective  tool 
for  analyzing  the  dynamic  response  of  deep  sea  ship's  moors.  A  surface  ship 
and  surface  buoy  modeling  capability  was  added  to  SEADYN  which  allows  the  treatment 
of  static  response  of  complex  moors  to  winds,  surface  and  subsurface  currents, 
and  various  working  loads.  A  frequency  domain  dynamic  response  solution  option 
was  added  which  allows  estimates  of  the  regular  wave  and  spectral  responses  to 
surface  waves.  The  frequency  domain  solution  uses  ship  response  characteristics 
from  a  file  generated  external  to  the  program  (e.g.,  the  NSRDC  Ship  Motion 
Program  [4])  and  also  adjusts  the  static  reference  state  to  represent  the  effects 
of  the  steady-state  portion  of  the  second  order  wave-induced  drift  forces.  In 
addition,  the  dimension  restrictions  were  removed  from  the  number  of  nodes  and 
elements  and  various  improvements  were  made  In  the  solution  routines.  The  work 
also  developed  a  unique  checking  option  which  allows  the  user  to  evaluate  the 
adequacy  of  the  cable  or  mooring  system. 

This  report  gives  a  brief  summary  of  the  technical  developments  supporting 
the  SEADYN/DSSM  program,  describes  the  solution  procedures  used,  and  details  the 
input  requirements.  Following  a  description  of  the  program  output  and  error 
returns,  a  few  sample  problems  are  presented  to  demonstrate  the  versatility  of 
the  program.  Details  of  the  program  structure  are  given  in  the  prograniner's 
reference  manual  [5J  and  specific  applications  of  the  program  to  the  deep  sea 
moor  problem  are  given  in  Reference  6. 


2.0 


PROGRAM  DESCRIPTION 


2.1  OUTLINE  OF  THEORY 

2.1.1  Basic  Modeling  Assumptions 

The  approach  taken  In  the  SEADYN/DSSM  computer  program  to  model  cable 
and  mooring  systems  can  be  described  as  a  discrete  element  approach.  It  can 
be  considered  as  a  combination  of  the  finite  element  method  and  the  lumped 
parameter  method. 

In  Its  classical  form,  the  finite  element  method  seeks  to  represent 
continuous  physical  systems  with  a  set  of  discrete  or  finite  elements  which  are 
formulated  by  assuming  the  character  of  the  element  response  In  terms  of  a  set 
of  Interpolating  functions.  In  Its  usual  form  It  Is  equivalent  to  a  Galerkin 
form  of  the  method  of  weighted  residuals  where  the  weighting  functions  are 
defined  Individually  on  each  element.  Viewed  from  another  perspective,  the  finite 
element  method  Is  a  form  of  the  Raylelgh-Rltz  method  In  which  the  trial  functions 
are  defined  only  on  Individual  subregions  (elements)  of  the  system.  The  basic 
equations  for  cables,  mooring  lines  and  hawsers  are  obtained  using  a  simple  finite 
element  In  the  form  of  a  straight  line.  The  element  Is  assumed  to  be  straight 
both  before  and  after  deformation  of  the  system,  but  no  restriction  Is  placed  on 
the  amount  of  stretch  and/or  rotation  the  element  Is  subjected  to.  It  Is  further 
assumed  that  bending  and  torsional  effects  are  negligible.  In  the  case  of  bending 
this  means  that  the  bending  stiffness  of  the  cable  has  negligible  Influence  on 
the  global  response  of  the  system.  Neglect  of  the  torsional  effects  does  not 
mean  that  twist  Is  unimportant.  It  simply  means  that  coupling  between  twist  and 
extension  Is  assumed  to  have  little  effect  on  the  overall  shape  and  response  of 
the  system.  An  obvious  situation  where  this  is  an  invalid  assumption  Is  in  low 
tension  conditions  where  a  twist  instability  may  result  In  kinking  or  hockling. 

The  only  deformable  components  in  the  system  are  assumed  to  be  the 
cable  elements.  The  material  is  assumed  to  be  hyperelastic,  i.e.,  nonlinear, 
time- Independent  with  loading  and  unloading  curves  coincident.  With  the  exception 
of  the  frequency  domain  solutions,  internal  damping  effects  are  Ignored. 


Any  component  of  the  system  which  cannot  be  modeled  as  an  individual 
line  element  or  a  set  of  line  elements  is  assumed  to  be  a  rigid  body.  These 
rigid  components  are  assumed  to  be  lumped  at  a  single  point  in  the  system.  They 
may  be  assumed  to  have  only  a  point  effect  or  to  act  as  a  rigid  connector  for 
arbitrarily  placed  lines  (e.g.,  a  ship). 

The  system  may  be  totally  immersed  in  a  fluid,  suspended  between  two 
fluids  (e.g.,  water  and  air)  or  fluid  effects  may  be  ignored.  The  treatment  of 
fluid  effects  makes  the  fundamental  assumption  that  the  fluid  and  structure 
problems  are  uncoupled.  This  means  that  except  for  specific  localized  effects 
the  overall  fluid  field  characteristics  are  unaltered  by  the  presence  of  the 
structure.  Thus  such  things  as  flow  alteration  due  to  structural  movement  and 
blockage  effects  are  not  dealt  with.  More  specifics  on  the  assumptions  and 
limitations  of  the  fluid  interaction  with  the  structural  system  are  discussed  in 
Section  2.1.7  and  Reference  2. 

The  Lagrangian  approach  is  taken  in  describing  the  motion  of  the  system. 

In  this  approach  all  physical  variables  are  expressed  in  terms  of  their  values 
at  an  initial  reference  state.  It  is  possible  to  change  the  reference  state  by 
employing  generalized  coordinate  transformations  which  account  for  distortions 
and  rotations.  Analytical  procedures  which  begin  from  a  reference  state  and 
never  change  that  reference  are  called  total  Lagrangian.  Updated  Lagrangian 
is  an  obvious  title  for  methods  which  periodically  move  or  update  the  reference 
state.  Either  procedure  can  be  used  and  the  results  obtained  should  be  equivalent. 
In  the  developments  which  follow  the  configuration  of  a  system  (or  an  element) 

p 

is  designated  by  the  capital  letter  C  and  a  pre-superscript.  The  symbol  C  means 
the  reference  configuration  while  means  the  configuration  at  some  time,  t. 

The  definition  of  quantities  like  stress,  strain  and  displacement  usually  involve 
two  configurations.  A  pre-subscript  is  used  to  denote  the  reference  configuration 
for  such  cases.  Thus  means  a  quantity  in  measured  relative  to  °C. 

Details  of  the  finite  element  method  applicable  to  cable  systems  are 
given  in  Reference  2.  Only  brief  summaries  of  the  results  pertinent  to  the 
SEADYN/DSSM  program  are  given  here. 


A  finite  element  which  has  the  form  of  a  line  (i.e.,  one-dimensional) 
and  uses  only  the  field  parameters  at  the  two  ends  of  the  line  In  the  Interpolating 
functions  Is  referred  to  as  a  one-dimensional  simplex  element  [3].  When  the 
element  Is  a  straight  line  In  30  space  and  the  field  parameters  are  the  nodal 
(end  point)  displacements  the  element  Is  called  a  truss  element  In  structural 
terminology. 

Consider  a  single  straight  element  which  Is  defined  by  the  position  of 
two  nodes  (one  at  each  end).  Select  a  local  coordinate  system  with  the  x  axis 
extending  from  the  first  node  to  the  second.  The  other  two  axes  may  be  chosen 
arbitrary  under  the  restriction  that  they  form  a  right  handed  cartesian  reference 
frame.  When  the  element  Is  In  Its  unloaded  state  It  has  a  length  °L.  Assume 
the  material  constitutive  relation  has  the  form 


t^  _  tp  t 
0  0  0 


(2-1) 


where 

qS  Is  the  2nd  Piola-KIrchhoff  stress  In 
qC  Is  the  Green's  strain  In 

qE  Is  a  nonlinear  material  modulus  which  may  be  a  function 
of  strain. 

The  Incremental  form  of  this  constitutive  relation  can  be  written  for  a  small 
strain  Increment, 

0 


t+Atp  ,  tp  .  tp  ,t 
0  0  0  T  0 


(2-2) 


where  ^Ej  Is  the  tangent  modulus  evaluated  at  ^C. 


Making  the  finite  element  assumption  that  the  displacement  at  any 
position  along  the  element  Is  a  linear  function  of  the  displacements  of  the 
nodes  one  can  write 


wr1 tten 


where 


^[M]  Is  the  element  mass  matrix 
^{f}  represents  the  external  nodal  forces  In 
^{g}  represents  the  nodal  reactions  In 
^{R}  Is  called  the  force  residual 

The  mass  matrix  for  the  straight  element  can  be  written  In  two  forms: 


Consistent  Mass  Matrix 


[M] 


R  R.  R, 
0  A  L 


i  ^3 


1  r  , 

7  ^3| 


(2-7) 


Lumped  Mass  Matrix 


CM] 


R  R.  R, 
P  ^  L 


(2-8) 


R  R  R 

where  p  Is  the  element  material  density  In  C,  and  A  1s  the  element  cross- 

p 

sectional  area  In  C.  With  the  assumption  of  conversation  of  mass,  the  element 
mass  matrix  does  not  change  with  deformation.  The  pre-sub  and  superscripts 
are  used  in  equation  (2-6)  since  this  Is  not  true  of  fluid  added  mass. 


The  consistent  mass  matrix  is  obtained  from  the  kinetic  energy  and 
equation  (2-3).  The  lumped  form  can  be  obtained  by  the  intuitive  process  of 
lumping  half  of  the  element  mass  at  each  node  or  by  summing  all  the  terms  on 
each  row  of  the  consistent  mass  matrix  and  assigning  the  sum  to  the  diagonal 
position. 


The  external  forces  may  be  due  to  point  or  distributed  loads.  Point 
loads  appear  as  specific  entries  in  the  global  equations  (see  Section  2.1.7). 
Distributed  loads  are  usually  from  gravity  effects  and/or  fluid  loading.  Fluid 
loading  effects  are  discussed  in  Section  2.1.4.  The  general  form  of  the  gravity 
loading  is 


-{f}  *  /  "[N]'  ■'{y}  dx  (2-9) 

0 

Mhere 

P 

{y}  represents  the  components  of  the  element  specific  weight 
(1n  fluid)  relative  to  the  local  coordinate  system 

^[N]^  is  the  transpose  of  *^[N] 

Substitution  from  equation  (2-3)  Into  equation  (2-9)  and  noting  the 
orientation  of  the  element  with  respect  to  the  direction  of  gravity  leads  to  the 
conclusion  that  these  forces  are  equivalent  to  placing  one  half  of  the  element 
weight  acting  In  the  gravity  direction  at  each  node.  It  should  be  noted  that 
equation  (2-9)  assumes  mass  is  conserved. 

The  internal  forces  of  equation  (2-6)  can  be  written  [2] 


t/n;  =  tc  R.  '•L 


(2-10) 


where  {x}  Is  a  unit  vector  In  the  direction  of  the  deformed  element.  The 
stress  term,  ^S,  can  be  written 


tr  . 

tL 


(2-11) 


where  P  Is  the  element  load  In  the  deformed  state.  This  allows  the  force 
residual  to  be  written 


^{R}  =  ^{f}  -  V 


An  Incremental  form  of  the  motion  equations  can  be  written 


E[M]  {Aq}  =  {Af}  -  J[K^]  {Aq}  -  J[C]  (Aq) 


(2-12) 


(2-13) 


p[K-j.]  is  called  the  tangent  stiffness  matrix 


^[C]  is  an  incremental  damping  matrix 

The  stiffness  matrix  can  be  obtained  from  consideration  of  the  second  variation 
of  the  strain  energy  in  ^C.  The  result  is  [2] 


J[Kt]  - 


(2-15) 


where 
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= 


®2  " 


t, 

t 

R“2  ■ 

r"i 

t.. 

t 

r''2  ■ 

r''i 

t . 

t 

r'^2  ' 

In  many  situations  (e.g.,  fluid  loading)  the  force  is  dependent  on 
the  deflection.  In  this  case 


Uf}-^  Uq} 

®  ^  3  J{q} 


=  {a7}  +  J[Kp]  {Aq} 

The  incremental  motion  equations  are  then  written 

truT  r.l’i  ^  tr/.T  r.li  1_  tr 


(2-16) 


where 


r[M]  {Aq}  +  J[C]  {Aq}  +  J[K^]  {Aq}  *  (Af) 


r[Kj]  =  p[K^]  -  JCKj^] 


(2-17) 


(2-18) 


Equation  (2-13)  or  (2-18)  can  be  used  to  model  small  displacement 
response  about  a  steady  state  deformed  configuration,  or  it  can  be  used  in 
nonlinear  dynamics  by  recalculating  the  stiffness  matrix  at  each  step.  It  should 
be  noted  that  both  equations  neglect  the  position  dependent  effects  in  the  mass 
matrix.  The  incremental  load  rotation  matrix,  is  nonsymmetric  and  causes 

some  problems  in  applying  equation  (2-18).  Its  effect  when  small  increments  are 
used  is  felt  to  be  minimal  and  is  ignored  in  SEADYN/DSSM.  The  incremental  damping 
matrix  may  be  difficult  to  obtain  in  the  more  general  situations.  A  simplified 
treatment  is  discussed  in  Section  2.2.3. 


/V  . 


An  alternative  form  of  the  incremental  equations  is  obtained  from 
equation  (2-6)  by  expanding  only  the  internal  loads  in  a  Taylor  series  and 
neglecting  higher  order  terms.  The  result  Is 


^R{q}  +  rCKj]  {^q} 


(2-19) 


This  equation  is  linearized  by  approximating  and  with  their 

values  at  t+At  while  remaining  in  the  orientation  defined  by  C.  Any  damping 
effects  are  assumed,  to  be  included  in 

2.1.3  Lumped  Bodies 

Two  forms  of  lumped  bodies  are  considered.  The  simplest  form  treats 
the  body  as  a  single  point  with  three  displacement  degrees  oi  Freedom.  The 
point  is  assumed  to  have  mass  but  no  rotational  inertia.  No  elasticity  effects 
are  attributed  to  the  body,  but  it  may  be  a  means  of  inducing  fluid  loads  into 
the  system  (see  Section  2.1.4).  If  the  mass  of  the  body  is  m,  then  the  mass 
matrix  is  simply 


u 

This  body  affects  only  one  node  in  the  system. 


(2-20) 


The  second  form  is  that  of  a  rigid  body  with  spatial  dimensions.  The 
mass  is  still  assumed  to  be  concentrated  at  a  single  point  but  that  point  has 
six  degrees  of  freedom.  The  mass  matrix  assumed  for  this  case  has  specific 
forms  only  in  the  case  of  mooring  buoys  and  surface  ships.  The  program  generates 
the  mass  matrix  for  mooring  buoys  and  assumes  the  mass  matrix  is  defined  by  input 
from  the  ship's  motion  file  for  surface  ships.  See  Sections  2.1.8  and  2.1.9 
for  details.  In  either  case  :he  mass  matrix  is  a  6  x  6  matrix. 

2.1.4  Modeling  Fluid  Loads 

The  primary  assumption  regarding  the  effect  of  fluid  immersion  on  the 
cable  system  as  stated  in  Section  2.1.1  is  that  the  fluid  and  structural  problems 
are  uncoupled.  The  independence  principle  [7]  is  assumed  in  the  treatment  of 


* '  .'«•  .*•  .V- 


-12- 


fluid  loading  on  the  cables.  In  brief,  the  independence  principle  asserts 
that  the  fluid  loading  can  be  treated  as  resulting  from  two  separate  flows: 
one  normal  and  one  tangential  to  the  cable. 


The  fluid  induced  loads  can  be  separated  into  a  part  involving  the 
relative  velocity  between  the  cable  and  the  fluid  and  a  part  involving  the 
relative  acceleration.  The  velocity  related  terms  can  be  written 


(w)  =  Wf^  (nl  +  w.p  {X}  (2-21) 

where 

w^  =  1/2  D  (2-22) 

w.^  =  1/2  0^  C.^  D  (2-23) 


{nl.  {X}  arii  unit  vectors  in  the  directions  of  the 
normal  and  tangential  components  of  the 
relative  velocity 

Cn|,  Cj  are  normal  and  tangential  drag  coefficients 
D  is  the  drag  diameter  of  the  cable 
is  the  density  of  the  fluid 


The  fluid  loading  vector  for  an  element  is  then 


^{f}  =  /  ^[N]^  {w}  dL 


(2-24) 


0 

with  {w}  approximated  as  constant  over  the  length  this  reduces  to  the  usual 
lumping  procedure  of  placing  one  half  of  the  load  at  each  node. 


The  acceleration  related  portion  of  the  fluid  loading  can  be  separated 
into  a  part  due  to  flow  field  acceleration  and  a  part  due  to  structural  acceleration 
The  flow  field  acceleration  part  is  neglected  in  the  cable  loading  and  the 
structural  acceleration  part  is  treated  as  an  added  mass.  The  added  mass  matrix 
for  a  straight  element  can  be  written 


-13- 


[M 


added 


]  = 


A  L 


where 


(2-25) 


is  an  added  mass  coefficient. 

It  should  be  noted  that  there  is  no  added  mass  tangential  to  the  cable  which 
makes  this  a  position  dependent  term. 


In  some  of  the  solution  procedures  economies  can  be  achieved  if  the 
mass  matrix  is  diagonal.  A  lumped  form  of  the  combined  cable  and  fluid  added 
mass  matrix  can  be  written:  , 


-  R.  R,  C^,  p-  '^A  \ 

[M]  =  Cig]  -  ^ - 


=  [Mq]  + 


"J 


(2-26) 


where 


-  _  R 


+  C, 


M  '’f 


(2-27) 


This  form  still  presents  some  difficulties  which  will  be  dealt  with  in  the 
discussion  of  the  solution  procedures. 


Fluid  loads  on  submerged  lumped  bodies  are  estimated  using  an  approach 
similar  to  that  used  for  the  cable  element.  Two  forms  for  lumped  bodies  are 
considered:  spherical  and  an  end-faired  cylinder.  The  drag  loading  on  a  sphere 
is  given  by 

'''sphere’  =  ''2  »f 'd  “  (2-281 

where 


D  is  the  diameter  of  the  sphere 
Cq  is  the  drag  coefficient 

V  is  the  relative  velocity  between  the  fluid  and  the  point 
where  the  sphere  is  located. 

is  a  unit  vector  in  the  direction  of  relative  velocity. 


The  added  mass  for  a  sphere  is 

Cu  P  JT  If  D 

["added!  •  -fL-l -  tij]  (2-29) 

The  end  faired  cylinder  loading  and  added  mass  is  assumed  to  have  the  same 
form  as  that  for  a  cable  element  except  that  the  total  effect  is  placed  at  a 
single  point  rather  than  being  distribured  between  two  nodes. 

There  are  specific  fluid  effects  peculiar  to  surface  buoys  and  ships 
which  are  discussed  in  Sections  2.1.8  and  2.1.9. 

2.1.5  Coordinate  Systems  and  Transformations 

Equations  for  the  cable  element,  as  well  as  ships  and  mooring  buoys, 
are  most  readily  developed  in  a  local  or  intrinsic  coordinate  system  which  is 
considered  to  move  with  the  element.  The  development  of  a  global  set  of  equations 
which  represents  the  behavior  of  the  assembled  system  of  elements  requires  a 
single  global  coordinate  system.  Transformations  between  the  local  coordinate 
system  and  the  global  system  in  a  reference  configuration  are  accomplished 

by  the  usual  rotation  of  coordinates.  The  general  form  for  the  components  of  a 
vector  at  a  point  is 

R^^hocal  "  ^“^global 

Columns  of  the  transformation  matrix  are  the  components  of  a  unit  vector  in  the 
direction  of  the  local  coordinate  axis  expressed  in  the  global  system.  Since 
only  right-handed  cartesian  coordinates  are  considered,  the  inverse  transformation 
is  obtained  with  the  transpose  of  p[T].  All  contributions  from  the  external  and 
internal  loads,  etc.,  must  be  transformed  to  the  global  system  before  they  are 
combined  with  other  components  in  the  system. 


It  should  be  noted  that  these  coordinate  transformations  apply  only 
to  quantities  expressed  relative  to  a  specific  reference  configuration.  As  long 
as  the  reference  cc nfiguration  remains  fixed  the  individual  coordinate  transforma¬ 
tions  remain  unchanged  regardless  of  how  much  deformation  occurs  between  *^C  and 


Coordinate  transformations  involving  nodal  point  displacements  for  an 
element  can  be  written 

R^'^^local  "  "Vi"  R^^’^global  “  R^^^  R^'^^global 

rL  '  ij 

The  form  for  nodal  forces  is  similar.  When  one  transforms  the  incremental  motion 
equations  from  the  local  to  the  global  system  the  equations  take  the  form 

JCM]  ^^R^giobal  R^^  ^^^^global  *  R^T^  ^^‘’^global  *  ^^^^global 

(2-32) 

where 

Im  =  rCt]^  Jcm]  ^[t] 

-  rCt]^  Jcc]  rCT] 

JrKr]  ■  r[t]'  ‘[k,]  ,tT] 

“"global  ■  rCT]'  “"local 

Thus  it  is  seen  that  coordinate  rotations  do  not  alter  the  form  of  the 
equations.  Unless  specific  emphasis  is  required,  no  further  distinction  between 
the  local  and  global  equation  forms  will  be  made.  It  is  assumed  that  the  equations 
are  written  in  a  homogeneous  system,  i.e. ,  all  displacements,  forces,  stiffnesses, 
etc.,  are  in  the  same  coordinate  system. 


2.1.6  Slave/Master  Transformations 

The  generalized  rigid  bodies  used  to  model  ships  and  mooring  buoys 
require  some  special  manipulations  to  connect  them  to  the  rest  of  the  system. 
Their  motions  are  assumed  to  be  described  by  the  six  degrees  of  freedom  of  a 
single  point.  The  attachments,  however,  will  not  connect  to  the  rigid  body  at 
that  node.  Since  the  body  is  assumed  to  be  rigid,  it  is  possible  to  express  the 


motion  of  any  point  on  the  body  In  terms  of  the  motion  of  a  single  point  and 
the  relative  positions  on  the  body.  The  node  used  to  model  the  body  Is  called 
a  master  node.  Any  other  point  on  the  body  Is  called  a  slave  node.  Given  the 
six  components  of  motion  at  the  master  node,  the  translation  components  of  a 
slave  node  can  be  written: 

1  0  0  0  A2  -Ay  (q) 

{qlsigve  =  0  1  0  •  -^2  0  Ax  (2-23) 

b  0  °J  master 

OR 

^'’^slave  "  ^‘^^master 

where  ax.  Ay,  Az  are  components  of  the  distance  between  the  two  points  measured 
from  the  master  to  the  slave,  I.e.,  ax  =  etc.  The  matrix 

[Tjl^|]  can  be  viewed  as  a  generalized  form  of  the  coordinate  transformation 
represented  by  equation  (2-30)  and  the  transformation  procedures  of  the  previous 
section  can  be  employed. 

It  should  be  noted  that  the  slave/master  transformation  Involves  an 
alteration  of  the  number  of  degrees  of  freedom.  When  an  end  of  a  cable  element 
connects  to  a  slave  node,  the  application  of  equation  (2-33)  In  the  transformation 
Indicated  by  equation  (2-32)  results  In  a  stiffness  matrix  (etc.)  which  Is  9  x  9 
Instead  of  6  x  6.  It  should  further  be  noted  that  the  slave/master  transformation 
form  assumes  small  displacements  and  Is  therefore  only  applicable  to  the  Incremental 
equations. 

Calculations  of  the  position  of  a  slave  node  once  an  estimate  of  the 
position  of  the  master  node  Is  obtained  Is  more  appropriately  made  using  the 
general  coordinate  transformation  rather  than  employing  equation  (2-33).  The 
form  Is 

^^^slave  ”  ^’‘^master  (  ^’'^slave  "  ^^^master)  (2-34) 

where  b[T]  Is  a  rotation  matrix  of  the  form  (2-30)  which  represents  the  total 
^  R  t 

angle  changes  from  C  to  C  at  the  master  node. 


2.1.7  Global  Equation  Forms 

The  contributions  from  each  of  the  elements  in  the  system  (cables, 
lumped  bodies,  and  rigid  bodies)  can  be  combined  in  a  very  simple  and  direct 
manner  once  they  have  been  generated  and  transformed  to  the  global  coordinate 
system.  This  is  done  element  by  element  by  accumulating  the  element  contributions 
in  the  appropriate  position  of  the  global  arrays.  An  ordering  of  the  degrees  of 
freedom  is  implied  in  this  procedure.  The  order  assumed  in  SEADYN/DSSM  is  simply 
the  three  displacement  components  (x,y,z)  stored  in  the  order  of  the  node  number. 
Thus  the  global  nodal  displacement  vector  assumes  the  x  component  of  node  number 
one  is  first,  the  y  component  of  node  z  is  fourth,  etc.  By  requiring  fixed  (no 
movement)  nodes  and  slave  (movement  defined  in  terms  of  another  node)  nodes  to 
be  numbered  after  nodes  which  have  active  degrees  of  freedom,  the  solution  book¬ 
keeping  is  greatly  simplified. 

The  assembled  global  equations  have  essentially  the  same  form  as  the 
element  motion  equations.  The  main  distinction  is  that  the  order  of  the  equations 
is  increased  to  include  all  of  the  active  degrees  of  freedom  in  the  system.  Noting 
this  the  total  nonlinear  equations  of  motion  can  be  written 

J[M]  *  ^f)  -  ^g}  =■  ^R}  (2-35) 

The  two  incremental  forms  are 


J[M]  {Aq}  +  J  [C]  {Aq}  +  {Aq}  »  {a7} 

(2-36) 

t+At[„]  t+At^qj  ^  ^  -  ^^{g} 

(2-37) 

It  should  be  emphasized  that  these  represent  the  assembled  equations  for  the 
system  and  that  it  is  assumed  that  the  constraints  Implied  by  the  boundary 
conditions  and  slave/master  conditions  are  accounted  for.  The  dynamic  equations 
reduce  to  the  static  equations  when  the  time  dependent  terms  are  dropped.  Thus 
the  nonlinear  static  equation  is 


and  the  Incremental  static  equations  are 


{Aq}  «  {A?} 

(2-39) 

rCKj]  {Aq}*  -  *{g} 

(2-40) 

In  the  static  case  the  parameter  t  Is  used  to  signify  a  load  step  rather  than  a 
time  step.  The  static  equations  presume  a  stable  physical  system  has  been 
described  by  Imposing  adequate  constraints  on  the  system.  If  this  Is  not  so, 
the  stiffness  matrices  are  singular  and  there  Is  not  a  unique  configuration  of 
the  system  which  will  satisfy  the  equation.  In  the  case  of  nonlinear  systems 
(particularly  those  with  surface  ships  and  mooring  buoys)  the  static  global 
equations  may  be  Ill-conditioned.  This  means  they  are  nearly  singular  and 
numerical  errors  in  the  solution  procedure  may  lead  to  apparent  singularities. 

This  will  be  given  more  attention  In  Section  2.2.1. 

In  some  situations  one  or  more  of  the  motion  components  at  a  node  may 
be  constrained,  but  since  other  components  at  the  node  remain  active  or  the 
constraint  is  a  function  of  the  position  It  Is  not  appropriate  to  list  the  node 
as  a  fixed  one.  In  those  cases  an  artificial  stiffness  can  be  added  to  the  system 
which  has  the  effect  of  constraining  the  degree  of  freedom.  This  Is  done  by 
adding  a  large  value  to  the  diagonal  element  of  the  stiffness  matrix  at  the 
position  corresponding  to  the  constrained  degree  of  freedom. 

2.1.8  Surface  Ships 

The  rigid  body  element  is  used  to  model  surface  ships.  A  single  node 
point  Is  used  to  define  the  position  of  the  ship.  Since  the  node  must  express 
the  angular  as  well  as  spatial  position  of  the  ship  it  is  required  to  have  six 
degrees  of  freedom.  (The  SEADYN/DSSM  program  uses  two  consecutive  nodes  of  three 
degrees  of  freedom  each  to  define  a  ship.)  Attachments  of  mooring  lines  and/or 
working  lines  are  handled  through  the  slave/master  transformation. 


In  static  analyses  the  steady  state  effects  of  winds  and  surface 
currents  acting  on  the  ship  are  treated  as  lateral  and  longitudinal  forces  and 
a  yaw  moment  which  are  assumed  to  act  at  the  ship's  reference  point.  The  values 
of  these  forces  depend  on  the  flow  velocities  and  the  angle  between  the  ship's 
heading  and  the  flow  direction.  Empirical  load  tables  giving  load  coefficients 
versus  heading  or  analytical  load  functions  may  be  used.  The  empirical  approach 
is  given  in  NAVFAC  OM-26  [8]  and  both  approaches  were  sunnarized  in  Reference  9. 
That  presentation  is  repeated  in  Appendix  B  for  the  convenience  of  the  user. 

The  structure  of  a  data  file  for  ship  loads  is  described  in  Appendix  C. 

The  dynamic  equations  for  the  response  of  surface  ships  to  waves  are 
usually  given  in  an  incremental  linearized  form.  These  equations  have  the  form 

‘[Hj  t  (Jj)  *  JCC3]  (Uj)  ♦  ‘[K3]  tuj)  .  (fj)  (2-41) 

where 

{U5}  represents  the  six  components  of  ship's  motion 
^  (surge,  sway,  heave,  roll,  pitch,  yaw) 

^[Mj]  is  the  ships  mass  matrix  including  rotational  inertial  terms 
is  the  added  mass  due  to  fluid  acceleration  effects 

^[Cj]  is  an  equivalent  linearized  damping  matrix 

^[Kj]  is  the  ship's  hydrostatic  restoring  matrix 

{f.}  are  the  point  equivalent  forces  representing  the  wave 
induced  exciting  forces 

In  order  to  obtain  this  linearization  it  is  usually  assumed  that  the 
ship  is  driven  by  a  simple  harmonic  wave.  With  this  assumption,  the  forces, 
added  mass  and  damping  are  frequency  and  heading  dependent.  In  addition,  the 
linearization  of  the  roll  damping  term  makes  it  dependent  on  the  magnitude  of 
the  roll  angle.  Equations  of  the  form  of  (2-41)  can  be  obtained  for  slender 
bodies  using  strip  theory  [4].  A  more  general  theory  is  required  for  other 


forms  [10].  The  SEADYN/DSSM  program  assumes  the  values  for  these  coefficients 
are  provided  through  a  data  file.  The  format  of  that  data  file  Is  described  In 
Appendix  D. 


The  element  equation  represented  by  equation  (2-41)  can  be  manipulated 
as  any  other  element  equation  and  combined  with  the  global  equations  of  the 
system.  No  new  concepts  are  Involved  In  these  operations. 

2.1.9  Mooring  Buoys 

Ship's  moors  often  Involve  surface  buoys  which  support  the  mooring  line 
and  are  connected  to  the  ship  through  a  hawser.  This  type  of  buoy  usually  remains 
on  the  surface  where  1t  Is  subjected  to  the  effects  of  winds,  currents  and  waves. 
The  general  form  of  the  buoy  motion  equations  linearized  to  represent  small 
excursions  from  a  static  reference  state  can  be  written: 


tCMse]  (u, 


SB 


}  + 


SB 


}  +  JCK, 


Sb]  (Ugg)  =  {fgg} 


(2-42) 


where  the  various  terms  follow  the  previously  established  pattern.  The  forcing 
term  represented  by  {f^g}  deals  only  with  wave  excitation.  Static  load  effects 
on  a  mooring  buoy  follow  the  same  form  as  described  previously  In  Section  2.1.3. 


In  order  to  avoid  a  great  deal  of  complexity,  it  Is  assumed  that  the 
buoy  is  spherical  in  shape  and  that  a  local  cartesian  coordinate  system  is 
selected  which  Is  vertical  in  the  z  direction  and  has  the  incident  wave  traveling 
1n  the  +x  direction.  No  loss  of  generality  is  incurred  with  this  choice  of 
coordinate  system  on  a  spherical  buoy  since  the  character  of  the  coefficients  do 
not  depend  on  the  orientation  or  attitude  of  the  buoy.  The  problem  is  further 
simplified  if  it  Is  assumed  that  the  buoy  Is  homogeneous  with  the  center  of 
gravity  at  the  geometric  center  of  the  buoy  and  that  the  geometric  center  is 
located  at  the  water  line  in  the  reference  state. 


Attachements  of  hawser  and  mooring  lines  to  the  buoy  can  be  readily 
handled  if  their  positions  relative  to  the  local  coordinate  system  are  known.  The 
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rigid  link  transformation  of  equation  (2-33)  is  used  for  this  purpose.  The 
positions  of  the  attachments  can  be  found  from  the  static  solution. 

Since  it  is  assumed  that  the  equations  use  as  the  reference  configuration, 
the  configuration  notation  will  be  dropped  for  this  discussion. 


Transformation  from  the  local  to  the  global  system  for  assembly  of  the 
buoy  equations  with  the  rest  of  the  moor  system  equations  is  a  straightforward 
process  which  follows  the  method  outlined  previously.  For  this  reason  only  the 
coefficients  in  the  local  coordinate  system  will  be  given  here.  It  should  be 
kept  in  mind  that  the  following  equations  represent  the  incremental  motion  equations 
for  a  surface  buoy  in  the  local  coordinate  system  just  described. 


Given  that  the  buoy  has  a  mass  designated  by  m  and  a  mass  moment  of 
inertia,  the  buoy  portion  of  the  mass  matrix  is 


f, 


m  I. 


[Mg] 


1 


^3' 


(2-43) 


where  I3  is  the  identity  matrix  of  order  3.  The  assumption  of  a  homogeneous 
sphere  should  be  recalled  at  this  point.  If  the  attachments  contribute  significant 
mass,  their  effects  can  be  treated  by  including  additional  lumped  masses  at 
those  nodes. 


The  added  mass  has  the  form: 


■VI 


XX 


[M,]  •  i 


(2-44)  ;  f': 
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The  mass  matrix  is  then  given  by 


The  wave  damping  matrix  has  a  form  similar  to  the  added  mass  matrix. 
Specific  values  for  the  added  mass  and  damping  coefficients  for  a  sphere  were 
given  by  Patton  [11].  His  values  were  obtained  by  curve  fitting  the  analytical 
results  presented  by  Kim  [12].  The  nondimensional  values  obtained  were  * 


1.089  +  0.529  a'  for  0  <  a'  <  0.74 

l/(-0.0318  +  0.954  a')  for  0.74  <  a*  <  3.4 

1.85  for  0  <  a'  <0.1 


1.02  a 


-0.256 


0 


-0.069  +  0.715  a' 
1.595 

0.126  +  1.7  a' 
1.18 


0.1  <  a*  <  3.4 
for  0  <  a'  <0.1 
for  0.1  <  a'  <  1.37 
for  1.37  <  a'  <  3.4 
for  0  <  a'  <  0.4 
for  0.4  <  a'  <  3.4 


where  a'  is  the  nondimensional  frequency  given  by 

I  2TTa  a 
'  ^  "  9 


a  =  radius  of  sphere 
X  «  wavelength 

ijj  =  circular  frequency  of  wave  encounter 


(2-46) 

(2-47) 

(2-48) 

(2-49) 

(2-50) 


damping 


The  added  mass  terms  are  nondimensionalized  by  the  factor  pa^  and  the 
terms  by  the  factor  pa  u,  thus 


*  See  the  Addendum  at  the  end  of  this  report  for  an  alternative  curve  fitting  which 
was  used  in  the  program. 
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The  roll,  pitch  and  yaw  added  mass  terms  arise  from  fluid  viscosity  and 
they  can  be  written 


A  3  A  3  A  3  P^~ 

''ee  ''h  3 


where 


f2^, 


1  +  6a 

- 

1  +  26a  +  26  a 


(2-53) 


(2-54) 


and  V  =  the  kinematic  viscosity  of  the  fluid.  Since  the  kinematic  viscosity 

-5  2 

of  water  is  of  the  order  of  10  ft  /sec,  the  rotational  added  mass  terms  will 
be  small  compared  to  the  buoy  inertia  terms.  Damping  due  to  rotational  motion 
is  very  small  and  will  be  neglected.  Thus, 


^90  “ 


(2-55) 


When  the  center  of  pressure  does  not  coincide  with  the  center  of  gravity 
of  the  buoy,  a  coupling  between  lateral  and  rotational  motions  exists.  These 
terms  can  be  written: 


With  the  origin  of  the  local  coordinate  system  at  the  geometric  center 

(also  center  of  gravity),  z  is  zero.  The  center  of  pressure  for  a  half 

Cy 

submerged  sphere  is  obtained  from 


The  damping  terms  presented  above  do  not  represent  the  effects  of  viscous 
drag.  The  viscous  terms  involve  the  square  of  the  relative  velocity  between  the 
buoy  and  the  fluid  and  are  therefore  nonlinear.  The  viscous  effects  are 
generally  of  less  importance  than  the  wave  damping.  Obviously,  this  is  not  the 
case  for  the  rotational  movement  since  those  terms  are  zero  for  wave  damping.  In 
a  free-floating  buoy  the  viscous  rotational  terms  would  play  an  important  part, 
but  in  a  mooring  system  where  the  hawser  and  mooring  leg  restrain  the  buoy  the 
rotation  is  limited.  The»«fore,  all  of  the  viscous  terms  will  be  neglected 
rather  than  attempting  to  linearize  them. 

The  only  nonzero  hydrostatic  restoring  force  on  a  half-submerged  spherical 
buoy  acts  in  the  heave  direction.  Its  value  for  small  displacements  is 

<22  =  09  (2-59) 

When  the  buoy  provides  a  connection  between  a  mooring  line  and  a  hawser  it 
develops  an  additional  stiffness  (resistance  to  motion)  due  to  the  tensile  force 
being  transmitted  across  it.  This  is  analogous  to  the  geometric  stiffness  term, 
[Kg],  seen  in  the  cable  element  stiffness  matrix  in  Section  2.1.2.  This  geometric 
stiffness  effect  is  automatically  taken  into  account  by  the  rigid  link  transforma¬ 
tion  of  the  elements  representing  the  attached  lines. 

The  right-hand  side  of  equation  (2-42)  represents  the  forces  due  to 
surface  waves.  Assuming  the  wave  is  harmonic  in  form,  Kim  [13]  shows  that  the 
wave  'nduced  forces  can  be  written 


where 


'SB' 


XX 

w 

0 

XX 

w 

A^ 

Z 

+ 

Z  + 

zz 

w 

zz 

w 

’60  w 


ftS.  -  ,3  mS 

_  3  .,S 

‘^xx  -  “  "x 


a!«  6.  +  cr,  e .+ i„  e 


66  W  2  W 


It  =  P9  /  6^  ri  ds  5  pg  /  cos  (a'x)  n,  ds 


lo  =  pg  /  e' 


,a'  (z+ix) 


-  zijj)  ds 


(2-60) 
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(2-62) 

(2-63) 

(2-64) 


m 


For  the  half  submerged  sphere 


=  I  =0 

^06  So  ^2  ^ 


(2-65) 


The  wave  pressure  component  of  the  heave  exciting  force,  I^,  is  closely 
approximated  by  the  pressure  at  the  water  surface  distributed  over  the  cross- 
section  at  the  water  surface.  A  plot  of  this  function  versus  the  nondimensional 
frequency  is  given  in  Figure  2-1. 

S  S  S  S 

The  values  for  and  for  a  half  submerged  sphere  were  given 

by  Kim  [13].  A  polynomial  curve  fit  of  those  functions  plus  the  curve  for  are 
summarized  in  the  following  table. 
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SUMMARY  OF  WAVE  EXCITING  FORCE  COEFFICIENTS 


FUNCTION 

^0 

“l 

^2 

^3 

^4 

^5 

®6 

^7 

X 

0. 

1 . 7586 

-8.2171 

12.0253 

-8.2882 

3.0081 

-0.5576 

0.0416 

1 . 7868 

1.0552 

-10.5792 

17.2729 

-13.2596 

5.3272 

-1.0781 

0.0866 

"x 

1.0833 

0.0833 

-1.4496 

0.7753 

-0.0705 

-0.0314 

0.0055 

0. 

0. 

4.2382 

-8.5367 

8.2743 

-4.6849 

1.5633 

-0.2823 

0.0211 

2 

I-j/Tra  og 

1.0 

-0.0004 

-0.1218 

-0.0026 

0.0069 

-0.0007 

0. 

0. 

where  the  function  Toriii  Is: 


N 

=  i-=^0 


(2-66) 


Thus  all  of  the  terms  necessary  for  treating  the  small  displacement 
behavior  of  a  restrained,  ha If- submerged  spherical  buoy  are  available.  The 
assumptions  employed  appear  to  be  reasonable  and  should  lead  to  a  good  approxima¬ 
tion  of  the  buoy  effects  in  a  deep  sea  moor.  Although  a  spherical  buoy  has  been 
assumed,  a  comparison  of  the  curves  presented  by  Kim  [12]  for  a  sphere  and  those 
presented  by  Garrison  [10]  for  a  ha If- submerged  cylindrical  buoy  with  an  aspect 
ratio  of  1.0  shows  that  the  added  mass  and  damping  coefficients  are  quite 
similar.  Therefore,  it  is  reasonable  to  expect  the  sphere  equations  to  give  at 
least  an  order-of-magnitude  approximation  of  a  cylindrical  buoy. 


The  small  displacement  assumption  deserves  some  further  comnent.  For 
wavelengths  of  the  order  of  the  buoy  diameter  one  would  expect  the  buoy  motion 
to  be  small.  However,  as  the  wavelep.wth  increases  the  buoy  motions  increase. 

Since  most  of  the  wave  energy  is  expected  in  the  longer  wavelengths,  the  buoy 
could  be  expected  to  see  large  motions  which  would  cause  these  equations  to  be 
inaccurate.  When  no  ship  is  in  the  system,  this  inaccuracy  could  be  serious. 
Fortunately,  the  ship  motion  becomes  a  significant  effect  in  the  longer  wavelengths 


and  the  buoy  motion  is  dominated  by  the  ship  movement  as  transmitted  through 
the  hawser  and  reacted  by  the  mooring  line.  In  this  case  the  contributions 
from  the  buoy  itself  (though  in  error)  would  generally  be  insignificant.  For 
the  shorter  wavelengths,  the  ship  appears  nearly  fixed  and  the  exciting  forces 
due  to  wave  action  on  the  buoy  become  important.  It  is  fortuitous  that  this  is 
the  range  in  which  the  buoy  equations  are  most  accurate. 

2.1.10  Deep  Sea  Ship's  Moors 

The  basic  components  required  in  modeling  deep  sea  mooring  systems 
for  surface  ships  are: 

a)  surface  ship 

b)  mooring  lines  (usually  submerged) 

c)  mooring  buoys 

d)  hawser  (usually  in  air  and  subjected  to  wind  loads) 

e)  floats,  sinkers  and  anchors 

Each  of  these  components  have  been  dealt  with  to  some  extent  in  the  previous 
discussion.  Unfortunately  it  is  not  possible  to  develop  a  fully  general  non¬ 
linear  analysis  of  ship's  mooring  dynamics.  The  primary  reason  for  this  is  the 
highly  complex  nature  of  the  interaction  between  the  sea  and  surface  bodies  such 
as  ships  and  mooring  buoys.  The  equations  presented  for  these  bodies  are 
linearized  equations  which  address  only  the  response  to  harmonic,  long-crested 
waves. 


The  theoretical  approach  used  in  dealing  with  mooring  systems  follows 
through  a  series  of  approximations.  The  first  step  is  to  obtain  a  description 
of  the  mooring  system  in  the  quiescent  state  where  only  gravity  loads  are  involved. 
This  is  called  a  dead  load  analysis.  There  are  some  pitfalls  in  this  step  which 
are  related  to  the  geometric  nonlinearity  of  the  problem.  One  does  not  usually 
know,  a  priori,  what  the  dead  load  configuration  of  the  system  is  and  the  terms 
in  the  static  equations  are  configuration  dependent.  This  is  the  so-called 
initial  configuration  problem  and  it  is  dealt  with  in  some  detail  in  Reference  2. 
Sections  2.2.7  and  3.5  also  discuss  the  problem. 
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The  active  loads  on  the  mooring  system  are  assumed  to  be  winds,  surface 
and  subsurface  currents  and  surface  waves.  The  effect  of  surface  waves  is 
primarily  a  dynamic  phenomenon  while  the  static  or  steady-state  effects  of  winds 
and  currents  predominate.  Therefore,  it  is  assumed  that  the  active  loads  on  the 
system  can  be  separated  into  a  static  effect,  which  is  primarily  due  to  winds 
and  currents,  and  a  dynamic  effect  from  surface  waves. 

The  next  step  in  the  analysis  is  then  a  static  analysis  to  obtain  the 
response  to  winds  and  currents.  This  is  referred  to  as  a  live  load  analysis  and 
it  may  include  point  loads  representative  of  imposed  work  loads.  Nonlinearities 
in  the  system  also  play  an  important  part  in  a  live  load  analysis.  The  geometric 
nonlinearity  is  still  present.  A  significant  new  nonlinearity  comes  from  the 
nonconservative,  position  dependent  loads.  The  flow  induced  loads  are  strongly 
dependent  on  the  orientation  of  the  various  elements  and  the  loads  change  direction 
and  magnitude  as  the  system  changes  shape  and  position. 

Material  nonlinearities  also  play  an  important  part  at  this  stage.  The 
most  pronounced  effect  comes  from  the  fact  that  the  lines  cannot  support  compressive 
loads.  Should  the  imposed  loads  cause  any  of  the  legs  to  go  slack  (usually  in 
taut  moors  with  neutrally  buoyant  lines)  the  material  stiffness  goes  to  zero  and 
the  group  of  elements  involved  with  the  leg  are  part  of  an  unstable  structure 
and  the  stiffness  matrix  becomes  singular. 

An  important  feature  in  slack  moors  (negatively  buoyant  legs)  is  their 
ability  to  resist  loads  by  changing  shape  and  by  lifting  line  or  laying  it  down 
at  the  bottom.  Modeling  this  interaction  with  the  bottom  using  reasonably  long 
cable  elements  produces  some  approximations  that  must  be  kept  in  mind  while  inter¬ 
preting  the  results. 

In  certain  situations  it  is  not  possible  (or  feasible)  to  obtain  a 
dead  load  configuration  and  then  proceed  to  the  live  load  analysis.  Direct 
solution  for  the  combined  dead  and  live  loads  is  usually  required  in  these  cases. 

One  example  is  the  solution  for  a  single  point  moor.  In  this  case  the  quiescent 
state  is  of  little  value  and  solution  for  the  combined  effects  is  usually  quite 


easy  to  obtain  if  one  has  a  good  estimate  of  the  total  horizontal  load  to  be 
supported  by  the  moor.  Another  situation  where  this  procedure  may  be  required 
Is  In  dealing  with  taut  moors  where  legs  gc:  slack.  In  this  case  the  slack  legs 
could  be  Ignored  (not  Included  In  the  model)  and  the  dead  load  Included  In  the 
live  load  analysis. 

The  analysis  of  wave  Induced  dynamics  begins  with  the  static  configura¬ 
tion  developed  by  the  wind  and  currents.  Of  necessity  this  is  a  frequency  domain 
analysis.  The  equations  presented  for  the  ship  and  mooring  buoy  dynamics  were 
obtained  by  assuming  the  excitation  was  from  a  harmonic  wave.  The  linearization 
process  renders  the  equations  frequency  dependent  and  limited  to  small  motion 
amplitudes. 

The  incremental  equation  (2-36)  with  the  load  rotation  term  neglected 
and  the  reference  state  taken  to  be  (where  is  the  static  reference  state) 
can  be  used  to  represent  the  motion  in  the  frequency  domain.  In  this  case  the 
only  nonzero  elements  of  {aT}  are  those  contributed  from  the  surface  buoys  and 
ship.  Assuming  these  forces  are  harmonic,  equation  (2-36)  can  be  written 

[M]  {ii}  +  [C]  {u)  +  [K.^]  {u}  =  {F}  e^“^  (2-67) 

Here  the  reference  indices  have  been  dropped  and  it  1s  understood  that  (u} 
represents  the  incremental  displacement,  {Aq},  from  ^C.  The  vector  {F>  is  a 
set  of  complex  force  amplitudes.  Assuming  the  equation  is  linear  and  transient 
effects  damp  out  the  response  has  the  form 

{u}  =  {U}  e’“^  (2-68) 

where  fU}  is  a  set  of  complex  response  amplitudes  and  u  Is  the  circular  frequency 
of  the  driving  function.  Substitution  of  equation  (2-68)  into  (2-67)  yields 


Solution  of  this  set  of  complex  simultaneous  linear  algebraic  equations  gives 
the  dynamic  response  amplitudes.  The  magnitude  of  the  response  in  each  degree 
of  freedom  is  given  by 


|U,.  I  =  U.  U>  (2-70) 

where  LI.*  is  the  complex  conjugate  of  the  ith  component  of  {U}.  The  phase  shift 
between  the  incident  wave  and  the  response  is  given  by 


(2-71) 


A  phase  shift  of  zero  corresponds  to  the  response  being  in  phase  with  the  incident 
wave  at  the  ship's  reference  point. 


The  dynamic  tension  response  in  each  element  is  obtained  by  adding  the 
displacement  increments  to  the  nodal  positions,  recalculating  the  element  tensions 
from  the  constitutive  relations  using  the  new  lengths,  and  subtracting  the  tensions 
in  the  static  reference  state. 


2.2  PROGRAM  IMPLEMENTATION 

2.2.1  Static  Solution  Procedures 

The  SEADYN/DSSM  program  offers  three  basic  solution  methods  for  static 
analyses.  The  various  features  of  each  of  them  are  discussed  at  length  in 
Reference  2  so  only  a  brief  description  will  be  given  here. 

The  first  method  is  a  sequence  of  linear  increment  (SLI  method).  The 
loads  are  divided  into  a  sequence  of  increments  and  the  basic  incremental 
equation  (2-39)  is  repeatedly  applied.  The  SLI  method  requires  the  regeneration 
of  the  incremental  stiffness  matrix  at  each  step  to  reflect  the  changes  in 
position  and  constitutive  relations.  It  has  the  undesirable  feature  of  drifting 
from  the  correct  solution  through  accumulating  errors  and  small  increments  are 
required  for  accuracy. 
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The  second  method  works  with  the  incremental  equation  (2-40).  The 
first  step  is  identical  to  the  SLI  first  step,  but  on  each  succeeding  step  the 
force  residual  from  the  previous  step  is  fed  back  as  a  corrector.  For  this 
reason  it  is  called  the  residual  feedback  method  (RFB).  Although  it  is  a 
non-iterative  method  it  tends  to  be  self-correcting.  On  responses  which  are 
reasonably  well  behaved  (particularly  monotonic  responses)  the  RFB  method  gives 
accurate  results  with  significantly  fewer  steps  than  does  the  SLI  method.  The 
RFB  method  costs  somewhat  more  per  step  because  of  the  residual  calculation. 

The  recalculation  of  the  incremental  stiffness  matrix  is  still  required  at  each 
step. 


A  more  general  method  which  employs  iterations  to  solve  equation  (2-38) 
is  known  as  the  modified  Newton-Raphson  method  (MNR).  It  begins  with  an  estimated 
configuration  and  uses  an  estimate  of  the  tangent  stiffness  matrix  to  obtain 
successive  displacement  increments  which  hopefully  lead  to  a  zero  force  residual. 
Being  an  iterative  method  it  is  by  far  the  most  accurate  of  the  three  methods. 

This  accuracy  is  not  free,  however.  In  the  first  place  the  method  is  not 
unconditionally  convergent  and  the  size  of  the  load  step  required  to  assure  con¬ 
vergence  is  not  easily  determined.  In  some  cases  it  may  be  extremely  small.  In 
ill-conditioned  systems  it  may  not  be  possible  to  get  convergence  without  some 
special  auxiliary  procedures. 


The  general  form  of  the  MNR  method  is 


[1^ 


(2-72) 


where  the  superscript  i  refers  to  the  iteration  step.  The  matrix  is  referred 
to  as  an  estimator  matrix.  If  [K]  is  the  incremental  matrix  [Kj],  and  it  is 
recalculated  at  each  iteration,  the  usual  Newton-Raphson  procedure  is  obtained. 

If  [K]  is  an  approximation  of  and/or  it  is  not  recalculated  at  each  iteration 
one  has  a  modified  Newton-Raphson  method.  When  [1<]  is  not  changed  and  the 
response  curve  is  monotonically  increasing  function  the  successive  estimates  will 


usually  oscillate  about  the  correct  solution.  If  the  step  size  Is  small  enough 
the  estimates  will  tend  towards  the  solution  (converge).  If  the  step  size  Is 
too  large  the  estimates  will  diverge. 

Once  oscillating  estimates  are  detected,  various  accelerating  procedures 
can  be  employed,  and  In  some  Instances  they  will  work  even  when  the  oscillations 
are  divergent.  The  program  detects  oscillations  by  monitoring  the  degree  of 
freedom  with  the  largest  Initial  response  (I.e.,  largest  component  of 
Oscillation  Is  signalled  by  repeated  sign  changes  of  the  critical  {Aq}  component. 
The  simplest  acceleration  scheme  averages  the  two  alternating  estimates  using 
the  sizes  of  the  critical  components  to  weight  the  average.  An  optional  procedure 
uses  a  one-dimensional  search  to  seek  a  close  estimate  of  {Aq}  which  crosses  the 
correct  solution.  This  search  Is  initiated  when  the  1th  Iteration  signals  a 
large  oscillation.  The  search  begins  at  the  (1-1)^^  position  and  takes  small 
steps  In  the  direction  of  the  1^^  1ncreit«nt.  A  new  step  Is  tried  at  each  position 
until  the  new  increment  reverses  direction.  At  this  point  the  last  two  positions 
are  averaged  to  get  a  new  starting  estimate  for  the  iteration.  The  size  of  the 
step  used  In  this  1-0  search  Is  controlled  by  Input.  The  Input  factor  is  the 
fraction  of  the  1^^  step  which  is  taken  as  the  first  step  of  the  search. 

Various  options  are  provided  to  measure  convergence  of  the  iterations. 
Reference  2  should  be  consulted  along  with  Section  3. 2. 3. 2  for  discussions  of  the 
convergence  criteria. 

Various  load  Incrementing  options  are  provided  with  the  static  analyses. 
In  general,  the  gravity  load  will  be  increased  in  increments  in  a  dead  load 
analysis.  Fluid  induced  loads  are  increased  in  increments  during  live  load 
analyses.  Both  analyses  allow  point  loads  which  are  incremented  from  zero  to 
the  specified  value,  held  constant,  or  reduced  from  the  value  to  zero.  It  is 
usually  assumed  that  gravity  loads  are  held  fixed  during  a  live  load  analysis,  but 
an  option  Is  provided  which  allows  gravity  load  to  be  incremented  during  a  live 
load  analysis. 
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When  ships  are  used  In  a  static  analysis  the  ship  is  constrained  to 
remain  on  the  surface.  This  normally  means  the  ship  is  fixed  in  heave,  roll 
and  pitch.  It  is  possible  to  input  stiffness  terms  to  allow  small  heave  and 
angular  motions  out  of  the  surface  plane.  No  checks  are  made  on  the  magnitude 
of  those  responses,  so  the  user  is  warned  to  evaluate  the  validity  of  his  results 
when  these  stiffnesses  are  used. 

2.2.2  Time  Domain  Dynamic  Solutions 

The  SEADYN/DSSM  program  provides  for  four  different  time  domain  solution 
methods.  Each  of  them  utilize  the  following  set  of  forward  difference  equations 
to  develop  the  solutions: 

=  ^{q}  +  At  ^{q)  +  ^  ^{q‘}  +  a  At  -  '‘^{ql) 

+  BAt^  -  ^q})  (2-73) 

=  ^{q}  +  At  ^q}  +  yAt  -  ^{q}) 

where  At  is  the  time  step  and  a,  B,  and  y  are  integration  parameters.  The 
usual  Newmark  formulas  are  obtained  with  a  =  0. 

When  equations  (2-73)  are  augmented  with  a  motion  equation  for  the 
system,  one  has  three  simultaneous  differential  equations  with  the  nodal  displace¬ 
ments,  velocities  and  accelerations  at  as  unknowns.  Specific  forms  of  the 
solutions  routines  for  this  set  of  equations  are  developed  in  Reference  2.  Only 
the  final  forms  will  be  presented  here. 

Direct  Iteration  Method  (DI) 

This  method  uses  equations  (2-73),  (2-35)  and  (2-26)  to  obtain  an 
iterative  formulation.  The  form  is 


'A'  “.■V  -I*  '•;.*r  •-’  V*  -L'  -L*  - '  -  r^r ; 
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'*»'«,)("*’)  =  '*«(,)(")  »  („,6)  At^  4(‘*''*(a5)<"*'>) 

t+At^"j{n+l)  (2-74) 


k» 


The  iteration  can  be  started  with  equations  (2-73)  with  o»0  =  Y“Oor  with  the 
residual  feedback  solution  to  be  described  below.  The  DI  method  has  been  shown 
to  be  an  accurate  and  cost  effective  approach.  Since  it  is  an  iterative  method 
which  retains  all  of  the  nonlinear  terms  it  readily  deals  with  all  of  the 
important  nonlinear  effects  such  as  position  dependent  loads  and  constraints  slack 
segments,  and  changing  geometry.  Of  particular  interest  is  the  ease  with  which 
problems  with  defined  motion  are  dealt  with.  When  problems  require  the  imposition 
of  a  known  motion  (e.g.,  cable  towing  problems)  the  moved  node  is  simply  held 
at  the  required  dynamic  state  during  the  iterations. 

Sequence  of  Linear  Increments  (SLI) 

This  method  starts  from  equations  (2-73)  and  uses  the  Incremental  form 
(2-36).  An  iterative  solution  is  also  employed,  but  it  has  been  shown  that  this 
method  costs  essentially  the  same  as  the  DI  method  and  is  much  less  accurate  [2]. 
It,  too,  is  capable  of  solving  the  moving  boundary  problem  but  an  augmented  global 
stiffness  matrix  involving  the  moved  degrees  of  freedom  is  required.  This 
increases  the  storage  requirements,  reduces  the  speed  and  increases  the  numerical 
error  potential. 

Residual  Feedback  Method  (RFB) 

Equations  (2-73)  are  Inverted  and  substituted  into  the  incremental 
form  (2-37).  The  result  is 


[K,„]  (AQ)  =  !%„) 


(2-75) 


Equation  (2-75)  Is  a  linear  algebraic  equation  which  can  be  solved  by  the  usual 
manner.  No  iteration  is  required  and  the  method  has  a  self-correcting  feature 
similar  to  the  static  RFB  method.  This  method  follows  more  closely  the  traditional 
Newmark  implicit  integration  form  [14].  It  has  the  unconditional  stability 
features  that  have  made  the  method  so  popular.  No  provision  is  made  for  solving 
the  moving  boundary  problem,  however.  The  fact  that  the  RFB  method  requires 
the  formation  of  a  global  stiffness  matrix,  the  evaluation  of  a  residual  and  the 
solution  of  simultaneous  equations  reduces  its  cost  effectiveness.  This  is 
usually  compensated  for  by  using  larger  time  steps. 


Modified  Newton-Raphson  Form  (MNR) 

This  is  an  iterative  method  based  on  a  procedure  similar  to  the  MNR 
static  solution.  The  form  is 

[I4]  a  ("“(Aq)'"*’)  •  '*“(«<">  -  (2-76) 
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The  iteration  is  started  with 


t+At^j^j{0) 
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where  ^''’^^{Aq}^^^  is  obtained  from  an  RFB  solution.  As  an  alternative  to 
equations  (2-77)  the  iterations  can  be  started  from  equations  (2-73)  with 

Q  =  3  =  Y  =  0. 


The  accuracy  of  this  method  has  been  demonstrated  but  its  cost  per  time 
step  does  not  appear  competitive  with  the  DI  method  [2].  Moving  boundary 
solutions  have  not  been  implemented  with  this  solution. 

It  would  appear  from  the  studies  done  in  Reference  2  that  the  best 
general  purpose  solution  method  of  the  four  is  the  DI  method.  On  some  problems 
the  method  may  encounter  some  difficulties,  however.  The  iterations  are  convergent 
only  for  step  sizes  of  the  order  it  takes  an  axial  wcve  to  traverse  the  shortest, 
stiffest  element.  When  the  imposed  forces  are  slowly  varying  the  routine  will 
attempt  to  increase  the  step  size  beyond  stable  behavior  limits  because  the 
motions  are  small  in  small  time  steps  and  the  iterations  converge  rapidly.  Control 
are  provided  to  the  user  to  prevent  step  size  increases  in  this  situation. 

In  some  situations  the  RFB  solution  may  be  found  very  cost  effective. 
Generally  this  will  be  in  nearly  linear  systems,  which  are  slowly  varying.  The 
user  is  cautioned  against  using  very  large  steps  with  the  RFB  method  since  gross 
errors  may  result  [2]. 


It  should  be  noted  that  the  diagonal  form  of  the  mass  matrix,  is 

used  in  each  of  the  above  methods.  The  advantages  in  storage  and  solution  effort 
are  obvious.  It  has  been  shown  that  with  proper  selection  of  the  integration 
parameters,  the  errors  introduced  by  the  lumped  mass  approximation  tend  to  com¬ 
pensate  for  those  induced  by  the  solution  algorithm  to  give  more  accurate  results. 
Experience  has  shown  that  a  =  0,  6  =  1/12  and  y  =  1/2  is  a  good  choice  of 
parameters.  Preliminary  observations  [2]  indicate  that  a  small  negative  value 
for  a  may  be  beneficial  in  controlling  amplitude  attenuations.  Values  of 
y  >  1/2  introduce  numerical  damping  of  the  higher  frequency  components  while  y  <  1/2 
introduces  negative  damping. 

No  internal  damping  terms  have  been  implemented  in  these  time  domain 
solution  forms. 

2.2.3  Frequency  Domain  Solution 
2.2. 3.1  Solutions  for  Regular  Waves 

The  response  amplitudes,  {Ul,  for  a  given  frequency  are  obtained  by 
solving  the  linear  simultaneous  algebraic  equations  represented  by  (2-69).  The 
coefficient  matrices  for  mass  and  damping  must  be  recalculated  for  each  frequency 
since  the  linearization  procedure  leaves  them  dependent  on  the  wave  frequency. 

The  damping  terms  present  some  further  difficulty.  In  addition  to  being 
frequency  dependent,  the  linearized  viscous  terms  are  dependent  on  the  amplitude 
of  the  response.  The  ship's  roll-damping  depends  on  the  roll  angle,  the  buoy 
rotational  damping  depends  on  the  rotation  angle,  and  the  cable  and  lumped  body 
damping  depend  on  the  lateral  displacement  amplitudes  (see  Appendix  E).  Thus 
it  is  seen  that  the  incremental  equations  are  not  strictly  linear. 

An  approximate  procedure  is  introduced  to  deal  with  this  problem.  This 
involves  iterative  solutions  of  the  equation  (2-69)  for  each  frequency.  The 
first  solution  at  a  given  frequency  is  calculated  assuming  a  ship's  roll  angle. 

The  roll  angle  obtained  from  the  solution  is  then  used  along  with  the  other 
pertinent  response  amplitudes  to  recalculate  the  damping  terms  and  obtain  another 
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solutlon.  This  procedure  is  repeated  until  two  successive  estimates  of  ship's 
roll  are  within  1®  of  each  other.  It  is  assumed  that  buoy  and  cable  damping  are 
less  important  than  ship's  roll  and  are  thus  converged  when  the  roll  has  converged. 
The  response  is  then  dependent  not  only  on  the  wave  frequency,  but  also  on  the 
wave  amplitude.  This  means  that  it  is  not  appropriate  to  assume  a  unit  amplitude 
for  a  given  wave  frequency  with  the  intent  of  obtaining  a  Response  Amplitude 
Operator  (RAO)  for  that  wave.  It  is  necessary  to  have  the  correct  wave  amplitude 
at  each  frequency.  Therefore,  the  sea  spectrum  must  be  used  in  the  calculation 
of  the  regular  wave  responses.  Once  the  steady  state  response  for  a  given  wave 
frequency  and  amplitude  is  obtained,  the  RAO  is  estimated  by  dividing  by  the 
wave  amplitude. 

The  program  allows  the  mass  matrix  to  be  formed  either  with  the  lumped 
or  consistent  form  of  the  element  mass.  The  lumped  form  is  used  for  the  fluid 
added  mass  terms. 

Provisions  for  internal  damr-'iig  effects  are  provided  in  the  proportional 
damping  form.  Thus 

[C]  *  a  [M]  +  6  [K.^]  (2-78) 

where  a  and  8  are  proportionality  constants.  The  damping  from  fluid  drag  effects 
on  the  cable  elements  must  be  linearized  before  it  can  be  used  in  equation  (2-69). 
This  linearization  is  described  in  Appendix  E. 

2. 2. 3. 2  The  Steady-State  Wave-Induced  Drift  Forces 

Whenever  waves  encounter  a  floating  body  there  arises  a  set  of  forces 
which  tend  to  move  the  body  in  the  direction  the  wave  is  traveling.  These  forces 
are  often  neglected  since  they  are  usually  small  in  magnitude.  These  are  the 
so-called  second  order  wave-induced  drift  forces.  They  are  generally  slowly 
varying  compared  to  the  frequency  of  the  incident  wave;  however,  they  have  an 
average  or  steady-state  component  which  may  be  significant  enough  to  cause  an 
adjustment  of  the  static  position  of  the  ship.  These  forces  are  directionally 
dependent  and  are  sensitive  to  the  amplitude  of  the  ship's  response  to  the  wave. 
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The  DTNSRDC  Ship's  Motion  File  provides  a  table  of  coefficients  which  can  be  used 
to  estimate  the  steady-state  drift  forces  after  the  ship's  dynamic  response  is 
obtained.  It  should  be  noted  that  these  forces  do  not  estimate  the  dynamic 
effects,  which  are  at  a  lower  frequency  than  ^e  incident  wave. 

The  specific  form  for  the  drift  forces  is  given  in  Appendix  D  in  the 
discussion  of  the  Ship  Motion  File.  The  result  is  a  set  of  forces  for  the 
lateral  and  longitudinal  directions  and  a  yaw  moment  acting  at  the  ships  reference 
point  for  each  wave  frequency.  It  is  assumed  that  the  drift  forces  are  cumulative 
for  the  various  wave  represented  in  a  wave  spectrum.  Therefore,  the  drift  forces 
for  each  of  the  regular  waves  are  accumulated.  It  is  felt  that  this  is  a  reason¬ 
able  approximation  since  the  wave  amplitude  indicated  by  the  wave  spectrum  is 
used  in  response  calculation. 

These  accumulated  drift  forces  can  then  be  used  as  an  additional  static 
loading  to  adjust  the  static  reference  state.  If  it  is  felt  that  this  adjustment 
will  affect  the  regular  wave  solutions  significantly,  the  user  can  request 
iterations  on  the  regular  wave  solutions  and  the  drift  force  adjustments  to  the 
static  reference  until  sufficiently  small  changes  are  found.  Either  the  MNR  or 
RFB  solutions  can  be  used  to  adjust  the  static  reference. 


2.2. 3. 3  Solution  Procedure  for  Random  Seas 

The  superposition  of  the  regular  wave  responses  to  represent  the  response 
to  random  seas  follows  the  well  established  methods  from  the  theory  of  random 
vibrations  [15].  The  sea  is  assumed  to  be  uni -directional  (long-crested)  and 
is  described  in  terms  of  a  generalized  spectral  energy  density  function  having 
the  following  form: 

S(u>)  =  A/u^  (2-79) 


Typical  values  for  the  parameters  A  and  B  are  given  in  Section  3.2.3.10. 


The  frequencies  to  be  included  in  the  set  of  regular  wave  calculations 
are  determined  by  specifying  a  frequency  increment.  Aw,  a  lower  bound  and 
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and  upper  bound,  Regular  wave  responses  are  then  calculated  for 

*  “uiin  J  *  (i-l)A(D 
1  max 

At  each  frequency  the  incident  wave  height  is  determined  from  the  sea  spectrum 
by  the  following: 

h(a.^)  =  (2-81) 

(Note  that  the  wave  height  is  twice  the  wave  amplitude.)  This  wave  height  is 
then  used  to  converge  on  a  steady-state  dynamic  response  using  the  methods 
described  in  the  previous  section. 

Let  H(ui^)  represent  the  response  of  one  of  the  quantities  (nodal  displace¬ 
ment  component  or  element  tension).  The  response  spectral  density  of  this  quantity 
is  then  given  by 

S^(«)  =  HH*  S{^)  (2-82) 


The  mean  square  of  the  response  is 


^00  CO 

E[x^]  =  ^  HH*  s(«)da,  =  ^  S^(w)da)  ’ 


(2-83) 


The  integral  is  evaluated  numerically  using  the  points  obtained  from  each  of 
the  regular  wave  solutions.  Thus 

E[x^]  I  S  (cojAu.  (2-84) 

i=l  ^  ^ 

where  N  is  the  number  of  regular  wave  components  used. 


When  the  response  quantities  represent  the  dynamic  excursions  relative 
to  the  static  reference  state,  then  their  expected  values  (i.e.,  their  means) 
are  zero.  The  magnitude  of  the  response  is  then  treated  as  a  static  part  plus 
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a  dynamic  part.  The  amplitude  of  the  dynamic  part  is  assumed  to  follow  a 
Rayleigh  distribution  with  a  mean-square  value  which  is  a  function  of  the  area 
under  the  response  spectrum  curve.  Assuming  the  sea  spectrum  used  is  based  on 
double  the  square  of  the  wave  height,  the  mean-square  of  the  response  amplitude 
is  the  value  obtained  by  equation  (2-84)  divided  by  eight.  It  is  possible,  then, 
to  make  statistical  estimates  of  the  maximum  response  by  making  statistical 
estimates  of  the  dynamic  part  and  adding  them  to  the  static  values  obtained  in  the 
updated  reference  configuration. 

2.2.4  Natural  Frequencies  and  Mode  Shapes 

The  incremental  equations  (2-36)  without  damping  or  external  forces 
has  the  form 

[M]  {u}  +  [Kj]  [u]  =  0  (2-85) 

If  it  is  assumed  that  simple  harmonic  motion  occurs  then 

{-J  [M]  +  [K^])  {u}  sin  ut  =  0  (2-86) 

since  {u)  sin  ut  =  0  only  in  the  trivial  case  of  {u}  =  0,  this  requires  that 

det  ([K^]  -  J  [M])  =  0  (2-87) 

Solutions  to  this  equation  lead  to  the  natural  frequencies  and  mode 
shapes  of  the  small  displacement  free  oscillations  about  the  configurations  used 
in  expressing  [M]  and  [K.^.].  An  iterative  solution  is  used  which  obtains  all  of 
the  natural  frequencies  and  mode  shapes  for  the  system. 

The  primary  assumption  made  in  getting  the  frequencies  and  mode  shapes 
is  that  the  equations  are  linear  and  no  significant  damping  exists.  Both  of  these 
assumptions  may  be  violated  to  some  extent  for  underwater  cable  structures.  The 
option  is  provided  since  the  ■f'requencies  and  mode  shapes  may  still  be  indicative 
of  the  structural  behavior. 


The  SEADYN/DSSM  program  provides  a  unique  feature  of  checking  the 
capacity  of  the  various  components  of  the  system  against  to  Imposed  loads. 

The  three  types  of  checks  provided  are: 

1.  Anchors  -  the  loads  Imposed  on  the  anchor  or  fixed  node  are 
summed  and  the  resultant  Is  compared  with  the  anchor  holding 
power. 

2.  Buoys  -  the  resultant  of  the  loads  In  the  lines  connecting 
to  the  buoy  Is  checked  against  the  buoy  capacity. 

3.  Lines  -  individual  cable  elements  are  checked  to  see  If 
they  are  loaded  beyond  their  capacities. 

The  component  checks  for  buoys  and  anchors  follow  the  procedures  outlined  In 
NAVFAC  DM-26  and  discussed  in  Reference  9.  The  component  capacities  can  be 
Input  or  obtained  from  the  Inventory  tables  developed  for  the  DESMOOR  program  [9]. 
The  Inventories  are  described  In  Appendix  C. 

All  of  the  adequacy  checks  rely  on  estimates  of  the  loads  at  the  end  of 
an  element.  The  one  dimensional  simplex  element  has  only  one  value  of  tension 
associated  with  It  regardless  of  how  long  It  Is,  what  its  weight  Is  and/or  how 
much  distributed  load  It  is  supporting.  The  tension  associated  with  the  element 
can  be  thought  of  as  being  at  the  midpoint  of  the  element.  The  procedure  for 
modeling  distributed  loads  must  be  recalled  to  see  how  to  estimate  element  end 
loads.  Equations  (2-9)  and  (2-24)  both  show  that  the  nodal  point  equivalent  loads 
for  distributed  loads  are  estimated  using  the  element  shape  functions.  These 
end  loads  represent  the  loads  applied  to  the  end  nodes  by  the  element  as  It 
supported  distributed  loads.  The  total  load  at  each  end  of  the  element  Is  then 
the  vector  sum  of  the  element  tension  and  the  element  loads.  The  direction  of 
the  resultant  load  at  each  end  gives  an  estimate  of  the  direction  the  cable  lies 
at  that  point.  In  symbolic  form  the  total  loads  applied  to  the  nodes  by  the 
element  are 
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Some  ambiguity  in  the  loads  applied  to  anchors  results  when  long  elements  are 
used  next  to  the  anchor.  If  the  line  actually  contacts  the  bottom  before  the 
anchor  connection  (i.e.,  between  the  two  nodes  which  define  the  element),  the 
SEADYN/DSSM  Program  will  not  detect  this  in  the  solution.  When  a  line  check 
is  made  the  resultant  of  the  forces  at  the  anchors  will  have  a  component  pulling 
down  indicating  the  line  was  hanging  below  the  anchor.  This  is  due  to  the  model 
not' being  able  to  sense  the  bottom  contact  and  transfer  the  weight  at  that  point. 
Detailed  and  accurate  modeling  of  the  line  interaction  with  the  bottom  requires 
short  elements  in  that  region.  Fortunately  the  lack  of  modeling  detail  at  the 
bottom  of  the  line  has  little  influence  on  the  response  at  the  top  of  the  line. 


2.2.6  Surface  and  Bottom  Constraints 

The  foregoing  discussion  introduces  the  problem  of  constraining  the 
model  of  the  system  to  responses  which  lie  between  the  natural  boundaries  imposed 
by  the  water  surface  and  the  bottom.  One  does  not  expect  buoys  or  lines  to  rise 
out  of  the  water  or  lines  and  anchors  to  go  below  the  bottom.  The  SEADYN/DSSM 
program  assumes  the  surface  and  bottom  are  flat  and  parallel .  Checks  are  made  at 
each  step  of  static  and  time  domain  analyses  to  see  if  nodes  of  the  system  are 
within  the  imposed  limits.  To  avoid  the  costly  operation  of  checking  all  nodes 
in  the  system  at  every  step  it  is  presumed  that  the  critical  points  are  where 
lumped  bodies  are  located.  Therefore,  only  those  nodes  are  checked.  When  one 
of  the  critical  nodes  is  within  a  certain  tolerance  distance  of  the  surface  or 
bottom  it  is  constrained.  For  buoys  or  floats  this  means  that  the  node  is  held 
fixed  in  the  vertical  direction  but  free  in  the  lateral  directions.  All  three 
components  are  fixed  for  anchors.  Whenever  the  vertical  resultant  of  all  of  the 
element  tensions  connecting  to  the  point  exceeds  the  sum  of  the  external  vertical 
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loads  at  that  point,  then  the  constraint  is  released.  The  external  loads  are 
assumed  to  include  the  distributed  loads  from  the  elements  and  the  weight  or 
buoyant  force  from  the  lumped  body. 

2.2.7  Initial  Configuration  Problems  and  Input  Generation  Schemes 

One  of  the  frustrating  features  of  cable  analysis  is  that  most  cable 
systems  do  not  have  rigidity  or  spatial  stability  unless  preloads  are  imposed. 
The  systems  are  usually  so  flexible  that  small  changes  in  preloads  cause  large 
shape  changes.  Noting  that  all  of  the  static  solution  methods  use  a  stiffness 
matrix,  one  is  faced  with  a  problem  of  getting  a  realistic  estimate  of  the 
stiffness  matrix  which  is  nonsingular.  In  many  situations  it  is  necessary  to 
have  very  accurate  estimates  of  the  initial  configuration  before  any  of  the 
solution  methods  will  work. 


Various  facets  of  this  problem  are  explored  in  Reference  2.  Some 
procedures  which  may  help  obtain  stable  starting  configurations  are  discussed  in 
Section  3.5. 


A  technique  that  has  proven  of  some  use  in  overcoming  ill -conditioning 
and  even  singularities  when  using  the  MNR  solution  is  the  use  of  numerical  damping. 
Felippa  [16]  shows  that  nonsingular  adjustment  to  the  estimator  matrix,  [1^  can  be 
generated  by  adding  a  matrix  of  the  form 

u  B  [1] 

where 

B  =  {R}  ^[K]  {R}/{R}  ^{R} 

and  u  in  a  user  specified  numerical  damping  coefficient.  This  additional  term 
tends  to  "stiffen"  the  estimator  matrix  and  increase  the  chances  for  convergence. 
Section  3.2. 3.2  gives  some  suggestions  for  selecting  u. 


The  program  also  provides  for  numerical  damping  to  be  used  with  the 
incremental  solutions.  This  feature  should  be  used  with  extreme  care  since  it 
alters  the  equations  of  equilibrium.  If  an  appropriate  value  can  be  selected 


it  would  be  possible  to  compensate  for  some  of  the  error  in  the  first  incremental 
step.  There  is  no  rational  way  available  to  estimate  how  much  damping  to  use 
in  this  case. 

The  program  provides  a  quick  and  convenient  way  of  getting  starting 
configurations  when  negatively  buoyant  lines  are  used.  If  a  line  between  two 
defined  points  is  a  catenary  which  reaches  the  bottom  with  a  horizontal  tangent 
somewhere  between  the  two  points,  then  nodes  can  be  generated  along  that  line.  The 
details  of  the  procedure  are  given  in  Section  3.2. 2. 5.  The  well  known  catenary 
equations  are  used  and  one  or  more  of  the  generated  nodes  are  assumed  to  be  on 
the  bottom  if  the  lower  defining  node  is  not  at  the  tangent  point.  The  element 
preloads  for  the  line  are  also  generated  and  can  be  assigned  to  the  line  elements. 


2.2.8  Condensed  Equation  Format 


Whenever  the  solution  procedure  calls  for  a  global  stiffness  matrix, 
a  reduced  form  of  the  matrix  is  generated.  First  it  is  assumed  that  the  matrix 
is  symmetric  so  that  only  terms  in  the  upper  triangle  need  be  generated  and 
stored.  Secondly,  the  equations  are  assumed  to  be  sparse  and  ordered  in  such  a 
way  that  the  nonzero  terms  lie  within  a  band  near  the  diagonal.  The  half-bandwidth 
(the  number  of  terms  from  the  diagonal  to  the  end  of  the  populated  band)  is  given 
by  three  times  the  largest  difference  between  active  node  numbers  on  any  element 
in  the  system.  The  matrix  is  then  stored  in  rows  which  begin  at  the  diagonal 
and  extend  to  the  end  of  the  band.  The  consistent  mass  matrix  and  damping 
matrix  are  also  generated  in  this  banded  form  for  the  frequency  domain  analysis. 


The  simultaneous  equation  solution  routines  are  formulated  to  operate 
on  the  matrix  in  this  condensed  form  and  the  full  matrix  is  never  generated. 

A  Choleski  decomposition  procedure  is  used  to  solve  the  equations.  The  necessary 
accumulations  of  products  associated  in  that  procedure  are  done  in  double  precision 
to  help  control  numerical  errors.  Solutions  of  the  frequency  domain  equations  are 
done  in  complex  arithmetic  with  the  real  and  imaginary  parts  summed  separately  in 
double  precision. 


The  compact  storage  and  special  solution  routines  provide  a  significant 
savings  in  program  size  and  the  number  of  operations  required  to  solve  the  equations. 


3.0 


INPUT  INSTRUCTIONS 


3.1  PROGRAM  OVERVIEW 

The  SEADYN/DSSM  Program  is  constructed  to  allow  a  high  level  of 
versatility  in  the  types  of  analysis  performed  and  the  detail  to  which  they  are 
pursued.  It  is  a  finite  element  program  which,  at  present,  has  only  one  element 
type  (the  nonlinear,  one-dimensional  simplex  element  in  three-dimensional  space) 
plus  provisions  for  surface  ships  and  buoys  and  lumped  bodies  (anchors,  weights, 
buoys,  etc.).  Various  types  of  analyses  can  be  performed,  and,  in  some  cases,  more 
than  one  solution  procedure  is  provided  for  a  given  type  of  analysis.  Which 
solution  procedure  should  be  used  is  very  dependent  on  the  problem  being  analyzed 
and  the  level  of  accuracy  desired.  The  variety  is  provided  with  the  hope  that  as 
the  user  becomes  more  familiar  with  the  program  and  the  nature  of  the  problem  he 
is  trying  to  solve,  he  may  be  able  to  "tune"  the  solution  procedure  to  the  problem 
and  obtain  the  desired  accuracy  at  a  minimum  cost. 

The  nonlinear  nature  of  many  cable  problems  and  the  complexity  of  the 
systems  and  the  loads  imposed  on  them  put  a  premium  on  versatility  and  flexibility. 
Rather  than  attempting  to  foresee  all  possible  analysis  situations  and  build  in 
general  purpose  logic  to  deal  with  them,  the  program  has  been  constructed  with 
some  basic  capability  and  then  provision  is  made  to  allow  the  user  to  write  his 
own  logic  when  the  built-in  capability  does  not  meet  his  needs.  Thus  in  the  case 
of  drag  coefficients  for  the  cables  and  lumped  bodies,  the  spatial  distribution 
of  fluid  flow  velocities,  and  the  time  variation  functions,  the  user  can  take 
the  limited  capability  provided  or  he  can  write  his  own  subroutine  to  generate 
the  data.  Ship's  characteristics  and  loading  data  are  assumed  to  be  provided  in 
the  form  of  data  files.  Mooring  component  inventories  used  in  adequacy  checks 
are  contained  in  a  subroutine  which  can  be  altered  by  the  user. 

Analysis  of  nonlinear  cable  systems  requires  some  basic  understanding 
of  the  behavior  of  such  systems.  It  is  not  the  intent  of  this  document  to  provide 
that  understanding,  but  some  basic  principles  can  be  stated.  The  major  nonlinearity 
in  cable  systems  is  usually  geometric  in  nature.  This  means  the  structure  changes 
its  characteristics  as  it  deforms  such  that  the  undeformed  configuration  cannot 


be  used  to  write  the  equations  of  equilibrium.  This  leads  to  the  requirement  for 
treating  the  loading  in  stages.  For  example,  it  is  usually  necessary  to  obtain 
a  solution  for  the  configuration  of  the  system  under  the  action  of  gravity  (Dead 
load  analysis)  before  one  can  analyze  the  effects  of  currents  or  imposed  motions, 
etc.  Indeed,  it  is  often  necessary  to  apply  a  given  load  in  a  series  of  increments 
in  order  to  get  a  convergent  solution. 

One  basic  principle  to  observe  is  that  the  structure  may  not  be  stable 
unless  adequate  preloading  is  provided.  Unfortunately,  it  is  not  usually  possible 
to  estimate  the  shape  of  a  system  under  preloading  with  sufficient  accuracy  to 
assure  a  stable  numerical  description.  In  the  case  of  an  initial  dead  loaded 
configuration  it  may  be  possible  to  start  from  an  unstable  "guess".  This  can  be 
done  with  numerical  damping  and/or  some  of  the  techniques  described  in 
Section  3.5. 

The  program  treats  a  given  problem  as  a  series  of  subanalyses.  This 
usually  follows  the  order  listed  below: 

DEAD  LOAD  -  determine  the  initial  static  configuration  with 
only  gravity  loads  acting.  (This  may  include 
specified  point  loads  to  simulate  reactions  at 
certain  "adjustable"  points.) 

LIVE  LOAD  -  change  of  configuration  due  to  current  and/or 
point  loads. 

DYNAMIC  LOAD  -  Dynamic  response  to  time  varying  currents, 
point  loads,  and/or  imposed  movement. 

OR 

Frequency  domain  response  of  mooring  systems 
to  surface  waves. 

The  configuration  at  the  end  of  one  subanalysis  is  presumed  to  be  the  starting 
configuration  for  the  subsequent  subanalysis.  These  subanalyses  can  be  skipped 
or  repeated  as  long  as  the  user  has  specified  consistent  starting  data  for  the 
next  subanalysis. 


A  RESTART  option  is  provided  which  can  be  used  to  save  data  for  various 
stages  of  the  analysis  for  later  re-use.  This  option  can  either  be  used  to 
resume  calculations  on  a  later  run  or  it  is  possible  to  return  to  a  previous 
configuration  within  the  same  run  and  initiate  a  new  sequence  of  subanalyses. 
Details  of  the  RESTART  option  are  given  in  Section  3.2.4. 

Another  feature  provided  is  an  eigenvalue  routine  which  can  be  called 
at  the  end  of  any  of  the'subanalyses  to  obtain  an  estimate  of  the  natural  modes 
and  frequencies  at  the  current  configuration. 

A  special  option  is  provided  which  allows  the  adequacy  of  the  mooring 
components  to  be  checked.  Component  inventories  may  be  referenced  for  the  capacity 
of  the  component  or  the  capacity  can  be  input  with  the  request  for  the  check. 

The  SEADYN/DSSM  Program  allows  the  user  to  choose  a  coordinate  system 
which  is  convenient  for  him.  There  are  only  two  general  restrictions.  First, 
the  global  system  must  be  a  right-handed  rectangular  cartesian  system.  Second, 
whenever  gravity  loads  are  used,  one  of  the  global  axes  must  be  parallel  to  the 
direction  of  gravity.  The  vertical  axis  may  point  either  up  or  down. 

The  program  develops  coordinate  transformations  from  the  global  system 
to  local  element  coordinate  systems.  The  local  system  for  a  cable  element  is 
defined  by  the  order  and  the  current  position  of  the  nodes  given  to  define  it. 

The  local  x  axis  runs  from  the  first  defining  node  through  the  second.  The 
remaining  two  local  axes  are  automatically  selected  to  j'it  the  purposes  of  the 
specific  calculations. 

The  local  systems  for  mooring  buoys  and  surface  ships  require  more 
careful  consideration  by  the  user.  In  each  of  these  cases  the  initial  orientation 
of  the  local  system  is  implied  by  nodal  input  data.  The  local  coordinate  system 
for  a  surface  ship  assumes  the  x  axis  is  parallel  with  the  longitudinal  axis  of 
the  ship  having  the  origin  at  the  ships  reference  point  (c.g.)  and  the  positive 
direction  pointing  aft.  The  local  y  axis  points  positive  starboard  and  the 
positive  z  direction  is  up.  The  initial  angle  between  the  global  system  and 
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the  local  ship  system  must  be  given  as  input  on  the  appropriate  node  card  (see 
Section  3. 2. 2. 5).  Ships  require  two  nodes  to  define  the  reference  point  (one 
for  position  and  one  for  angles)  plus  one  additional  node  for  each  point  where 
a  line  is  attached  -  or  motion  output  is  desired.  These  attachment  points  are 
designated  as  slaves  to  the  reference  point. 

Mooring  buoys,  or  riser  buoys,  are  characterized  by  a  rigid  body  interposed 
between  the  two  cable  "e-lements  and  constrained  to  ride  on  the  surface.  Mooring 
buoys  are  distinct  from  lumped  buoys  since  they  require  four  node  numbers  to  de¬ 
fine  them  while  a  lumped  buoy  requires  only  one.  The  four  nodes  are  used  as 
follows:  one. for  each  of  the  two  attachments  and  two  for  the  buoy  center  point 
and  its  orientation.  The  local  coordinate  system  for  a  mooring  buoy  is  generated 
from  the  positions  of  the  center  and  the  attachments.  The  attachment  nodes  must 
be  slaved  to  the  center  node.  The  dynamic  equations  used  in  the  frequency  domain 
analysis  of  mooring  buoys  presumes  that  the  attachment  nodes  are  co-planar  with 
the  center  node.  Since  no  checks  are  made  for  this  situation,  the  user  must  be, 
careful  to  check  that  himself. 

There  is  no  system  of  units  implied  in  the  program's  calculations.  There¬ 
fore,  it  is  necessary  to  specify  all  input  data  in  a  consistent  set  of  units.  All 
spatial  measures  (length,  diameter,  area,  etc.)  must  be  given  in  the  same  units. 
Likewise,  all  force  data  must  be  consistent.  For  example,  if  the  nodal  coordinates 
are  given  in  feet  and  the  applied  loads  are  given  in  pounds,  then  element  diameters 
must  be  given  in  feet  and  the  load/strain  data  must  be  given  in  pounds.  Further- 
more,  specific  weights  would  then  be  given  in  Ibs/ff;.  If  the  time  unit  used  is 

seconds,  then  the  acceleration  due  to  gravity  would  have  the  units  of  ft/se<^and 
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mass  would  have  the  units  Ib-ssc  /ft.  Input  of  conversion  factors  is  provided  to 
allow  use  of  load  files,  motion  files  and  inventories  which  have  different  units 
than  those  assumed  in  the  run. 

Each  analysis  performed  by  the  program  requires  three  sets  of  cards: 

ANALYSTS  TITLE  CARD:  An  alpha-numeric  card  used  to  label  the- output  and 

identify  the  analysis. 

Revised  JULY  1978 
Revised  SEPT  1378 


SYSTEM  DESCRIPTION  DATA:  A  set  of  cards  describing  the  system  being 

analyzed  and  giving  the  initial  configuration. 
Provides  data  on  the  element  materials  and 
connectivity,  ship's  characteristics,  lumped 
body  data,  boundary  conditions  and  fluid 
media  characteristics.  These  cards  may  be 
superseded  by  a  RESTART  procedure  which  obtains 
this  data  from  a  previously  saved  file. 

SUBANALYSIS  OPTION  DATA:  A  set  consisting  of  two  or  more  cards  which 

identifies  the  analysis  to  be  performed  and 
the  procedures  to  be  used.  In  addition,  the 
specific  data  required  (e.g.,  loads,  imposed 
motion,  spectral  data,  etc.)  are  given.  Any 
number  of  subanalysis  data  sets  may  be 
associated  with  one  system  description  set. 
Remember  the  need  for  sequentially  consistent 
subanalysis  states. 

The  specific  forms  these  data  sets  take  are  described  in  the  next 
section.  Section  3.3  describes  the  optional  user-supplied  subroutines.  Problem 
size  limitations  are  discussed  in  Section  3.4.  The  expressions  used  in  calculating 
the  default  drag  coefficients  (those  assumed  by  the  program)  are  sunmarized  in 
Section  3.6.  A  description  of  the  component  inventories  is  provided  in  Section 
3.7. 

3.2  INPUT  CARD  DESCRIPTIONS 

3.2.1  Analysis  Title  Card  (12A6)* 

VARIABLE 

COLS.  NAME  DESCRIPTION 

1-72  HED  Any  descriptive  title  (see  RESTART  option 

for  special  usage). 


♦FORTRAN  format  specification 


3.2.2  System  Description  Data 
3. 2. 2.1  Master  Control  Card  (1415) 


1-5 

NN 

Total  nimber  of  nodal  points 

6-10 

NFIX 

Number  of  fixed  nodes 

11-15 

MVB 

Number  of  moved  nodes 

16-20 

NSLAVE 

Number  of  slave  nodes 

21-25 

NE 

Number  of  cable  elements 

26-30 

MATT 

Number  of  cable  material  types  (10  max) 

31-35 

IBUO 

Number  of  lumped  bodies  (buoys/anchors)  (30  max) 

36-40 

NSHIPS 

Number  of  ships  (3  max) 

41-45 

JDLD 

Dead  load  flag  (see  note  3) 

46-50 

JDYN 

Dynamic  flag  (see  note  4) 

51-55 

IBG 

Output  flag  j^O  prints  nodal  masses,  gravity 
loads,  element  slopes,  initial  and  unstrained 
element  lengths 

56-60 

NGEN 

Number  of  node  generation  cards  provided 
(see  Section  3. 2. 2. 6) 

61-65 

NSFILE 

Number  of  ships  on  data  file 

>0  means  read  data  to  create  the  file 
<0  means  file  already  for  use. 

66-70 

INVY 

Component  inventory  flag 

0  inventory  acceptable  or  not  used 

1  inventory  unit  conversion  required 

NOTES: 

1.  The  node  numbering  scheme  is  not  arbitrary.  The  following 

numerical  order  must  be  observed: 

NODES  WITH  ACTIVE  DEGREES  OF  FREEDOM 

these  are  the  basic  unknown  degrees  of 
freedom.  Any  node  which  is  not  a  SLAVE 
and  has  one  or  more  motion  component, 
to  be  solved  for. 

MOVED  NODES 

nodes  with  motion  completely  described  in 
a  dynamic  analysis.  These  nodes  are 
assumed  fixed  in  static  analyses. 


FIXED  NODES 

boundary  nodes  which  are  completely  fixed 
at  all  times  in  all  subanalyses. 

SLAVE  NODES 

nodes  with  movement  specified  as  a  function 
of  another  node  (MASTER  NODE)  through  a 
rigid  link  transformation.  The  MASTER  NODE 
is  one  of  the  active  nodes. 

In  most  solution  procedures  a  bandwidth  dependent  scheme  is 
employed.  Active  nodes  should  be  numbered  in  such  a  way  that 
bandwidth  is  minimized.  This  numbering  scheme  is  left  up  to 
the  user  tto  define. 

The  dead  load  flag  controls  the  use  of  gravity  loads.  Greater 
than  zero  means  the  gravity  load  is  incremented  from  zero  to 
its  full  value  in  a  DEAD  LOAD  subanalysis  and  has  its  full 
value  constant  in  any  LIVE,  DYN  or  FREQ  subanalysis.  Zero 
means  the  system  is  weightless  (DEAD  subanalysi's  not  to  be 
called  for).  Less  than  zero  means  the  gravity  loads  will  be 
incremented  along  with  the  currents,  etc.,  in  a  LIVE  sub¬ 
analysis;  the  full  gravity  loads  will  be  used  in  any  DYN  or 
FREQ  subanalyses. 

The  dynamic  flag  controls  the  generation  of  the  system  mass 
matrix.  Zero  means  no  mass  matrix  will  be  formed  (DYN  or 
FREQ  subanalysis  not  allowed).  Non-zero  means  a  lumped  mass 
matrix  will  be  generated. 

Up  to  three  ships  can  be  involved  in  the  system,  except  for  the 
FREQ  subanalysis.  Only  one  ship  can  be  involved  in  that  case. 


3. 2. 2. 2  Parameter  Card  (215,  7E10.0) 


COLS. 

VARIABLE 

NAME 

1-5 

IDIR 

6-10 

INDRAG 

11-20 

GRAV"" 

21-30 

G1 

31-40 

G2. 

41-50 

VAIR 

51-60 

GAIR 

61-70 

SURFCE 

71-80 

BOTTOM 

DESCRIPTION. 

Direction  of  gravity  (signed  integer 
indicating  which  of  the  global  axes 
gives  the  gravity  direction) 

Drag  model  option 

fO  calls  for  subroutine  DRAGCO  (see  Section  3.3) 

Magnitude  of  gravitational  acceleration  (LT  )* 
Typical  value  32.17  ft/sec^ 

Kinematic  viscosity  of  fluid  (L^  T”^) 

Typical  value  1.77  x  10“5  ft^/sec 

Specific  weight  of  fluid  (FL"^) 

Typical  value  64.0  Ib/ft^ 

2  -1 

Kinematic  viscosity  of  air  (L  T  ) 

Typical  value  1.68  x  10"^  ft^/sec 

Specific  weight  of  air  (FL  ) 

Typical  value  0.0765  Ib/ft^ 

Surface  Coordinate  (position  on  gravity  axis)  (L) 
Bottom  coordinate  (position  on  gravity  axis)  (L) 


NOTES 

T.  It  is  assumed  that  the  global  axes  are  selected  such  that  one  of  them 
coincides  with  the  direction  of  gravity.  Thus  IDIR  =-2  means  that 
gravity  acts  in  the  negative  y  direction  (i.e.,  +y  is  up) 

2.  INDRAG  can  be  used  to  select  user-specified  drag  coefficients, 
since  it  is  passed  to  the  DRAGCO  subroutine  whenever  it  is  called. 

If  it  is  zero  it  overrides  all  other  drag  coefficient  specifications, 
and  the  program  will  use  the  built  in  or  default  values  described 

in  Section  3.6. 

3.  If  either  SURFCE  or  BOTTOM  is  non-zero  the  motion  of  any  node  where 
there  is  a  lumped  body  will  be  limited  to  stay  between  SURFCE  and 
BOTTOM.  When  the  body  is  on  the  surface  it  has  only  the  vertical 
component  fixed.  When  on  the  bottom  it  may  have  all  three  components 
fixed  or  only  the  vertical  component.  (See  JanCR  flag  of  section 

3. 2. 2.9.)  The  booy  remains  constraineo  to  tne  surface  or  bottom  until 
the  element  forces  are  sufficient  to  move  it  away.  (See  note  8  of 
section  3. 2. 2. 9.) 

*■(  )  gives  units  expected 


REVISED  SEPT  1978 


4.  Input  for  the  fluid  and/or  air  terms  is  required  only  when  the 
effects  of  the  fluids  are  to  be  considered. 

3.2.2. 3  Inventory  Conversion  Factors  (3E10.0) 

Read  only  when  INVY/O  (see  MASTER  CONTROL  CARD) 


VARIABLE 

NAME 

DESCRIPTION 

1-10 

FRCVY 

Inventory  force  conversion  factor 
multiplies  weights  and  strengths 

11-20 

FLNVY 

Inventory  length  conversion  factor 
multiplies  buoy  diameters  and  lengths 

21-30 

FDIVY 

Inventory  diameter  conversion  factor 
multiplies  line  diameters 

iOTES 

1.  A  zero  input  or  blank  field  means  the  factor  is  1.0  (i.e.,  no 
conversion  required). 


In  each  case  the  value  obtained  from  the  inventory  is  multiplied 
by  the  given  factor  to  get  it  into  the  units  implied  by  the  rest 
of  the  input  data. 
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3. 2. 2. 4  Ship  Load  File  Data 

This  set  of  data  consists  of  a  variable  number  of  cards  which  describe  the 
ships  and  their  load  function  tables  which  are  to  be  placed  on  the  ship's  reference 
file.  These  cards  are  read  only  when  NSFILE  >  0  on  the  CONTROL  CARD.  The 
following  cards  will  be  input  for  each  of  the  ships  indicated  by  NSFILE.  The 
number  of  ships  that  can  be  input  is  controlled  only  by  the  amount  of  file  space 
available.  If  NSFILE  =  1  the  data  is  read  and  remains  in  core  (the  ship's 
reference  file  is  not  set  up  in  this  case).  The  program  assumes  the  reference 
file  is  on  logical  unit  10. 

A  complete  ship  description  consists  of  the  following: 

SHIP  LOAD  TITLE  CARD 
UNIT  LABEL  CARD 
SHIP  PARAMETER  CARD 
WIND  CARD 

WIND  HEADING  CARD(S) 

WIND  FORCE  COEFFICIENT  CARDS 
SURFACE  CURRENT  CARD 
SURFACE  CURRENT  HEADING  CARD(S) 

SURFACE  CURRENT  FORCE  COEFFICIENT  CARDS 


SHIP  LOAD  TITLE  CARD  {12A6) 

VARIABLE 

COLS.  NAME  DESCRIPTION 

1-72  SHPCAP  Any  descriptive  title 

UNIT  LABEL  CARD  (A6.  4X,  A6,  4X.  A6,  4X,  A6) 

VARIABLE 

COLS.  NAME  CESCRIPTION 

1-6  WLBL,  Wind  Force  Label  (e.g.,  "TONS",  "POUNDS") 

11-16  CLBL  Current  Force  Label 

21-26  LLBL  Length  label  (e.g.,  "FEET".  "METERS") 

31-36  VLBL  Velocity  label  (e.g.,  "KNOTS,  "FT/SEC") 

NOTE:  These  labels  are  output  with  the  ship  data  as  a  reminder  of  the 
units  used.  They  are  used  for  no  other  purpose. 


SHIP  PARAMETERS  (8E10.0) 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-10 

TSLT 

Total  ship  length 

11-20 

TSAE 

End  projected  wind  area 

21-30 

TSAS 

Side  projected  wind  area 

31-40 

TSWL 

Water  line  lengti) 

41-50 

TSB 

Beam  at  mid  ships 

51-60 

TSD 

Draft  at  mid  ships 

61-70 

TSDSP 

Volume  displacement 

71-80 

TSAP 

Propeller  projected  area 

WIND  CARD  (215,  6E10.0) 

1-5  NWIND 

6-10  NTHETW 

11-20  SCALE 

21-30  WNDVEL(l) 


No.  of  wind  velocity  tables,  maximum  Is  5. 
No.  of  heading  In  each  table,  maximum  Is  20 
Test  S  ' lie  factor  (A  means  1/Ath  scale). 
First  wind  velocity  (smallest) 


61-70  WNDVEL(5) 


Fifth  wind  velocity 


(Repeat  as  required  to  get  NTHETW  entires.) 


COLS. 

1-10 
n-20 

NOTES: 

1.  Headings  should  be  between  0°  and  360°,  listed  from  the  smallest  to 
the  largest. 

2.  If  the  largest  value  is  180°  the  loading  functions  are  assumed  to 
be  symmetric  about  180®  for  the  end  forces  and  skew- symmetric  for 
the  side  forces  and  yaw  moments. 

3.  The  angle  is  measured  relative  to  the  ship's  local  coordinate  system. 


VARIABLE 

NAME  DESCRIPTION 

WNDHEO(I)  First  wind  heading  (Degrees) 

etc. 
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WIND  FORCE  COEFFICIENT  CARDS  {3E10.0) 


One  card 

for  each  heading  repeated  for  each  velocity. 

COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-10 

WNDCOE  (I, 

1.  J) 

End  force  coeff.  for  Ith  heading  and  Jth 
velocity. 

11-20 

WNDCOE  (I, 

2.  J) 

Side  force  coeff.  for  Ith  heading  and  Jth 
velocity. 

21-30 

WNDCOE  (I, 

3.  J) 

Moment  force  coeff.  for  Ith  heading  and  Jth 
velocity. 

(I  varies  before  J). 

SURFACE  CURRENT  CARD  (215. 

6E10.0) 

1-5 

NCRNT 

No.  of  current  tables. 

6-10 

NTHETC 

No.  of  headings  in  each  current  table 

11-20 

TDEPTH 

Test  water  depth  ■  (L) 

21-30 

CURVEL(l) 

First  current  velocity  (smallest)  (LT”^) 

€1-70 

CURVEL{5) 

Fifth  current  velocity 

(Repeat  as  required  to  get  NTHETC  entries.) 

VARIABLE 

COLS.  NAME  DESCRIPTION 

I-IO  CURHED(l)  First  Current  Heading  (Degrees) 

n-20  etc. 

(See  notes  for  WIND  HEADING  CARO.) 

SURFACE  CURRENT  FORCE  COEFFICIENT  CAROS  (3E10.0) 

(One  card  for  each  heading  repeated  for  each  velocity.) 

1-10  CURC0E(I,  1,  J)  End  force  for  Ith  heading  and  Jth  velocity 

11-20  CURC0E(I,  2,  J)  Side  force  for  Ith  heading  and  Jth  velocity 

21-30  CURC0E(I,  3,  J)  Moment  force  for  Ith  heading  and  Jth  velocity 

(I  varies  before  J.) 


3. 2. 2. 5  Nodal  Point  Cards  (215,  3E10.0,  315,  5X,  A6) 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-5 

I 

Node  number 

6-10 

KSKP 

Node  generation  code  (see  notes) 

11-20 

X0(3*I-2) 

X 

21-30 

X0( 3*1-1) 

y  Global  Coordinates  of  Node  (I) 

31-40 

X0(3*I) 

z 

41-45 

N0DFIX(3*I-2) 

X  Constraint  codes  for  components 

46-50 

N0DFIX( 3*1-1) 

y  0  =  free  component 

51-55 

N0DFIX(3*I) 

1  =  fixed  component 

z  -I  =  means  buoy  at  this  node  is  on 

61-66 

DFLAG 

the  surface  and  this  component 
attached  to  surface  until  buoy 
pulled  under 

-M  where  M  ^  I  means  this  node  is 
slave  to  master  node  M 

"DEGREE"  when  nodal  coordinates 

are  angles  to  give  headings  for 
ships  and  the  angles  are  in  degrees. 
Otherwise  the  angles  are  assumed  to 
be  in  radians. 

NOTES 

1.  Input  of  nodes  need  not  be  in  numerical  order,  but  the  last  node 
(highest  node  number)  terminates  the  reading  of  nodes  and  must 
be  input  last.  Omitted  nodes  can  be  generated  using  a  straight 
line  with  uniform  spacing  with  the  aid  of  the  KSKP  parameter. 

If  KSKP  ^  0  on  a  node,  then  that  node  is  designated  as  the  last 
node  on  the  line  and  the  one  input  preceding  it  is  designated 
as  the  first  node  on  the  line.  Nodes  are  generated  evenly  along 
the  line  between  these  two  points  with  node  numbers  incremented 
by  KSKP  from  the  first  to  the  last  node.  The  difference  between 
the  node  numbers  at  the  ends  of  the  line  must  be  an  integer 
multiple  of  KSKP,  and  KSKP  cannot  be  negative.  The  generated 
nodes  will  have  the  same  constraint  codes  as  the  first  node  on 
the  line.  More  general  node  generation  schemes  are  available 
in  Section  3. 2. 2. 6. 


DSSM-2 


( 

I 
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2.  All  of  the  NN  nodes  must  be  accounted  for  in  the  combined 
specifications  for  nodes  and  generation  schemes.  See 
MASTER  CONTROL  card  for  notes  on  numbering  restrictions. 

3.  Nodes  which  define  the  position  of  ships  or  mooring  buoys 
must  allow  for  two  nodes  in  succession.  The  first  node 
defines  the  spatial  position  and  the  node  immediately 
following  gives  the  headings  of  the  local  coordinate  systems 
for  the  bodies.  (See  Section  3.1  for  a  discussion  of 
coordinate  systems.)  Since  both  the  ships  and  mooring 
buoys  can  introduce  very  soft  components  in  combination 
with  very  stiff  components  it  is  recommended  that  these 
nodes  be  numbered  at  the  end  of  the  active  nodes  to 
minimize  ill -conditioning  effects  in  the  solution. 

4.  All  attachment  points  on  a  ship  must  be  slaved  to  the  node 
defining  the  position  of  the  ship.  The  program  imposes  no 
limit  to  the  number  or  relative  location  of  these  attachments. 

5.  Mooring  buoy  frequency  domain  equations  presume  only  two 
attachments  in  addition  to  the  center  node.  Both  of  these 
attachment  nodes  must  be  slaved  to  the  node  defining  the 
position  of  the  center  node.  The  node  numbers  for  these 
attachment  points  are  also  required  on  the  input  for  the 
mooring  buoy  (Section  3. 2. 2. 8).  It  is  possible  to  slave 
more  attachments  to  the  buoy  center  node,  but  two  attachments 
must  be  identified  as  primary  ones  and  their  node  numbers 
provided  on  the  buoy  cards  if  a  frequency  domain  solution  or 
component  adequacy  check  is  requested. 

6.  The  initial  orientation  of  a  ship  relative  to  the  global 
coordinate  system  must  be  input.  The  angle  required  in  this 
definition  can  be  seen  in  the  following  sketch.  (The  global 
axes  are  x,  y,  z  and  the  local  ships  axes  are  1,  2,  3.) 
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7,  Displacement  constraint  codes  (except  master  node  flags)  are 
intended  for  use  only  on  ACTIVE  nodes,  i.e.,  those  node  numbered 
before  MOVED,  FIXED  and  SLAVE  nodes.  (See  note  1  Section  3. 2. 2.1.) 
They  produce  displacement  fixity,  but  they  do  not  reduce  the  number 
of  equations.  Globally  fixed  nodes  can  be  numbered  as  FIXED  nodes 
to  remove  them  from  calculations.  MOVED  nodes  have  an  implied  FIXED 
status  in  all  but  DYN  subanalyses.  Nodes  with  lumped  bodies  may 
change  their  constraint  codes  as  surface/bottom  interaction  is 
treated. 


g.  Slave  nodes  need  not  have  elements  connected  to  them.  This  allows 
the  investigation  of  the  response  of  specific  points  on  a  ship  at 
locations  other  than  the  attachments  of  mooring  and  working  lines. 


Node  Generation  Cards  (715,  5X,  2E10.0) 

Required  only  if  NGEN  f  0  (Master  Control  Card) 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-5 

NUMN 

Number  of  nodes  to  be  generated 

6-10 

N1 

Node  at  beginning  of  line 
(lowest  node  on  a  catenary) 

11-15 

N2 

Node  at  end  of  line 
(highest  node  on  a  catenary) 

16-20 

KSKP 

Node  number  increment 

21-25 

NBC 

Boundary  Condition  Option 

0  -  do  not  copy  constraint  codes  but 
set  all  generated  nodes  to  zero 

1  -  copy  from  node  N1 

2  -  copy  from  node  N2 

26-30 

NFIRST 

First  generated  node 
(default  =  N1  +  KSKP) 

31-35 

NCAT 

Generation  Code 

41-50 

GCAPW 

0  -  straight  line 

1  -  Catenary  with  lower  limit 

2  -  Catenary  without  lower  limit 
Weight/ length  of  generated  catenary  i 

51-60 

GCATH 

Horizontal  component  of  tension  (F) 

Revised  JULY  1978 


Revised  SEPT  1978 
Revised  OCT  1978 
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NOTES 

1.  Both  nodes,  N1  and  N2,  must  have  been  defined  on  previous  node 
cards  or  generation  cards. 

2.  The  node  number  increment  may  be  either  negative  or  positive. 
Incrementing  starts  with  the  first  node  after  N1  (NFIRST). 

N1  and  N2  need  have  no  numerical  relation  to  the  generated 
nodes,  only  the  spatial  relation  is  significant. 

3.  A11  generated  nodes  are  evenly  spaced  along  the  generated  line. 
3.2.2. 7  Cable  Element  Cards  (615,  ElO.O,  15,  5X,  2E10.0,  15) 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-5 

I 

Element  Number 

6-10 

IT(I,1) 

Node  number  for  1st  end  of  element 

11-15 

IT(I,2) 

Node  number  for  2nd  end  of  element 

16-20 

MAT(I) 

Material  number  used  for  this  element 

21-25 

KOMP(I) 

Compression  Code 

0  =  no  compression  stiffness 

1  =  can  support  compression 
-n  =  means  get  element  preload  from  n^ 
generated  catenary,  then  set  to  zero 

26-30 

KSKP 

Element  Generation  Code  (default  »  1) 

31-40 

SIG(I) 

Element  preload  force,  +  tension  (F) 

(see  notes  for  options) 

41-45 

ISIGO 

Start  up  procedure  flag  -  on  1st  element  only 
(see  notes  and  Section  3.5) 

51-60 

DSO(I) 

Unstretched  length,  required  only  when  ISIGO  »  +2 

61-70 

EAFAC 

Preload  factor,  used  only  on  the  first  element  when 
ISIGO  <  0.  (Default  value  =  0.001) 

71-75 

MEDIUM(I) 

Fluid  medium  flag 

0  means  element  is  in  water 
1  means  element  is  in  air 


1.  Elements  must  be  1n  ascending  element  number  order  and  the  first 
and  last  ones  must  be  input.  If  any  other  elements  are  omitted 
they  will  be  generated  by  using  the  KSKP  obtained  from  the  element 
card  following  the  omitted  elements  to  increment  the  node  numbers 
given  for  the  element  preceding  the  omitted  ones.  All  other  element 
data  will  be  assumed  to  be  the  same  as  given  by  the  element  given 
previous  to  the  omitted  ones. 

2.  The  startup  procedure  flag  (ISI60)  has  four  basic  options: 

ISIGO  =  0 

The  preloads  are  presumed  to  represent  an  equilibrium  state.  The 
element  lengths  obtained  from  the  nodal  positions  are  used  with 
these  preloads  and  the  constitutive  relations  to  calculate  the 
unstretched  lengths.  DSO(I)  and  EAFAC  are  not  required. 

ISIGO  =  1  (or  -1) 

The  configuration  represented  by  the  nodal  positions  is  an 
unstretched  compatible  arrangement  of  cables.  The  preloads 
(SIG(I))  are  merely  estimates.  The  unstretched  lengths  are 
obtained  from  the  nodal  positions.  When  ISIGO  =  -1  the  preload 
estimates  are  obtained  by  EAFAC  times  the  initial  value  of  EA  in 
the  material  constitutive  relation.  In  this  case  SIG(I)  is  not 
input. 

ISIGO  =  2  (or  -2) 

The  configuration  represented  by  the  nodal  positions  is  a  compatible 
estimate  of  the  deformed  configuration.  The  values  of  the 
unstretched  lengths  are  given  on  the  element  cards.  Element  loads 
are  calculated  from  the  unstretched  lengths,  stretched  lengths 
represented  by  the  nodal  positions,  and  the  material  constitutive 
relations.  EAFAC  is  used  to  estimate  a  preload  only  when  the 
estimated  configuration  results  in  a  zero  preload  in  an  element. 

ISIGO  =  3  (or  -3) 

The  nodal  positions  are  correct  only  for  the  constrained  nodes 
(fixed  and  moved).  They  do  not  represent  a  compatible  configuration 
of  the  system.  Input  of  the  unstretched  lengths  is  required  and  the 
element  preloads  may  be  estimated  by  input  or  from  EAFAC  (i.e.,  with 
ISIGO  =  -3).  The  option  does  not  appear  to  work  and  its  use  is  not 
recotinended. 


3.  When  a  string  of  elements  lies  along  a  catenary  curve  that  has  been 
generated  by  a  node  generation  card  (Section  3. 2. 2. 6),  the  element 
load  for  each  element  can  be  obtained  from  the  catenary  generation 
by  setting  KOMP(I)  for  each  element  on  the  line  to  -n,  where  n  is 
the  number  of  the  generation  card.  The  element  load  will  be  taken 
to  be  the  catenary  tension  half  way  between  the  two  nodes  defining 
the  element.  KOMP(I)  will  be  set  to  zero  (no  compression  stiffness) 
by  the  program. 


3.2. 2. 8  Cable  Material  Data  (215,  3E10.0,  2I5,E10.0/(2.E10.0)) 
Two  or  more  cards  required  for  each  material. 


COLS. 


1-5 

6-10 


11-20 

21-30 


31-40 

41-45 


VARIABLE 

NAME 


I 


IDRG(I) 


DIAM(I) 

G3(I) 


CAMC(I) 

ME(I) 


MED(I) 

TENULT(I) 


DESCRIPTION 


Material  number  (10  max) 

Drag  Coefficient  Code 

j'O  and  INDRAG  f  0  means  drag  coefficients 
will  be  obtained  from  subroutine  DRAGCO. 
Otherwise  the  default  coefficients  of 
Section  3.6  will  be  used. 


Cable  diameter,  used  only  for  fluid  load 
calculations  (L) 


Weight  per  unit  length  (FL”  )  negative 
means  lighter  than  reference  fluid 
(see  cols.  46-50) 

Added  mass  coefficient  (default  =*  1.0) 


Number  of  points  in  tens ion/s train  curve 
(if  zero  an  exponent  form  is  assumed) 


Reference  medium  code  for  this  material  data. 


0  means  water 
1  means  air 


Ultimate  Tension  capability  of  this  material 
(See  note  3  for  usage) 


If  ME(I)  j^.O  reads  ME(I)  cards  in  the  order  of  increasing  tension 

1-10  TT(J,I)  value  of  tension 

11-20  STR(J,I)  0^^  value  of  strain 

If  ME(I)  =  0  reads  one  card 

1-10  n(l.I) 

11-20  STR(1,I) 

^2 

where  T  =  c 

1.  When  point  input  is  used  and  data  is  required  outside  of  the 
given  values,  the  slope  of  the  last  (first)  segment  is  assumed 
to  be  continued. 

2.  Input  must  be  in  increasing  tension  order.  Thus  materials  with 
compression  stiffness  must  begin  by  listing  the  largest  compressive 
load  first  and  progressing  to  the  largest  tensile  load. 

3.  The  ultimate  tension  input  is  used  only  by  the  frequency  domain 
dynamic  solution.  When  TENULT(I)  is  greater  than  zero  the  random 
response  estimates  for  tensions  will  be  factored  and  compared 

to  the  ultimate  tension.  See  Section  3.2.3.10  for  further 
discussion  of  the  factors  and  usage  of  TENULT(I). 


3. 2. 2. 9  Buoy/Anchor  Data  (215,  3E10.0,  2I5y  2E10.0,  215) 


One  or 

two  cards  for 

each  body  (30  max) 

COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-5 

IM(I) 

Node  where  body  is  located 

6-10 

IDRB(I) 

Drag  Coefficient  Code 

11-20 

ADM(I) 

/  0  ^  INDRAG  7^  0  means  drag  coefficients 
will  be  obtained  from  subroutine  DRAGCO. 

Otherwise  the  default  coefficients  of  Section 

3.6  will  be  used. 

Buoyancy  (negative  means  heavier  than  the  fluid)  (F) 

21-30 

OBU(I) 

Body  diameter  (L) 

31-40 

BLEN(I) 

Body  length  (L) 

41-45 

JSLP(l.I) 

0  for  sphere  or  lump 
>0  for  in-line  cylinder 

Code  A  Notes  for  use 

46-50 

0SLP(2,I) 

Code  B 

51-60 

BAMC(r) 

Added  mass  coefficient 

61-70 

SBAMP(I) 

Default  =  0.5  for  sphere 
=1.0  for  cylinder 

Amplitude  of  surface  motion  (L) 

71-75 

IBS(I) 

Absolute  value  is  the  time  function  number 

76-80 

used  in  subroutine  TVARY  when  buoy  is  on 
surface  in  time  domain  solution 

Negative  sign  signals  this  body  is  a  mooring  buoy 

Continuation  code,  read  extra  data  card  when 

this  is  non-zero. 


Extra  Data  (4E10,0,I5) 

Read  only  when  continuation  code  is  non-zero.  These  values  assumed  to 
be  zero  when  card  not  read. 


1-10 

BWNO(I) 

11-20 

BSCD(I) 

21-30 

BMOM(I) 

31-40 

RELFAC(I) 

41-45 

JANCR(I) 

Wind  drag  coefficient  for  surface  buoy 

(Cp 


-D  *  V  «-  > 


Surface  current  drag  coefficient 


(Cp  *  A^)  (L") 


Mooring  buoy  mass  moment  of  inertia  (FIT  ) 

Holding  factor  for  dynamic  surface/bottom  constraint 

,  .  ^  ,  (default  =  1.001 

Anchor  fixity  flag(not  used  on  buoys) 

0  means  fix  vertical  only  ,  „ 

1  .eans  fix  an  three  d.o.f. 


DSSM-2 


NOTES 
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1.  Lumped  bodies  (not  mooring  buoys)  may  be  either  an  approximate 
sphere  or  an  in-line  cylinder.  Cylinders  are  assumed  to  have 
faired  ends  and  their  orientation  is  obtained  by  averaging  the 
directions  of  the  elements  indicated  by  Code  A  and  Code  B.  If 
only  Code  A  is  given  the  slope  of  that  element  is  used. 

2.  Mooring  buoys  require  4  nodes  in  their  definition.  IM(I)  for  a 
mooring  buoy  give  the  node  number  which  defines  the  position  of  the 
buoy  center.  Node  IM(I)  +  1  must  be  reserved  for  buoy  angular 
response.  Code  A  and  Code  B  give  the  numbers  of  the  two  attachment 
nodes. 


3.  When  a  buoy  is  on  the  surface  the  values  on  the  extra  card  will  be 
used  to  calculate  buoy  loads.  If  a  time  domain  analysis  is  being 
performed  the  buoy  can  be  given  a  vertical  movement  defined  by 
SBAMP(I)  and  TVARY  whenever  it  is  on  the  surface. 

4.  Surface  buoys  are  constrained  in  the  vertical  direction  but  are 
free  to  move  laterally.  This  constraint  remains  in  force  until 
the  line  tensions  are  sufficient  to  submerge  the  buoy. 

5.  Anchors  have  all  components  constrained  on  the  bottom  until  the 
line  tensions  are  sufficient  to  support  the  anchor  weight. 

6.  No  more  than  ten  lines  can  be  attached  to  any  buoy/anchor 
(including  mooring  buoys). 


7.  Generation  of  nodes  on  a  catenary  may  cause  the  generation 
of  up  to  tv/o  dummy  anchors  per  catenary  when  the  line  is 
tangent  to  the  bottom  before  the  bottom  node  is  reached. 


RELFAC  represents  a  dynamic  release  factor.  It  is  the  ratio  of 
the  vertical  load  it  takes  to  overcome  the  surface/bottom 
constraint  and  the  weight  or  buoyant  force  of  the  body.  It  is 
used  only  in  DYN  subanalyses.  All  other  subanalyses  use  a  factor 
of  1 .0. 


Revised  SEPT  1973 
Revised  OCT  1978 


3.2.2.10  Ship's  Data  (215,  4E10.0/8E10.0/8E10.0) 
Three  cards  are  required  for  each  ship. 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1st  Card 

- 

1-5 

ISHIP(I) 

Node  where  ship  is  located 

6-10 

LSHP(I) 

Load  function  option 

-1  search  file  for  equivalent  ship 

0  use  analytical  load  functions 

n>0  use  nth  ship  from  loading  file 

11-20 

CPROP(I) 

Propeller  resistance  coefficient  (default  *1.0 

21-30 

CR(I) 

Longitudinal  resistance  coefficient  for  hull 
to  be  used  only  for  analytical  functions 
(default  =  program  calculates  one) 

31-40 

CS(I) 

Hull  wetted  surface  coefficient  to  be  used  only 
for  analytical  functions  (default  =  2.70) 

41-50 

CMS 

Amidships  coefficient  to  be  used  to  calculate  Cn 
for  analytical  functions  ^ 
(default  *  0.98) 

2nd  Card 

1-10 

SLT(I) 

Total  length  of  ship  (L) 

11-20 

SAE(I) 

End  projected  wind  area  (L^) 

21-30 

SAS(I) 

Side  projected  wind  area  (L^) 

31-40 

SLWL(I) 

Water  line  length  (L) 

41-50 

SBEAM(I) 

Beam  at  midships  (L) 

51-60 

SDRFT(I) 

Draft  at  midships  (L) 

61-70 

SDSP(I) 

Volume  Displacement  (L^) 

71-80 

APROP(I) 

Propeller  projected  area  (L^) 
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COLS. 


VARIABLE 

NAME 


DESCRIPTION 


3rd  Card 


1-10 

FSFRWd) 

Load  table  wind  force  conversion  factor 
(default  =  1.0) 

11-20 

FSFRC(I) 

Load  table  current  force  conversion  factor 
(default  =  1.0) 

21-30 

FSLEN(I) 

Load  table  length  conversion  factor 
(default  =  1.0) 

31-40 

FSVEL(I) 

Load  table  velocity  conversion  factor 
(default  =  1.0) 

41-50 

SHIPK(1,I) 

22 

Heave  restoring  coefficient  (default  =  10  ) 

51-60 

SHIPK(2,I) 

22 

Roll  restoring  coefficient  (default  =  10  ) 

61-70 

SHIPK(3,I) 

22 

Pitch  restoring  coefficient  (default  =  10  ) 

71-80 

SHIPK(4,I) 

Heave/Pitch  restoring  coefficient  (default  =  0.) 

1.  Linearized  ship  restoring  coefficients  may  be  input  if  desired. 
Otherwise  the  ship  will  be  assumed  fixed  in  heave,  roll,  and 
pitch  during  static  analyses.  During  frequency  domain  analyses 
the  restoring  matrix  is  obtained  from  the  ship  motion  file. 

2.  The  second  card  may  be  left  blank  If  LSHP(I}  =  n  and  the  moored 
ship  is  the  same  as  on  file  n  of  the  ship  loading  file  (i.e., 
no  similarity  scaling  required). 

3.  The  conversions  factors  on  the  third  card  multiply  the  values 
of  force,  length,  and  velocity  from  the  ship  loading  file  to 
get  values  consistent  with  the  units  implied  by  the  rest  of  the 
input  data. 


When  analytical  loading  functions  are  used  the  longitudinal 
resistance  coefficient  may  be  Input  or  calculated  using  a 
table  look  up.  The  total  coefficient,  Cj^,  consists  of  three 
parts: 


residuary  coefficient 
frictional  resistance  coefficient 


=  fouling  resistance  coefficient 
The  value  for  Is  always  calculated  In  the  program  using 


-  a  0.456  1700 

(logioRe)  e 


where  is  Reynold's  Number. 

When  CR(I)  is  Input  it  is  taken  to  be  +  aC^.  When  It  Is 
zero  or  no  Input  Is  given,  then  1s  obtained  from  a  table 
and  aC^  Is  assumed  to  be  O.OOOS. 


This  set  of  cards  is  used  to  specify  the  analysis  activities,  select 
the  solution  routines  and  provide  the  required  input  parameters.  Various 
sequences  of  subanalyses  are  possible.  The  general  rule  is  that  any  subanalysis 
can  be  performed  if  an  appropriate  reference  state  is  defined  at  the  time  the 
analysis  is  requested.  Subanalyses  are  executed  in  a  sequential  manner  using 
the  state  from  input,  preceding  subanalysis,  or  RESTART  as  a  starting  point.  The 
various  options  are: 

SUBANALYSIS  TYPE  DESCRIPTION 

DEAD  Static  analysis  with  gravity  effects  only 

(possibly  including  point  loads) 

LIVE  Static  analysis  with  currents,  wind,  point 

loads,  and  gravity. 

DYN  Time  domain  dynamic  response  to  time  variable 

currents,  point  loads  and/or  boundary  motion 

MODE  Obtain  undamped  natural  frequencies  and  mode 

shapes. 

FREQ  Get  steady-state  dynamic  response  to  regular 

waves  and  superimpose  to  get  random  response 
estimates  using  a  wave  spectrum.  (Limited 
to  responses  driven  by  surface  ships  and/or 
mooring  buoys.) 

CHEK  Check  the  adequacy  of  lines,  buoys  and  anchors. 

NEW  Terminate  subanalysis  sequence  and  define  a 

new  analysis  case,  i.e.,  an  analysis  title 
card  is  read  next. 

END  Terminate  run. 

When  the  NEW  option  is  selected  a  completely  new  case  can  be  input  or 
the  RESTART  option  Initiated  to  return  to  a  previous  configuration  and  start  a 
new  subanalysis  sequence. 

A  variety  of  solution  procedures  are  provided.  They  can  be  classed  as 
linearized  Incremental,  self-correcting  incremental  or  iterative/incremental .  The 
following  is  a  list  of  the  solution  options: 
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OPTION  NO. 
0 


1 


2 


(3.4) 

5 


10 


DESCRIPTION 


Iterative/incremental  method  using  a  Modified 
Newton-Raphson  (MNR)  technique.  Iteratively 
solves  the  total  equations  at  each  load/time 
step.  For  dynamic  analyses  the  RFB  predictor 
is  used. 

Purely  incremental  procedure  using  a  Sequence 
of  Linear  Increments  (SLI).  For  dynamic 
analyses  the  linearized  incremental  motion 
equations  are  solved  iteratively. 

A  self-correcting  incremental  procedure  using 
the  Residual  Feedback  (RFB)  technique.  Dynamic 
solutions  obtained  with  implicit  algebraic 
solution. 

Options  Deactivated  —  do  not  use. 

Iterative  solution  of  total  equations  using  the 
Direct  Iteration  (DI)  Technique.  Applies  only 
to  DYNAMIC  LOAD  analyses.  The  approach  is 
essentially  a  single  step  predictor-corrector 
solution.  When  the  option  is  given  a  negative 
sign  the  RFB  method  is  used  as  a  predictor. 

A  special  starting  procedure  for  static  analyses. 
The  first  load  step  is  applied  in  a  specified 
number  of  sub-steps  using  the  RFB  method  then 
each  of  the  following  steps  using  the  MNR 
method.  In  dynamic  analyses  this  gives  a  MNR 
solution  with  an  explicit  predictor. 


Each  of  the  dynamic  time  domain  analyses  assumes  some  form  of  the  Newmark 
difference  equations.  Option  5  is  limited  to  dynamic  analyses  and  option  10  has 
special  limitations  as  noted.  Each  of  the  others  can  be  used  for  both  static 
and  dynamic  problems.  Moving  boundary  (moved  nodes)  dynamic  problems  cannot  be 
solved  using  Options  0,  10  or  2.  The  most  efficient  solution  method  for 
dynamics  appears  to  be  option  5,  while  for  static  problems  some  combination  of 
options  0  and  2  (or  perhaps  10)  should  be  considered.  Option  1  usually  requires 
very  small  increments  for  accuracy. 

Nonlinear  analyses  are  much  more  complicated  than  linear  ones.  One 
of  the  complicating  factors  is  the  possibility  of  incrementally  updating  the 
frame  of  reference  for  the  governing  equations  involving  large  displacements. 


> 
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SEADYN/DSSM  operates  with  three  configurations  of  the  structural  system  in  order 
to  solve  the  equations.  They  are: 

p 

Reference  Configuration  (  C):  The  configuration  used  to  describe  the 
element  positions  and  coordinate  transformations.  The  stresses,  strains 
and  displacements  are  described  relative  to  this  configuration.  It  may 
be  changed  periodically  in  the  analysis. 

Deformed  Configuration  (^C):  The  configuration  at  some  specified  time 
or  load  level.  This  is  the  initial  configuration  in  an  incremental 
procedure.  The  stiffness  matrix  is  written  for  this  configuration 

p 

relative  to  C. 

Incremental  Configuration  The  configuration  at  the  end  of  a 

time  or  load  step. 

A  fourth  configuration  which  represents  each  element  in  its  unstrained  state  is 
implied  but  not  identified  as  a  specific  configuration  of  the  system  since  it  may 
not  be  possible  to  satisfy  nodal  compatibility  with  the  unstrained  lengths. 


R  t 

Various  combinations  of  updates  of  C  and  C  can  lead  to  various  solution 

R  t 

procedures.  For  example,  if  both  the  configurations  C  and  C  were  held  at  the 
starting  configuration  the  incremental  procedures  reduce  to  a  simple  small  displace¬ 
ment  linear  analysis.  The  usual  nonlinear  incremental  procedure  is  obtained  by 
updating  (moving)  after  each  step.  This  updating  involves  the  regeneration  of 

P 

the  incremental  stiffness  matrix.  With  C  held  fixed  the  procedures  are  called 
Total  Lagrangian.  In  this  case  the  Local -to-Global  coordinate  transformations 

p 

remain  unchanged  throughout  the  analysis.  When  C  is  updated  these  transformations 

P 

must  also  be  updated.  Updating  or  moving  C  at  each  step  leads  to  a  less 
complicated  form  for  the  stiffness  matrix  but  adds  the  requirement  of  recalculating 
the  element  transformations  at  each  step.  Little  difference  in  the  results  should 

D 

be  expected  from  different  schemes  of  updating  C.  The  major  differences  come 
from  different  round-off  accumulations.  In  the  case  of  large  strain  of  highly 
nonlinear  materials  there  will  be  some  differences  due  to  the  approximation  of 
the  Eulerian  (updates  at  each  step)  incremental  constitutive  relation  using  the 
transformed  Lagrangian  form. 


Each  subanalysis  Is  initiated  with  a  Subanalysis  Option  (SAO)  Card. 

The  format  of  this  card  is  the  same  for  each  subanalysis  but  various  interpretations 
are  made  of  the  data  on  the  card  depending  on  the  subanalysis  type.  Various 
additional  cards  may  be  required  as  the  result  of  the  entries  on  the  SAO  Card. 

More  specific  discussion  of  the  additional  cards  is  given  following  the  description 
of  the  SAO  card. 


3. 2. 3.1  Subanalysis  Option  Card  (A4,  IX,  ISIS) 


VARIABLE 
COLS.  NAME 


DESCRIPTION 


1-4  TYPE  Analysis  Type,  DEAD,  LIVE,  etc. 

6-10  INP(1}  DEAD,  LIVE:  estimated  number  of  load  steps 

DYN :  not  used 

MODE:  mode  shape  order  flag 

0-11  St  mode  shapes  in  order  of 
increasing  frequency 
1-list  in  order  calculated 
FREQ:  number  of  load  steps  to  be  used  on  the 
first  of  the  drift  force  updates, 
(default  *  10) 

CHEK:  dynamic  configuration  record  number  to  be 
checked  (0  if  previous  SAO  not  FREQ) 

11-1S  INP(2)  DEAD,  LIVE,  DYN,  (FREQ):  number  of  steps 

between  printing 
MODE:  mode  shape  output  flag 
0-print  all  mode  shapes 
n-print  n  mode  shapes 
-n-print  n  mode  shores  and  also  write 
them  on  logical  unit  20. 

CHECK:  component  inventory  flag 
0-inventory  not  used 
1 -otherwise 


16-20  INP(3)  DEAD,  LIVE,  DYN:  Solution  option  number  (see 

section  3.2.3) 

FREQ:  Solution  option  for  drift  force  updates 
(only  0  or  2  allowed) 

(See  Note  3) 

CHEK:  component  inventory  print  flag 
0-do  not  print  inventory 
1 -otherwise 
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INP(4) 


DEAD.LIVE.DYN,  MODE:  Update  Option  flag 
establish  at  beginning  and 
0-do  not  update  (Total  Lagrangian) 

-update  •'C  at  step  Intervals  indicated 
-1-use  *^0  from  preceding  subanalysis  and 
do  not  update. 

FREQ:  iteration  option  flag 

0-start  from  original  static  reference 
1-start  from  last  configuration  found 

INP(5)  DEAD,  LIVE:  Start  up  option 

n>l-divide  first  step  into  n  subintervals 
(see  section  3.5) 

DYN:  not  used 

FREQ:  drift  force  iteration  flag 

0-do  not  update  configuration  for  drift  forces 
1-iterate  configurations  and  regular  wave 
solution 

-1-update  configuration  once 
+2 -same  as  +1  except  drop  the  sway  and  yaw 
forces  and  use  only  the  surge  force 
(for  single  point  moors) 

INP(6)  DEAD, LIVE, DYN:  number  of  point  loads  to  be  read. 

negative  means  use  the  loads  from  the 
preceding  subanalysis. 

FREQ:  Component  check  flag 

><0-form  dynamic  response  file  for  CHEK 
on  logical  unit  04 

INP(7)  LIVE, DYN:  fluid  load  option 
0-no  current 

-2-read  uniform  current  card 
-1-get  fluid  velocity  only  once  from 
subroutine  CURREN 

>^l-get  fluid  velocity  from  CURREN  each 
time  loads  are  calculated. 

FREQ:  number  of  wave  headings  (default  »  1) 

INP(8)  OYN:  number  of  payout/reel-in  ends 
FREQ:  mass  matrix  flag 
0-lumped  mass 
1 -consistent  mass 

INP(9)  DYN:  dynamic  initial  condition  flag 

0-static  start  or  continuation  of  a  dynamic 

analysis,  leave  nodal  velocities  as  is 
1-read  initial  velocity  card  and  set 
all  nodes  to  the  velocities  read 
FREQ:  free  ship  flag 

^0  means  calculate  unrestrained  ship  response 

INP(IO)  Optional  debug  output  flag 

(IBG)  (See  Section  4.2) 


-79- 


56-60  INP(n)  DEAD,  LIVE,  DYN:  restart  file  flag 

O-do  not  save  file 

-n-rewind  and  save  every  nth  output  record 
n^l -start  with  present  file  position  and 
save  every  nth  output  record. 

61-65  INP(12)  DEAD,  LIVE,  DYN(FREQ):  Step  Size  Control  Number, 

step  size  (options  0  and  5}  Is  Increased 
when  Iterations  converge  with  this  number 
of  iterations  or  less  (default  *  2) 

66-70  INP(13)  DEAD,  LIVE,  DYN(FREQ):  Iteration  Limit, step  size 

(options  0  and  5)  Is  decreased  If  con¬ 
vergence  not  obtained  In  this  number  of 
iterations  (default  “  20) 

71-75  INP(14)  DEAD,  LIVE,  DYN(FREQ):  Number  of  Trials, 

subanalysis  Is  aborted  if  step  size 
is  decreased  this  number  of  times 
without  convergence  (default  *  3) 

76-80  INP(15)  Newton-Raphson  Update  Interval.  Used  on 

solution  option  0  to  signal  how  often  to 
recalculate  the  tangent  stiffness  matrix. 

This  is  a  packed  word: 

•  Units  and  tens  digits  give  number  of  steps 

^  before  new  K^. 

•  Hundreds  and  above  digits  give  the  number 
of  alternating  sign  trials  on  a  given 
step  before  new  (default  =  205,  I.e., 

.  2  alternating  tries  and  5  steps) 


NOTES 

1.  If  TYPE  is  "NEW"  the  next  card  read  is  a  new  Analysis  Title  Card. 
If  TYPE  is  "END"  the  job  is  terminated.  Each  of  the  other 
subanalysis  options  are  listed  below  with  a  summary  of  the  cards 
that  must  or  may  be  read  in  association  with  the  option.  The 
additional  cards  are  listed  in  the  table  in  the  order  they  must 
appear  in  the  deck  if  they  are  used. 


-80- 


SUBANALYSIS  OPTION 


CARD 

SECTION 

DEAD 

LIVE 

DYN 

MODE 

FREQ 

CHEK 

Solution  Parameter 

3.2. 3.2 

R 

R 

R 

N 

R 

N 

Point  Loads  and 
Variation  Codes 

3.2. 3. 3 

0 

0 

0 

N 

N 

N 

Fluid  Velocity 

3. 2. 3.4 

N 

0 

0 

N 

N 

N 

Wind  and  Surface 

Current  Data 

3.2. 3.5 

N 

0 

N 

N 

N 

N 

Time  Data 

3.2. 3.6 

N 

N 

R 

N 

N 

N 

Moving  Boundary 

3.2. 3. 7 

N 

N 

0 

N 

N 

N 

Payout/Reel -in  Data 

3.2. 3.8 

N 

N 

0 

N 

N 

N 

Initial  Velocity 

3.2. 3.9 

N 

N 

0 

N 

N 

N 

Frequency  Domain  Data 

3.2.3.10 

N 

N 

N 

N 

R 

N 

Component  Check  Data 

3.2.3.11 

N 

N 

N 

N 

N 

R 

R  =  Required 
0  =  Optional 
N  =  Not  used 
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2.  Dynamic  time  domain  analysis  involving  surface  ships  and  mooring 
buoys  is  very  restricted.  The  program  does  not  generate  coefficients 
for  the  dynamic  response  of  these  bodies,  it  merely  assumes  they  are 
rigid  connectors  of  the  mooring  lines  and  a  place  where  external 
loads  can  be  applied.  The  mass  for  mooring  buoys  is  included  but 

no  mass  matrix  for  ships  is  generated.  Loading  from  wind  and  surface 
currents  can  only  be  defined  in  a  LIVE  subanalysis.  Although  the 
program  will  act  on  a  request  for  time  domain  responses  with  ships 
and  mooring  buoys  this  option  should  be  regarded  as  non-opera tional 
at  present. 

3.  The  solution  for  frequency  donain  responses  to  waves  (FREQ)  has 
built  into  it  a  LIVE  analysis  with  the  MNR  solution  (option  0)  or 
with  the  RFB  solution  (Option  2).  It  allows  the  choice  of  solution 
method  with  user-selected  step  size  on  the  first  configuration 
update  at  each  wave  heading.  When  iterative  updates  are  requested 
(INP(5)  >  0)  the  value  given  for  INP(4)  may  alter  the  solution 
option  used  on  updates  following  the  first  one.  If  INP(4)  ■  0 

the  same  procedure  is  used  on  each  update  with  the  initial  static 
reference.  If  INP(4)  =  1  then  the  configuration  from  the  previous 
update  is  used  as  a  starting  point  with  the  MNR  solution  in  one 
load  step.  Whenever  the  MNR  solution  is  used  parameters  in  INP(12) 
through  INP(15)  are  utilized  as  well  as  the  convergence  parameters 
from  the  Solution  Parameter  Card.  This  option  presumes  the  ship 
and  moor  have  been  moved  into  a  quasi -static  reference  configuration 
by  the  preceding  subanalysis  and  the  loads  (point  loads,  currents, 
wind,  gravity)  from  that  solution  remain  active.  No  load  related 
data  should  be  given  on  a  FREQ  SAG  card.  The  drift  force  updates 
presume  the  drift  forces  are  in  load  set  3;  therefore,  the  previous 
static  solution  cannot  use  load  set  3  if  drift  force  iterations  are 
requested. 

4.  The  FREQ  option  presumes  there  is  only  one  ship  and  the  motion 
coefficients  are  obtained  from  the  Ship  Motion  File. 


3. 2. 3. 2  Solution  Parameter  Card  (6E10.0) 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-10 

DMU 

Numerical  damping  factor  (see  note  1) 

11-20 

RERR 

Residual  norm  error  bound,  used  in  the  MNR 
solution  only  (see  note  2)  (default  ®  .001) 

21-30 

DERR 

Displacement  norm  error  bound,  (see  note  2) 
(default  »  .001) 

31-40 

DALPHA 

Internal  damping  multiplier  of  mass  matrix 
(FREQ  only) 

41-50 

DBETA 

Internal  damping  multiplier  of  stiffness 
matrix  (FREQ  only) 

51-60 

SRCHFC 

1-D  search  factor 

0.0  no  1-D  search  on  alternating  estimates 

<0.0  no  alternating  estimate  check 
>0.0  the  1-D  search  initial  guess  factor 

61-70 

PARMT 

MNR  Extrapolation  Parameter  (default  *  0.5) 

NOTES 

1.  The  numerical  damping  factor  is  used  to  avoid  problems  with 
singular  or  ill-conditioned  stiffness  matrices.  It  has  no 
influence  on  option  5.  It  should  be  used  with  caution  with 
options  1  and  2  since  it  alters  the  stiffness  matrix,  hence 
the  ec.uilibrium  equations,  in  these  two  cases.  Possible 
values  are 

DMU  ^1.0  very  heavy 

1.0  2  DMU  ^0.1  heavy 

0.1  i  DMU  i  0.01  moderate 

0.01  ^  DMU  >0.001  light 

0.001  ^  DMU  very  light 

-8 

If  DMU  <  10  and  a  singularity  is  encountered  then  DMU  is  set  to 
0.001  and  the  step  is  tried  again.  On  a  repeated  singularity 
DMU  is  increased  to  0.1  and  the  step  is  tried  once  more.  If  the 
singularity  persists  the  calculation  is  aborted. 
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The  error  bounds  are  used  to  test  for  convergence  of  the  iterative 
solutions.  REftR  is  used  only  in  the  MNR  solutions  (option  0).  DERR 
is  used  in  options  0  and  5.  The  conWrgence  criteria  is 

RNORM  £  /RERR/ 

DNORM  <  /DERR/ 

Where  IrnORM  and  DNORM  represent  norms  of  the  nodal  force  residual 
and  displacement  increments,  respectively.  Both  of  these  inequalities 
must  be  satisfied  for  convergence  of  the  MNR  iterations  while 
only  the  displacements  are  checked  in  the  DI  method  (option  5). 

Some  flexibility  in  the  form  of  these  norms  is  available.  If  the 
value  of  RERR  is  input  greater  than  zero,  then 

N  1/2 

RNORM  =(_^r//N)  /T^^^  ;  ,  - 

where 

R.  =  the  i^^  component  of  the  nodal  residual 

N  =  the  total  number  of  nodal  degrees  of  freedom 

is  the  maximum  element  tension 
max 

If  the  value  of  RERR  is  negative,  then  .  • 


RNORM  = 

where 

M  =  the  total  number  of  nodal  components  including  the  fixed 
nodes,  i.e.,  the  reactions  are  included  in  the  denominator. 
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If  the  value  of  DERR  is  greater  than  zero,  then 
N‘ 

DNORM  =  (  A(Aq.f/N) 

where 


^A(Aq^.)  =  the  component  of  the  change  in  the  displacement 
increment. 

If  DERR  is  negative,  then 

.  N  0,1/2 


DNORM  = 


I  A(Aq  f 
i=l  ' 

N  2 


See  the  discussion  of  norms  in  Section  2. 

3.  The  terms  DALPHA  and  DBETA  can  be  used  to  specify  internal  damping 
proportional  to  the  mass  and/or  stiffness  matrix.  The  assumed  ( 
form  of  the  damping  matrix  is: 

[C]  =  a[M]  +  6[K] 

\ 

At  the  present  this  is  available  only  on  the  frequency  domain 
solution. 

4.  The  1-D  search  factor,  SRCHFC,  can  be  used  with  the  MNR  solution 
(option  0)  in  an  attempt  to  enhance  a  poor  ini tiaT con-figuration. 

See  the  discussion  of  the  MNR  method  in  Section  2.  . 

5.  PARMT  is  an  extrapolation  param.eter  used  in  incremental  MNR  solutions. 
The  starting  estimate  for  the  displacement" on  all  but  the  first  step 
is  the  accumulated  displacements  from  the  previous  step  plus  PARMT 
times  the  change  in  displacement  calculated  from  the  preceding  step. 
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MICROCOPY  RESOLUTION  TEST  CHART 
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3.2. 3. 3  Point  Loads 


Expects  INP(6)  load  cards  plus  a  variation  code  card 
Load  Cards  (215,  3E10.0) 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-5 

I 

Node  where  load  Is  applied 

6-10 

J 

Load  set  number  (3  max,  default 

11-20 

FP(3*I-2,J) 

X  component  of  load 

21-30 

FP(3*I-1.J) 

y  component  of  load 

31-40 

FP(3*I,J) 

z  component  of  load 

NOTE 

If  the  node  I  Is  used  for  the  angular  position  of  a  ship  or  mooring 
buoy  the  load  1s  assumed  to  be  a  moment  about  the  specified  axis  with  units  of 
(FL).  Otherwise  the  loads  are  point  forces  with  units  of  (F). 

Variation  Code  Card  (315) 


1-5 

ILF(l) 

Load  set  variation 

code 

for  set  1 

6-10 

ILF(2) 

Load  set  variation 

code 

for  set  2 

11-15 

ILF(3) 

Load  set  variation 

code 

for  set  3 

NOTE 

For  DEAD  or  LIVE 

The  variation  codes  signal  the  variation  form  during  the  subanalysis. 
1  -  Increment  the  load  set  from  0  to  full  value 
0  -  hold  load  set  at  full  value 
-1  -  decrease  the  load  set  from  full  value  to  zero. 

FOR  DYN 

0  hold  load  set  constant  at  full  value 

fO  use  the  code  in  subroutine  TVARY  to  obtain  the  time  variable 
factor  (see  Section  3.3) 
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3.2. 3.4  Fluid  Velocity  Card  (3E10.0) 

One  card  used  only  when  INP(7)  ■-2 

VARTARI  r 

COLS.  NAME  DESCRIPTION 


1-10 

VI(1) 

11-20 

VI(2) 

Vy  velocity  components  of  uniform  flow  field 

21-30 

VI(3) 

Vz 

Wind  and  Surface  Current  Data  (6E10.0} 

One  card  used  only  when  NSHIPS  f  0  on  LIVE  subanalysis. 

1-10 

WIND 

Wind  velocity  (LT“^) 

n-20 

WAD 

Initial  global  heading  of  wind  (degrees) 

21-30 

CURNT 

Surface  Current  Velocity  (LT“^) 

31-40 

CAD 

Initial  global  heading  of  current  (degrees) 

41-50 

HEDINC 

Heading  Increment  (degrees) 

51-60 

HEDEND 

Final  wind  heading  (degrees) 

NOTES 

1.  See  Note  6  of  Section  3. 2. 2.5  for  global  heading  convention. 

2.  When  HEDINC  is  not  zero  a  sequence  of  static  configurations 
will  be  generated  until  the  wind  heading  moves  from  WAD  to 

>  HEDEND.  The  difference  between  WAD  and  CAD  is  held  fixed. 
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3.2. 3.6  Time  Data  Card  (6E10.0) 

One  card  required  for  each  DYN  subanalysis. 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-10 

DT 

Time  Increment 

11-20 

DTU 

Update  time  Increment 

(used  to  signal  the  update  of  ^C 

see  Section  3.2.3) 

21-30 

TMAX 

Limit  time 

31-40 

GAMNEW 

Gamma  Integration  parameter 

41-50 

BETNEW 

Beta  Integration  parameter 

51-60 

ALPNEW 

Alpha  Integration  parameter 

NOTES 

1.  If  DTU  »  0  or  blank  OTU  Is  set  *  TMAX,  I.e.,  no  update  of  ^C. 

2.  If  DT  _<  0  the  program  will  estimate  the  time  step.  If  the 
Input  value  of  DT  =  -A  then  the  estimated  value  Is  multiplied 
by  A.  If  DTU  =  -B  then  OTU  is  set  to  B  times  the  DT  estimate. 

3.  The  integration  parameters  are  those  from  the  generalized  Newmark 
difference  equations  (see  Section  2.2.2).  The  Alpha  parameter 
is  not  used  in  the  MNR  solution  routine  (solution  option  0). 


■s*s**4 
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3. 2. 3. 7  Moving  Boundary  Data  Cards  (15,  5X,  3(215,  E10.0)} 

One  card  for  each  moved  node  (see  Master  Control  Card) 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-5 

I 

Moved  Node  No. 

11-15 

JB(3*(I-1)+1) 

X  motion  code  (see  Notes) 

16-20 

IMTMFN(.+1) 

X  time  function  number 

21-30 

UB(.+1) 

X  motion  amplitude 

31-35 

JB(.+2) 

y  motion  code 

36-40 

IMTMFN(.+2) 

y  time  function  number 

41-50 

UB(.+2) 

y  motion  amplitude 

51-55 

JB(.+3) 

z  motion  code 

56-60 

IMTMFN(.+3) 

z  time  function  number 

61-70 

UB(.+3) 

z  motion  amplitude 

NOTES 

1.  The  motion  codes 

1  =  displacements  are  specified 

2  =  velocities  are  specified 

3  *  accelerations  are  specified 

2.  Non-zero  time  function  numbers  are  passed  to  the  TVARY  subroutine 
to  obtain  the  time  variable  factor  (see  Section  3.3).  When  the 
time  function  number  is  zero,  the  motion  component  Is  held  fixed 
at  the  value  given. 


3.2. 3. 8  Payout/Ree1-1n  Data  (2I5.2E10.0.5I5} 

One  card  for  each  payout/reel-1n  end  (5  max) 


VARIABLE 

NAME 

DESCRIPTION 

1-5 

JOP(I) 

Node  where  payout/reel -In  occurs 

6-10 

JPELT(I) 

Initial  element  number  for  payout/reel  ' 

n-20 

PAYV(I) 

Payout  rate(+  payout,  -  reel-ln) 

21-30 

AMAXL(I) 

Mitosis  Length 

31-35 

NGROW(I) 

Number  of  nodes  available  for  growth 

36-40 

NSHRINK(I) 

Number  of  nodes  that  can  be  reeled  In 

41-45 

NNPOI(T) 

Node  number  Increment  (default  »  1) 

46-50 

NELPOI(I) 

Element  nimiber  Increment  (default  *  1) 

51-55 

NPOVRY(I) 

Number  for  time  variation  function 

0  means  constant  rate 

means  get  time  variation 
from  TVARY  subroutine. 

NOTES 

1.  Payout/reel -In  option  not  provided  under  Contract  N62477-76-C-0002. 
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3. 2. 3. 9  Initial  Velocity  Card  (3E10.0) 
One  card  if  INP(9)  ■  1. 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-10 

VIO) 

11-20 

VI(2) 

''y 

components  of  uniform  initial  velocity 

21-30 

VI(3) 

V, 

NOTES 

1.  This  card  is  used  in  systems  such  as  towed  bodies  where  the  initial 
state  is  a  quasi-static  condition  where  all  nodes  in  the  system 
are  moving  with  the  same  velocity  without  accelerating.  This 
option  assumes  the  static  configuration  is  defined  by  input  or 
was  obtained  in  a  previous  LIVE  analysis  with  a  uniform  flow 
field  having  the  same  magnitude  as  that  given  here  but  the 
opposite  sense. 


3.2.3.10  Frequency  Domain  Data  Set 

This  set  of  data  cards  is  read  only  for  a  FREQ  subanalysis.  The  data 
set  consists  of  the  following; 

SPECTRUM  CARD 

SHIP  MOTION  FILE  CONVERSION  FACTORS  CARD 

repeat  the  following  for  each  heading  requested  on  FREQ  SAO  card 

WAVE  HEADING  CARD 

CALCULATION  OPTION  CARD 

REGULAR  WAVE  RESPONSE  CARDS  (Optional) 

RANDOM  RESPONSE  CARDS  (Optional) 

The  Ship  Motion  File  is  assumed  to  be  available  on  logical  unit  08 
whenever  there  is  a  ship  in  the  system.  There  must  be  a  ship  or  at  least  one 
surface  mooring  buoy. 


SPECTRUM  CARD  (6E10.0) 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-10 

SPECA 

Spectrum  coefficient  A,  (L^  T'^) 

11-20 

SPECB 

Spectrum  coefficient  B,  (T"^) 

21-30 

DOMG 

Aw,  frequency  increment  (T~^) 

31-40 

OMGMN 

“min’  bound  on  frequency  (T”^) 

41-50 

OMGMX 

uppe>^  bound  on  frequency  (T’^) 

51-60 

AMPMN 

Cut-off  amplitude  for  waves,  (L) 

(default  =  0.0001,  see  Note  2) 


NOTES 

1.  The  wave  spectrum  is  assumed  to  have  the  form 

r/  \  _  rt/  5  -B/u) 

5  ( w  J  -  A/  w  e 

where  S(u)  is  based  on  twice  the  square  of  the  wave  height.  Any 
spectrum  having  this  general  form  can  be  input.  Values  for 
common  spectra  are  listed  below: 
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SPECTRUM 

Piers  on-Moskowitz 

Bretschneider 

I.S.S.C. 


A 

135.0 

4200  H  ^/T  ^ 
s  s 

2760 


B 

9.7  X  loVV|^^ 

1050/Tj^ 

690/t/ 


where 

V|j  =  wind  speed  (knots) 

*  significant  wave  height  (ft) 

Tj  “  significant  wave  period  (sec) 

u  “  circular  frequency  (radians/sec) 

2 

S(uj)  =  spectral  value  (ft  -sec) 


2.  Regular  wave  responses  are  calculated  for  waves  with  frequencies 

between  u  and  u  which  have  wave  amplitudes  (based  on  the 
nt  1  n  max 

spectrum)  greater  than  AMPMN.  Incrementing  proceeds  beyond 
until  a  wave  amplitude  less  than  AMPMN  is  encountered. 

Ship  Motion  File  Conversion  Factors  (3E10.0) 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-10 

FRCFAC 

force  conversion  factor 

11-20 

FACLEN 

length  conversion  factor 

21-30 

TIMFAC 

time  conversion  factor 

NOTES 

1.  These  factors  are  used  as  multiplifers  of  the  data  on  the  Ship 
Motion  File.  Each  of  them  have  a  default  value  of  1.0. 

Wave  Heading  (ElO.O) 


1-10 


GHED 


Wave  heading  in  global  system  (degrees) 
(see  Note  6  of  Section  3.2. 2.5  for  global 
heading  convention.) 


VARIABLE 

NAME  DESCRIPTION 

CALTYP  "RAND"  for  random  response 

"REGb"  for  regular  wave  response 
"DONE"  to  terminate  option 

1.  The  word  given  signals  what  set  of  data  is  to  be  read  next.  When 
CALTYP  =  "REGb"  the  next  cards  read  are  requests  for  regular  wave 
responses  at  specific  nodes.  When  CALTYP  *  "RAND"  the  next  cards 
read  are  requests  for  random  response  calculations  for  various 
components. 

2.  These  option  sets  can  be  repeated  as  often  as  desired  for  a  given 
wave  heading.  The  sequence  is  terminated  with  a  DONE  Card. 

Regular  Wave  Response  Data  (215,  2E10.0,  215) 


1-5 

NODE 

Node  for  which  response  is  requested.  Zero 
means  no  more  data  and  a  new  calculation  option 
card  will  be  read  next. 

6-10 

LOCAL 

0-produce  output  in  global  system 

1-produce  output  in  ship's  local  system 

11-20 

WFRQ 

Circular  frequency  of  response  (T"^) 

21-30 

WAMP 

Wave  amplitude  for  response  (L)  (default  =  1.0) 

31-35 

lOCODE 

0-output  is  displacements 

1 -output  is  position 

36-40 

NUM 

Number  of  divisions  per  cycle  (default  =  30) 

1.  These  cards  will  generate  the  displacement  versus  time  or  position 
versus  time  for  all  three  nodal  components  through  one  cycle  of  moti 

2.  The  local  coordinate  system  will  not  be  used  if  position  output 
is  requested. 

3.  There  is  no  limit  to  the  number  of  requests  that  can  be  made.  The 
option  is  terminated  by  a  blank  card  (NODE  =  0). 


COLS. 

1-5 


Random  Response  Data  (A4,  IX,  215.  5X,  4E10.0) 

VARIABLE 

COLS. 

NAME 

DESCRIPTION 

1-4 

VTYPE 

"NODE".  "SHIP",  "TENS",  "DONE"  (see  notes) 

6-10 

NUMC 

Node  or  element  number 

11-15 

NDRT 

If  VTYPE  =  "NODE"  gives  the  global  component 
di recti on 

1  *  X 

2  =  y 

3=2 

21-30 

SFSTAT 

Static  Load  Factor  ("TENS"  only) 

31-40 

SF3 

(1/3)  Load  Factor  ("TENS"  only) 

41-50 

SFIO 

(1/10)  Load  Factor  ("TENS"  only) 

51-60 

SFIOO 

(1/100)  Load  Factor  ("TENS"  only) 

1.  The  spectral  response  of  the  ship,  element  tensions,  and  any  of 
the  nodal  displacement  components  can  be  estimated.  VTYPE  signals 
which  one  is  desired.  Any  number  of  these  cards  may  be  provided. 

The  input  is  terminated  by  a  card  with  VTYPE  *  "DONE".  The  response 
data  calculated  represent  the  average  of  the  1/3,  1/10,  1/100 
highest  responses. 

2.  VTYPE  =  "SHIP"  requires  no  other  entries  on  the  card. 

3.  VTYPE  =  "TENS"  requires  the  specification  of  the  element  number 
in  NUMC.  If  values  are  given  for  SFSTAT  through  SFIOO  and  if  the 
element  material  has  the  ultimate  tension  specified,  then  estimates 
of  the  factored  maximum  loads  will  be  printed. 


4.  VTYPE  =  "NODE"  requires  the  specification  of  the  node  number  in 
NUMC  and  the  component  direction  in  NDRT. 


5.  One  must  include  a  TENS  card  for  each  element  that  will  be  involved 
in  a  subsequent  component  adequacy  check.  For  example,  if  a  check 
of  an  anchor  capacity  is  to  be  made  including  dynamic  effects, 
every  element  which  connects  to  the  anchor  (or  fixed  node  where 
the  anchor  is)  must  be  called  out  on  a  TENS  Card. 

6.  The  wave  heading  card  and  the  associated  random  response  data 
cards  must  be  repeated  for  each  heading  requested  on  the  FREQ 
SAO  card  as  INP(l). 

3.2.3.11  Component  Check  Data  (A4,  IX,  315,  3E10.0) 

Read  only  for  TYPE  =  CHEK  on  SAO  card. 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-4 

CTYPE 

"ANCH",  "BUOY",  "LINE",  "DONE" 

6-10 

NELD 

Node  or  element  number 

11-15 

I  CODE 

Component  code 

16-20 

JCODE 

Bottom  factor  code 

21-30 

CAPY 

Component  capacity 

(default  =  use  inventory  with  SAFAC) 

31-40 

CMP  ID 

Component  Identifier  for  Inventory 
Anchor  Weight 

Chain  Size 

Line  Diameter 

Buoy  O.D. 

(input  in  inventory  units) 

41-50 

SAFAC 

Safety  Factor  (default  =  1.0) 

NOTES 

1.  CTYPE  identifies  the  type  of  component  to  be  checked.  Only 
anchors,  buoys,  or  lines  are  allowed.  Any  number  of  cards 
can  be  given.  Input  is  terminated  by  CTYPE  =  "DONE". 

2.  Negative  I  CODE  with  CTYPE  =  ANCH  means  the  vertical  load  will 
be  used  in  the  capacity  check. 


3.  The  CHEK  SAO  can  be  initiated  following  any  other  subanalysis'.  If 
the  previous  one  was  not  FREQ,  the  present  state  is  evaluated.  If 
it  was  FREQ  and  the  dynamic  response  file  was  generated,  the 
dynamic  tensions  will  be  included  in  the  check  for. those  element 
on  the  file.  (See  Note  5  of  Section  3.2.3.10) 

4.  If  CTYPE  is  "ANCH"  or  “BUOY"  then  NELD  is  the  node  where  the  anchor 
or  buoy  is  located.  Anchor  checks  can  be  made  for  fixed  nodes 

as  well  as  Tor  active  nodes  where  bodies  have  been  defined. 

5.  No  more  than  ten  lines  can  be  attached  to  any  anchor  or  buoy  at 
an  active  node  (see  Note  6,  Section  3. 2. 2. 9).  Anchor  checks  at 
fixed  nodes  allow  up  to  twenty  lines  to  be  attached  to  the  fixed 
node. 

6.  Component  Codes 

ANCHORS 

1.  Navy  standard  stockless  with  stabilizers 

2.  NAVSHIPS  lightweight 

3.  NAVFAC  stato 

4.  Imbedment  (no  inventory) 

5.  Stock  (Admiralty)  (no  inventory) 

6.  Mushroom  (no  inventory) 

BUOYS 

1.  Bar  riser  chain 

2.  Spherical  or  other  (no  inventory) 

LINES 

0.  Chain 

(R) 

1.  Samson  2-in-l  ^braided  nylon 

2.  Samson  2-in-l^  power  braid 

3.  Samson  2-in-l  stable  braid 

4.  Samson  Blue  Streak 

5.  Other  (no  inventory)  .  •  ■ 


7.  Bottom  factor  codes  for  anchors 


1 .  Compacted  Sand 

2.  Stiff  dense  clay 

3.  Sticky  clay  of  medium  density  (cohesive) 

4.  Soft  mud  (fluid),  loose  coarse  sand,  gravel 

5.  Hard  bottom  (rock,  shale,  boulders) 

3.2.4  Restart  Capability 

It  is  possible  to  save  the  contents  of  COMMON  at  each  output  interval 
for  use  as  restart  data.  If  a  save  of  the  data  is  requested  the  contents  of 
COMMON  (labeled  and  unlabeled)  are  written  on  tape  in  a  single  logical  record 
in  the  binary  mode. 

The  first  record  on  the  tape  is  the  title  card  from  the  analysis  which 
creates  the  tape.  The  first  word  in  this  title  may  be  used  as  an  identifier  for 
the  tape.  The  COMMON  restart  records  are  then  numbered  sequentially  after  the 
title  record. 

The  restart  files  can  be  saved  during  a  solution  sequence  by  entering 
a  value  for  IMP(ll)  on  the  SAO  Card.  The  appropriate  values  are: 

0  Do  not  save  file. 

-n  Rewind  File  and  Write  Title  REcord  -  Then  save  COMMON. 

n  Begin  saving  COMMON  from  present  position  on  the  file 
(Title  Record  is  written  if  the  file  counter  is  zero.) 

where  every  nth  output  record  will  be  saved. 

The  program  recognizes  the  following  file  codes  for  saving  the  data: 


FILE  CODE 

GENERATED  SY 

01 

Dead  Load  Analysis 

02 

Live  Load  Analysis 

03 

Dynamic  Analysis 
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The  restart  capability  has  two  functions.  First,  it  allows  the  restart 
of  the  job  on  a  subsequent  run  in  which  the  contents  of  the  file  are  to  be 
preserved.  This  is  accomplished  by  attaching  the  appropriate  file  and  providing 
the  RESTART  input  cards.  The  second  function  is  to  allow  the  user  to  return  to 
the  results  of  a  previous  calculation  (e.g.,  multiple  live  loads  from  the  same 
dead  load  configuration).  In  either  case  the  analysis  is  treated  as  a  new  data 
case  and  the  three  cards  given  below  are  required.  The  normal  data  sequence  can 
be  followed  beyond  the  point  of  the  restart. 

3. 2. 4.1  Title  Card  (12A6) 

VARIABLE 

COLS.  NAME  DESCRIPTION 

1-72  HED  Any  title  (see  below  for  special  usage) 

3. 2. 4. 2  Restart  Control  Card  (415) 

1-5  NFILE  Minus  the  number  of  the  record  to  be  read  from 

the  file  (see  output  heading) 

0  means  use  the  last  record  written  in  present  run 

6-10  NTAPE  File  Code  Number  (usually  01,  02,  03,  or  04) 

11-15  Identification  Check  Flag 

1  -  Read  the  title  record  from  the  file  and 
compare  1st  6  characters  with  the  1st 
6  characters  in  HED  (above).  Abort  if 
they  are  not  the  same 

0  -  No  I.O.  check 

16-20  IRST  Restart  Save  Flag 

(same  function  as  INP(ll)  on  SAO  Card) 


3. 2. 4. 3  Restart  Option  Card  (A4,  IX,  15,  6E10.0} 


COLS. 

VARIABLE 

NAME 

DESCRIPTION 

1-4 

Analysis  Type  (DEAD,  LIVE,  DYN,  blank) 

If  blank  -  this  signals  the  contents  of  the 
COMMON  record  are  to  be  the  initial 
conditions  for  a  new  load  case.  The  SAO 

Card  is  read  next. 

If  DEAD.  LIVE,  or  DYN  the  COMMON  record  must 
contain  data  consistent  with  the  type  given 
or  the  run  is  aborted. 

6-10 

NIPR 

New  value  for  number  of  steps  between  printing. 
If  this  is  >  0  and  not  equal  to  the  value 
given  previously  as  INP(2),  the  output  interval 
will  be  changed  by  the  ratio  NIPR/INP(2). 

11-20 

on 

New  Step  Size 

0  New  load  increment  for  DEAD  or  LIVE 

0  New  small  time  step  (DT)  if  DYN 

0  Blank  if  no  change 

21-30 

OTUT 

New  Update  Time 

31-40 

TMAXT 

New  Time  Limit 

41-50 

TGAM 

New  Value  for  Newmark's  y 

51-60 

TBET 

New  Value  for  Newmark's  b 

61-70 

TALP 

New  Value  for  Alpha 

NOTE:  Leave  columns  21-70  blank  if  not  DYN  or  DTT  =  0. 


: 


-101- 


USER  SUPPLIED  SUBROUTINES 


Three  subroutines  are  available  to  be  written  by  the  SEADYN/DSSM  user 
and  are  called  If  signalled  by  Input  data.  Each  of  them  are  described  below. 

3.3.1  Subroutine  DRAGCO 

This  subroutine  allows  the  user  to  specify  various  sets  of  drag  coefficients. 
Calls  to  DRAGCO  are  Initiated  when  INORAG  j*  0  on  the  Parameter  Card  (Section 
3. 2. 2. 2)  and  either  the  cable  material  data  (Section  3. 2. 2. 8)  or  the  lumped  body 
data  (Section  3. 2. 2. 9)  provide  a  non-zero  drag  coefficient  number.  The  calling 
sequence  for  the  subroutine  Is:  ^ 

SUBROUTINE  DRAGCO  (INDRAG, IBOD,IDR,RE,RET,CN,CT) 

where 

INDRAG  >  Option  call  number 

IBOO  =  1  for  spherical  buoy 

2  for  cylindrical  buoy 

3  for  cable 

IDR  »  drag  coefficient  number  specified  by  lumped  body  or  cable 
Input  data 

RE  =  Reynolds  number  based  on  the  normal  component  of  the 
relative  fluid  velocity 

RET  =  Reynolds  number  based  on  the  tangential  component  of  the 
relative  fluid  velocity  (not  given  on  spheres) 

CN  =  Return  variable  for  the  calculated  normal  drag  coefficient 

CT  =  Return  variable  for  the  calculated  tangential  drag 

coefficient  (not  used  on  spheres). 

INDRAG,  IBOD  and  IDR  can  be  used  as  indices  to  select  the  drag  coefficients  from 
user  defined  catalogs  of  functions  if  that  Is  desired. 
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3.3.2  Subroutine  CURREN 

This  subroutine  specifies  the  spatial  and  temporal  characteristics  of 
the  fluid  flow  field.  It  will  be  called  once  in  each  subanalysis  where 
INP(7)  »  -1  (Section  3.2. 3.1)  or  every  time  the  fluid  loads  are  evaluated  if 
INP(7)  >  0.  The  subroutine  provides  the  three  components  of  the  fluid  flow 
velocity  in  the  global  coordinate  system  for  every  node  on  each  call.  The 
calling  sequence  is: 

SUBROUTINE  CURREN  (T,X,N3,V,NFL) 

where 

T  =  Time  at  which  the  current  is  evaluated 

X  =  Nodal  position  vector 

X(3*(I-1)+J)  gives  the  Jth  component  at  the  Ith  node 

N3  =3  times  the  number  of  nodes 

V  =  The  calculated  flow  velocity  components 

V(3*(I-1)+J)  gives  the  Jth  component  at  the  Ith  node 

NFL  =  The  flow  field  number  (INP(7)  >  0) 

3.3.3  Subroutine  TVARY 

This  subroutine  specifies  the  time  variation  for  each  of  the  point  loads 
and/or  moved  nodes.  It  is  called  once  for  each  time  variable  component  at  T  =  0. 
and  at  each  time  step.  The  calling  sequence  is: 

SUBROUTINE  TVARY  (T.FNT.IFN) 

where 

T  =  time  for  the  functional  evaluation 
FNT  =  the  calculated  value  of  the  function 
IFN  =  time  function  number 
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3.4  PROBLEM  SIZE  LIMITS  AND  FILE  USAGE 

The  SEADYN/DSSM  computer  program  uses  Implied  dimensioning  of  the 
COMMON  storage  area  to  make  It  relatively  problem  size  Independent.  The  problem 
size  as  Indicated  by  the  total  number  of  nodes  and  elements  and  the  type  of 
analyses  requested  Is  checked  to  see  If  the  COMMON  block' Is  large  enough  to 
complete  the  analysis.  If  It  Is,  the  analysis  proceeds.  If  It  Is  not,  a 
message  Is  printed  and  the  rest  of  the  analysis  Is  skipped  by  searching  for  the 
next  NEW  or  END  SAO  Card. 

The  size  of  the  COMMON  block  can  be  set  In  the  main  program  which  Is 
simply  a  starter  routine  for  the  SEADYN/DSSM  program.  The  main  program  Is: 

COMMON  A(xxxxx) 

NCOM  “  xxxxx 
CALL  SEADYN  (NCOM) 

STOP 

END 

The  minimum  required  values  for  NCOM  for  the  various  subanalyses  can 
be  calculated  from  the  following  expression. 

NCOM  i  NBASE  +  NSUBAL 

where 

NBASE  =  66*NN  +  30*NE 

NN  =  Total  number  of  nodes 

NE  =  Total  number  of  cable  elements 

The  various  values  for  NSUBAL  are: 

SAO  NSUBAL 

DEAD,  LIVE,  DYN  NF3*IBEND 

MODE  2*(NF3+1)*NF3 

FREQ  2*NF3*IBEND 

CHEK  0 


NF3  =  3*(NN-NFIX-NSLAVE-MVB*IDYN) 

IBEND  »  half-bandwidth  of  stiffness  matrix 
(0  for  DYN  with  solution  option  5) 

NFIX  =  number  of  fixed  nodes 

NSLAVE  =  number  of  slaved  nodes 

MVB  «  number  of  moved  nodes 

I DYN  =  0  for  DYN  option 

1  for  DEAD  and  LIVE  options 

The  value  of  NF3  used  by  the  MODE  and  FREQ  options  is  the  NF3  from  the  immediately 
preceding  subanalysis. 

The  program  prints  two  values  for  the  half-bandwidth  when  the  element 
input  data  is  processed.  The  largest  value  corresponds  to  the  size  required  by 
the  SLI  solution  (option  1)  when  MVB  f  0  and  a  DYN  SAO  is  used.  The  smaller 
value  pertains  to  all  other  situations  using  a  stiffness  matrix.  The  two  values 
will  be  identical  when  MVB  *  0. 

The  program  does  place  size  limits  on  some  of  the  special  features  and 


operations.  These  limits  are  summarized  below: 

Maximum  number  of  material  types  =10 

Maximum  number  of  points  in  any  Tension/Strain  =  20 
Table 

Maximum  number  of  moving  boundary  nodes  =  5 

Maximum  number  of  payout/reel -in  points  =  5 

Maximum  number  of  ships  in  the  system 

DEAD, LIVE, DYN  =  3 

FREQ  =  1 

Maximum  number  of  anchor/buoys  =  30 

Maximum  number  of  lines  attached  to  an  active 

node  where  anchor/buoy  is  located  =  10 

Maximum  number  of  lines  connected  to  a  fixed 
node  for  CHEK  SAO  =  20 

Maximum  number  of  generated  catenaries  =  20 


These  limits  are  due  to  assumed  array  dimensions  and  can  be  altered  with  some 
minor  reprogramming  effort. 

The  program  may  require  the  use  of  various  storage  files  depending  on 
the  options  and  procedures  selected.  A  number  of  files  are  cjirectly  referenced 
by  the  program,  while  in  the  case  of  the  RESTART  option  the  file  code  is  read  in 
with  the  input  data.  In  the  latter  case  the  file  may  be  any  valid  numeric  file 
code  that  doesn't  interfer  or  conflict  with  other  files  used  in  the  run.  The 
following  is  a  list  of  the  files  which  may  receive  direct  reference  by  the 
program.  In  all  cases  except  the  system  I/O  files  they  are  assumed  to  be 
sequential  binary  files. 


FILE  CODE 
01 
02 
03 
04 

05 

06 

08 

09 

10 

11 

12 

20 


FUNCTION 

Save  RESTART  data  for  DEAD  analysis 

Save  RESTART  data  for  LIVE  analysis 

Save  RESTART  data  for  DYN  analysis 

Element  tension  dynamic  response  file  from 
FREQ  to  be  used  in  CHEK  option 

System  input  file 

System  output  file 

Ship  Motion  File 

Scratch  storage  for  MODE  and  EREQ  options  to 
save  the  data  from  the  previoi^  SAO.  In 
addition,  the  FREQ  option  writes  configuration 
data  for  the  drift  force  iterations. 

Ship's  Static  Load  File 

Scratch  file  used  in  FREQ  to  store  frequency 
independent  portions  of  motion  equations. 

Scratch  file  used  in  FREQ  to  store  the  response 
amplitude  operators  for  the  frequency  span. 

Optional  file  generated  by  MODE  containing 
natural  frequencies  and  mode  shapes. 


21 


Scratch  file  for  MODE 


The  files  other  than  05  and  06  are  referenced  only  when  the  options 
selected  imply  or  directly  request  the  use  of  a  file.  Those  designated  as 
scratch  files  (09,  11,  12)  are  used  anytime  the  indicated  options  are  called 
and  they  may  be  overwritten  in  the  progress  of  the  solution.  The  Ship  Motion 
File  (unit  08)  is  used  only  as  an  input  file  which  is  reference  only  for  FREQ 
when  there  is  a  ship  in  the  system.  The  Ship's  Static  Load  File  (unit  10)  may 
be  generated  with  input  from  SEADYN/DSSM  or  used  as  an  input  file.  The  remaining 
files  are  referenced  only  upon  direct  request  through  an  input  parameter.  Unit  01 
is  written  only  when  INP( 11)^0  on  a  DEAD  SAO.  Similarly  units  02  and  03  are 
written  only  when  INP(11)?‘0  on  LIVE  and  DYN  SAO's,  respectively.  Unit  04  is 
written  only  when  INP(6)^0  on  a  FREQ  SAO.  It  is  expected  to  be  read  when  CHEK 
immediately  follows  FREQ  and  INP(1)^0  on  the  CHEK  SAO  card.  Unit  20  is  written 
only  when  INP(2)<0  on  a  MODE  SAO  card. 

3.5  THE  INITIAL  CONFIGURATION  PROBLEM 

The  use  of  stiffness  methods  to  obtain  the  initial  configuration  of 
cable  systems  presents  some  special  problems.  Generally  the  initial  lengths  and 
the  connectivity  of  the  cables  are  known  but  a  set  of  nodal  coordinates  and 
consistent  element  tensions  are  not.  In  order  to  have  a  nonsingular  stiffness 
matrix  (which  is  required  in  all  of  the  SEADYN  static  solutions)  "reasonable" 
estimates  of  the  configuration  and  element  loads  are  required.  One  approach  is 
to  arrange  the  unstretched  elements  in  a  compatible  configuration  and  assign  an 
arbitrary  tension  to  each  of  the  elements  or  specify  numerical  damping.  The 
tension  or  numerical  damping  is  required  since  most  cable  systems  are  unstable 
without  preloads.  Whenever  only  two  unloaded  cable  elements  are  connected  to  a 
node  in  three  dimensions,  there  is  a  singularity.  Whether  or  not  there  is  a 
singularity  with  more  than  two  elements  at  a  node  depends  on  the  spatial  arrange¬ 
ment.  By  specifying  some  preload  in  the  elements  the  stiffness  matrix  is  made 
nonsingular  because  of  the  geometric  stiffness  term.  The  introduction  of  numerical 
damping  (Section  3. 2. 3. 2)  also  eliminates  the  singularity  by  adding  a  non-zero 
quantity  to  each  diagonal  element  of  the  stiffness  matrix.  Some  caution  is 
necessary  when  using  either  an  assumed  tension  or  numerical  damping  with  the  SLI 
method  (Solution  Option  1)  since  both  of  them  alter  the  equations  of  equilibrium 
and  may  result  in  significant  solution  errors.  The  RFB  incremental  solution 


(Solution  Option  2)  Is  also  affected  but  to  a  lesser  extent  because  of  the  self- 
correcting  nature  of  the  method.  The  MNR  solution  (Solution  Option  0)  will  not 
experience  any  solution  errors  from  this  procedure  but  the  rate  of  convergence 
and  solution  stability  will  be  affected. 

Often  it  is  not  possible  to  develop  a  compatible  arrangement  of  elements 
and  maintain  their  unstretched  lengths.  This  generally  occurs  in  prestressed 
networks.  In  such  cases  a  rational  approximation  of  an  equilibrium  configuration 
can  be  used  in  conjunction  with  the  constitutive  relations  and  a  list  of  the 
unstretched  lengths  to  estimate  the  preloads  and  generate  the  stiffness  matrix. 

If  this  approximate  configuration  results  in  a  singularity  the  previous  device 
of  introducing  some  artificial  tension  in  critical  elements  or  using  numerical 
damping  could  be  used.  In  this  approach  the  MNR  procedure  is  the  only  one  that 
is  appropriate. 

Each  of  the  above  procedures  can  be  initiated  by  appropriate  input  on 
the  cable  element  cards  (Section  3. 2. 2. 7).  A  start  with  an  unstretched  compatible 
configuration  (represented  by  the  input  nodal  positions)  using  numerical  damping 
to  avoid  a  singularity  is  obtained  by  setting  ISIGO  »  0  on  the  first  element 
card,  leavinn  all  of  the  element  tensions  blank,  and  specifying  the  desired 
numerical  damping  coefficient  on  the  Solution  Parameter  Card  (Section  3. 2. 3. 2). 

If  the  tension  specification  is  used,  ISIGO  =^1  and  the  tensions  areinput  or  estimated 
by  the  product  of  the  EAFAC  and  the  initial  £A  from  the  element  constitutive 
relations. 

The  third  approach  is  signalled  by  setting  ISIGO  =  -2  on  the  first  element 
card  and  inputting  the  unstretched  lengths  on  each  element  card.  Tensions  will 
then  be  estimated  from  the  lengths  implied  in  the  nodal  input.  When  an  element 
tension  estimate  is  zero  the  EAFAC  is  used  as  described  above. 


An  additional  specification  for  ISIGO  is  +^3.  This  presumes  no  estimate 
of  the  configuration  has  been  made  and  the  only  correct  nodal  input  is  the 
locations  of  the  fixed  and  moved  nodes.  The  other  nodal  positions  are  arbitrary. 

The  element  cards  then  must  give  the  estimated  tensions  (EAFAC  used  with  ISIGO  =  -3) 
and  the  unstretched  lengths.  Only  the  MNR  solution  should  be  used  with  this 
option.  The  performance  of  this  technique  has  generally  been  unsatisfactory. 
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It  should  be  noted  that  numerical  damping  can  be  specified  In  addition 
to  any  of  these  procedures. 

Another  device  which  can  be  employed  with  incremental  analyses  Is  made 
available  through  the  SAO  Card  (Section  3.2.3.1).  The  Input  parameter  INP(5) 
allows  the  first  load  step  to  be  subdivided  Into  a  given  number  of  non-uniform 
steps.  With  INP(5)  =  n  the  size  of  each  sub-step  is  given  by 
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n(n+l)  (2n+l) 
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where  P^  is  the  size  of  the  first  load  step.  The  procedure  represents  an  attempt 
to  minimize  the  errors  induced  by  the  first  step  and  the  starting  estimate  and 
It  can  be  used  with  solution  Options  1  and  2.  Solution  Option  10  uses  this  same 
procedure  with  the  RFB  method  until  the  sub-steps  are  completed  and  then  shifts 
to  the  MNR  solution. 


Initial  configuration  determinations  with  the  MNR  method  on  soft  systems 
quite  often  leads  to  divergence  or  very  slow  convergence  requiring  very  small 
steps.  In  such  cases  it  has  been  found  more  cost  effective  to  use  either  of  the 
first  two  starting  procedures  (ISIGO  ■  0  or  1)  in  conjunction  with  the  RFB 
solution.  Once  all  of  the  load  steps  are  completed  the  final  configuration  is 
then  used  as  the  initial  configuration  for  a  single  step  MNR  correction  to  the 
result.  This  then  means  that  the  first  two  subanalyses  would  be  DEAD  LOAD 
solutions;  the  first  using  Solution  Option  2  with  from  10  to  100  steps  and  the 
second  with  Solution  Option  0  with  1  step. 


3.6  DEFAULT  DRAG  COEFFICIENTS 

In  the  event  that  fluid  loading  is  required  and  the  drag  coefficients 
are  not  given  by  the  subroutine  DRAGCO  the  program  will  select  the  drag  coefficients 
from  a  set  of  built-in  functions  referred  to  as  default  coefficients.  The  default 
coefficients  are  used  whenever  INDRAG  =  0  (Parameter  Card,  Section  3. 2. 2. 2), 
regardless  of  the  Drag  Coefficient  Codes  on  the  cable  materials  and/or  the  lumped 
bodies.  If  INDRAG  /  0  the  default  coefficients  will  be  used  for  those  components 
which  have  zero  Drag  Coefficient  Codes. 
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3.7  COMMENTS  ON  THE  DYNAMIC  ANALYSIS  OPTIONS 

The  various  possibilities  for  dynamic  time  domain  analyses  range  from 
traditional  small  displacement  linear  analyses  to  fully  nonlinear  iterative 
solutions  of  the  total  equations.  There  are  three  input  parameters  which  control 
the  specification  of  the  solution  procedure.  These  are  the  Solution  Option 
(INP(3)),  the  Update  Option  (INP(4))  and  the  Update  Time  (DTU).  These  parameters 
control  the  selection  of  the  solution  method  (equations  used)  and  the  update  of 
^C  and  ^C.  The  update  of  ^C  occurs  at  DTU  intervals  and  *^C  is  updated  at  INP{4)*DTU 
intervals.  Since  the  tangent  stiffness  matrix  is  recalculated  only  when  ^C  is 
updated  it  is  possible  to  do  small  displacement  analyses  with  INP{3)  =  1  or  2, 

INP(4)  =  0  and  DTU  =  TMAX.  An  incremental  solution  which  takes  a  few  linear 
steps  before  recalculating  the  stiffness  matrix  would  result  from  INP(3)  =  1  or 
2,  INP(4)  =  n  ^  0,  and  DT  <  DTU  <  TMAX.  Fully  nonlinear  solutions  are  obtained 
with  solution  options  0  and  5.  Both  represent  iterative  solutions  to  the  total 
equations.  Option  0  (the  MNR  method)  requires  the  generation  of  the  tangent 
stiffness  and  the  selection  of  INP(4)  and  DTU  will  alter  the  procedures  somewhat 
in  regard  to  the  computational  effort  and  some  differences  in  cost  and  round-off 
errors  may  be  seen  with  variations  in  these  parameters.  Option  5  does  not  use 

p 

the  stiffness  matrix  and  thus  there  is  no  advantage  to  be  had  in  updating  C 
and  ^C.  If  the  updates  are  used  with  option  5  one  might  expect  some  slight 
differences  in  round-off  errors  because  of  the  way  the  displacements  are  summed. 

The  iterative  solutions  allow  two  types  of  predictors.  An  explicit 
predictor  (Newmark  Beta  =  0)  is  used  for  Solution  Option  5.  An  RFB  predictor 
will  be  used  with  the  DI  method  with  Solution  Option  =  -5.  In  this  case  the 
moving  boundary  calculations  may  not  be  correct  since  the  RFB  predictor  does 
not  deal  with  moving  boundaries  consistently.  The  MNR  solution  uses  the  RFB 
predictor  for  Solution  Option  =  0  and  the  explicit  predictor  for  Solution  Option 
=  10. 


4.0 


OUTPUT  DESCRIPTION 


The  SEADYN/DSSM  Program  provides  various  types  and  amounts  of  output 
depending  on  the  option  selected  and  the  value  given  for  various  parameters. 

Three  general  levels  of  output  are  provided.  The  first  is  a  reiteration  of  the 
input  data  augmented  with  supporting  data  to  allow  the  user  to  check  his 
specification  of  the  problem.  The  second  is  the  results  of  the  requested  computa¬ 
tions.  These  two  constitute  the  normal  output  for  the  program.  The  next  section 
describes  these  two  levels  of  output  in  some  detail.  The  third  level  consists 
of  various  amounts  of  intermediate  calculations  which  can  be  used  as  an  aid  in 
input  correction,  solution  optimization  and/or  program  debug.  This  output  is 
described  briefly  in  Section  4.2. 

4.1  NORMAL  OUTPUT 

The  problem  input  data  is  output  with  labels  and  additional  data 
generated  by  input  request.  The  final  results  of  node  and  element  generations 
appear  in  the  output  rather  than  the  Individual  cards  that  caused  the  generation. 
When  the  system  involves  ships  and  data  for  the  ship's  loading  is  obtained  from 
a  predefined  Ship  Load  File,  the  title  and  unit  Tables  are  printed  to  indicate 
which  set  of  load  data  is  being  used.  In  a  few  instances  the  input  data  is 
augmented  with  additional  data  to  aid  in  checking  the  input  (e.g.,  the  mass  per 
unit  length  of  cable  is  calculated  and  output  with  the  material  property  input 
data). 


Output  for  the  various  subanalyses  varies  with  the  type  of  analysis 
performed.  In  addition  to  listing  the  pertinent  subanalysis  input  data,  the 
DEAD,  LIVE  and  DYN  subanalyses  give  the  spatial  position  of  each  node.  Values 
of  the  nodal  velocity  are  also  given  for  DYN  subanalyses.  The  element  tensions 
are  listed  at  the  right-hand  side  of  the  page.  If  there  are  more  elements  than 
nodes  the  excess  element  tensions  are  listed  across  the  page  under  the  nodal 
output.  This  set  of  data  is  repeated  for  each  requested  load  or  time  step.  The 
input  parameters  INP{1)  and  INP(2)  (Section  3.2. 3.1)  control  the  frequency  of 
output  for  DEAD  and  LIVE  subanalyses.  INP(l)  gives  the  number  of  steps  used  to 


apply  the  load  and  INP(2)  indicates  how  many  steps  are  taken  before  output  is 
printed.  The  time  step,  DT,  replaces  INP(l)  in  DYN  subanalyses.  In  all  cases 
the  state  at  the  beginning  of  the  subanalysis  is  printed.  The  results  of  every 
step  are  printed  regardless  of  changes  in  step  size  when  INP(2)  »  1.  With 
INP(2)  greater  than  one,  the  output  interval  is  held  approximately  the  same 
even  though  the  step  size  may  change.  For  example,  if  the  original  time  step 
is  0.01  seconds  and  INP(2)  =  5,  then  output  will  occur  at  intervals  of  0.05 
seconds  even  though  the  step  size  may  be  altered  by  the  solution  procedure.  If 
the  step  size  is  changed  such  that  it  will  no  longer  generate  data  at  the  0.05 
second  interval  the  output  will  occur  as  nearly  as  possible  at  the  0.05  second 
intervals.  The  results  for  any  of  the  output  records  can  be  saved  on  a  file 
for  subsequent  restart  usage.  This  is  controlled  by  INP(ll).  When  INP{11)  =  i  n 
then  contents  of  the  COMMON  storage  areas  are  output  on  the  appropriate  file  for 
every  nth  output  record.  (See  Section  3.2.4  for  more  details.)  When  an  output 
record  is  saved,  a  message  is  printed  to  indicate  the  record  number  on  the  file. 

Output  for  the  MODE  subanalysis  consists  of  the  natural  frequency 
and  period  and  the  mode  shape  for  each  of  the  natural  modes  represented  in  the 
model.  The  solution  procedure  does  not  solve  for  the  modes  in  any  specific 
order  but  output  may  be  requested  in  the  order  of  increasing  natural  frequency. 

An  option  is  provided  to  save  the  frequency  and  mode  shape  data  on  an  output 
file  for  use  external  to  the  SEADYN/DSSM  program.  The  file  format  is  described 
in  the  Program  Maintenance  Manual  [5]. 

The  output  from  the  FREQ  subanalysis  consists  of  reiteration  of  the 
subanalysis  input  data,  identification  data  from  the  ship  motion  file,  regular 
wave  solution  data,  random  response  results  and  regular  wave  results.  The 
identification  data  printed  from  the  ship  motion  file  are  the  title,  unit  labels, 
ship's  parameters,  nondimensional  wavelengths,  headings  and  roll  angles  (in 
radians)  as  they  were  obtained  from  the  file.  The  output  data  for  the  regular 
wave  solutions  gives  the  wavelength,  amplitude  and  slope  for  each  wave.  If  the 
amplitude  is  below  the  specified  cutoff  value,  a  message  is  printed  and  the 
calculation  proceeds  to  the  next  frequency.  If  the  frequency  is  above  the  given 
upper  bound  and  the  wave  amplitude  is  below  the  cutoff,  the  calculation  is 
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assumed  to  be  completed.  When  a  ship  is  In  the  system,  the  ship  response  at 
each  frequency  is  listed  in  the  local  coordinate  system  (surge,  sway,  heave, 
roll,  pitch,  yaw)  normalized  to  a  unit  wave  amplitude.  This  represents  the 
constrained  ship's  response  after  the  iterations  on  roll  damping  are  completed. 

The  data  given  are  the  amplitude  and  phase  angles  (degrees)  for  each  component. 

The  steady  state  drift  force  components  for  the  stated  wave  amplitude  are  also 
given.  If  INP(9)  is  not  zero  the  ship  response  data  are  given  for  the  unrestrained 
ship  at  each  stage  of  the  roll  damping  iterations.  This  output  is  preceded  by 
the  complex  components  of  the  unrestrained  ship  response  (again  in  the  local 
coordinate  system). 

When  all  frequency  data  has  been  processed,  the  total  combined  drift 
forces  from  all  of  the  regular  waves  are  printed.  The  values  given  are  for 
surge,  sway  and  yaw.  If  updates  in  the  static  reference  state  are  to  be  calculated, 
the  output  for  a  LIVE  subdnalysis  will  follow.  Depending  on  the  option  selected, 
the  regular  wave  solutions  will  be  repeated  until  the  static  reference  state 
converges  or  calculation  will  proceed  to  the  next  stage. 

At  the  completion  of  the  frequency  scan  and  any  reference  configuration 
updates,  the  user  has  the  option  of  requesting  random  response  calculations  and/or 
detailed  data  for  the  regular  wave  response  of  any  node.  Output  for  the  random 
responses  consists  of  the  points  in  the  input  and  response  spectra  listed  for 
each  frequency  and  the  statistical  estimates  for  the  average  of  the  1/3,  1/10 
and  1/100  highest  responses.  All  six  components  of  ship  responses  are  listed  in 
local  coordinates.  Tension  responses  give  the  reference  static  value  and  the 
combined  static  and  dynamic  values.  If  the  ultimate  strength  is  provided  on 
the  material  data  card  (Section  3. 2. 2. 8)  and  load  factors  are  given  on  the 
random  response  request  (Section  3.2.3.10)  then  the  tension  output  will  also 
include  a  listing  of  the  ultimate  strength  and  the  sum  of  the  factored  static 
and  dynamic  components.  Random  response  data  for  nodes  consists  of  the  statistical 
values  (1/3,  1/10,  1/100)  for  the  individual  dynamic  displacement  component  in 
the  global  coordinate  direction.  The  position  of  the  node  in  the  static  reference 
state  is  also  printed. 


Regular  wave  response  output  for  specified  nodes  consists  of  the  wave 
frequency  and  amplitude  for  which  the  results  are  given,  the  amplitude  and  phase 
angles  (degrees)  of  each  of  the  nodal  components  and  the  values  of  the  response 
through  one  full  cycle.  The  output  may  be  either  the  dynamic  displacements  In 
the  global  or  ship's  local  coordinates  or  the  nodal  position  In  global  coordinates. 

At  the  completion  of  the  calculations  for  an  individual  wave  heading 
the  calculations  may  be  repeated  for  other  wave  headings  or  the  FREQ  subanalysis 
may  be  terminated.  If  a  request  has  been  made  to  save  random  dynamic  response 
data  for  subsequent  adequacy  checks,  a  message  will  be  printed  which  Identifies 
the  record  number  on  the  file  which  saves  the  results. 

The  output  for  the  CHEK  subanalysis  has  three  variations  on  one  general 
form.  The  general  form  is  a  list  of  the  input  parameters,  a  printout  of  the 
Imposed  loading,  and  messages  to  signal  when  capacities  are  exceeded  or  unusual 
situations  are  encountered.  If  component  inventories  are  used  to  get  the  capacity 
of  the  Item  being  checked,  the  inventory  parameter  which  is  closest  to  the  value 
given  in  the  Input  will  be  printed  followed  by  the  calculated  capacity.  This 
calculation  takes  the  capacity  In  the  Inventory,  converts  It  to  consistent  units 
and  divides  by  the  safety  factor. 

The  capacity  checks  are  made  on  the  state  defined  when  CHEK  is  called. 

If  the  preceding  subanalysis  was  DEAD,  LIVE  or  DYN  it  deals  with  the  state  at  the 
end  of  the  subanalysis.  If  the  preceding  subanalysis  was  FREQ  and  a  nonzero 
record  number  is  given  (INP{1)),  the  static  reference  state  will  be  evaluated 
first  and  then  the  increased  values  for  the  spectral  responses  will  be  printed. 

Output  for  checks  of  anchor  adequacy  give  the  total  load  and  vertical 
component  applied  to  the  anchor.  The  angle  in  degrees  between  the  total  load 
resultant  and  the  horizontal  is  also  given.  A  message  is  printed  If  the  angle 
implies  the  load  resultant  acts  down  on  the  anchor.  This  is  an  indication  that 
the  line(s)  actually  touch  bottom  before  reaching  the  anchor  and  the  model  was 
not  detailed  enough  to  sense  this.  Recall  the  discussion  in  Section  2.2.5.  If 
the  angle  to  the  resultant  load  exceeds  3°  and  the  inventory  has  been  used,  a 
message  will  be  printed.  If  spectral  response  data  has  been  provided,  the 


output  for  the  reference  state  is  followed  by  a  listing  of  the  resultant  of 
the  dynamic  components  which  represent  the  average  of  the  highest  1/3,  1/10,  1/100 
responses.  These  are  followed  by  similar  values  representing  the  sum  of  the 
static  and  dynamic  parts.  These  are  followed  by  similar  lists  for  the  vertical 
load  components. 

Output  for  adequacy  checks  on  buoys  is  similar  to  that  for  anchors. 

The  total  resultant  load,  its  vertical  component  and  angle  between  the  resultant 
and  the  horizontal  are  printed.  The  resultant  load  should  be  vertical  and  the 
angle  should  be  90**.  This  is  because  the  buoy  is  assumed  to  have  no  external 
constraints  in  any  direction  other  than  vertical  and  the  attached  lines  must 
therefore  balance  the  horizontal  loading.  If  the  component  Inventories  are  used, 
the  buoy  O.D.  from  the  Inventory  which  is  closest  to  the  input  value  will  be 
printed.  This  will  be  followed  by  the  calculated  capacity  which  is  obtained  by 
dividing  the  upper  limit  inventory  value  by  the  safety  factor.  If  the  vertical 
component  of  the  resultant  exceeds  the  buoy  capacity  a  message  will  be  printed. 
Spectral  response  data  follows  the  same  pattern  as  in  the  anchor  case. 

Line  capacity  checks  give  the  nominal  line  load  (the  element  tension 
determined  in  the  analysis)  and  the  element  end  loads  which  are  obtained  from 
the  procedures  outlined  in  Section  2.2.5.  The  global  components  of  the  end  loads 
are  also  given.  If  the  component  inventory  is  used  the  line  diameter  from  the 
inventory  which  is  closest  to  the  input  value  will  be  printed  followed  by  the 
calculated  capacity  (the  inventory  value  divided  by  the  safety  factor).  A 
message  will  be  printed  if  the  line  capacity  is  exceeded  by  the  largest  of  the 
end  loads.  Spectral  response  data  given  for  lines  are  the  dynamic  tensions 
(1/3,  1/10,  1/100)  and  the  combined  tensions  obtained  by  adding  the  largest  end 
load  to  the  dynamic  tension  estimates. 

The  CHEK  subanalysis  also  provides  for  the  printing  of  the  component 
inventories.  The  output  will  follow  the  form  given  in  Appendix  C. 


4.2  OPTIONAL  DEBUG  OUTPUT 

The  integer  field  in  columns  51-55  on  the  Master  Control  Card 
(Section  3. 2. 2.1)  and  on  the  Subanalysis  Option  Card  (Section  3.2.3. 1)  is  used 
to  signal  requests  for  extra  output.  In  general,  this  output  was  intended  for 
use  in  debugging  and  checkout  of  the  program.  It  has  also  been  found  useful  in 
finding  difficult  input  errors  and  analyzing  the  efficiency  of  various  choices 
of  solution  methods  and  parameters.  The  output  is  not  usually  requested  and 
little  effort  has  been  expended  in  generating  labels  to  identify  the  output.  It 
is  assumed  that  when  the  higher  levels  of  output  detail  are  requested  that  the 
user  is  familiar  enough  with  the  program  structure  and  the  solution  methods  to  be 
able  to  utilize  the  data  intelligently.  Indiscriminant  requests  for  extra  output 
is  strongly  discouraged  since  it  greatly  adds  to  the  output  volume  and  increases 
run  costs. 

The  Fortran  variable  IBG  is  the  item  read  to  control  the  extra  output. 
Generally,  low  levels  of  output  are  obtained  by^small  values  for  IBG.  Thus  when 
IBG  is  less  than  7  only  a  small  amount  of  specific  output  is  produced.  Values 
of  IBG  greater  than  6  will  give  significant  amounts  of  output  on  iterative 
solutions.  The  following  is  a  list  of  the  values  of  IBG  recognized  in  the 
subanalyses  and  a  brief  description  of  the  output  they  generate.  (Frequency 
codes:  I  =  each  iteration,  S  =  each  step,  F  =  each  wave  frequency,  H  •=  each 
wave  heading) 

IBG  FREQUENCY  EXTRA  OUTPUT 

2.  1  S  The  number  of  iterations  required  to  get  convergence 

of  iterative  solutions  (DEAD,  LIVE,  DYN). 

p 

A  message  whenever  the  reference  configuration,  C, 
is  updated. 

2  I  A  message  whenever  an  element  is  in  compression 

with  no  compression  stiffness  (i.e.,  slack). 

A  message  whenever  a  buoy  or  anchor  changes  its 
status  relative  to  the  surface  and  bottom  limits. 

The  values  of  the  net  vertical  forces  applied  to 
buoys  and  anchors. 


EXTRA  OUTPUT 


The  values  of  nodal  displacement,  velocity  and 
acceleration  In  the  nodal  displacement  and 


the  nodal  displacement  and 


acceleration  at  the  Incremental  reference  state. 


the  nodal  positions  and  forces  In 
node. 


for  each 


Payout  update  data  (not  Implemented  yet) 

Ship's  static  loads  In  local  and  global  coordinates. 

For  WR  Solution:  five  values  are  printed  — 
residual  norm,  displacement  norm  (the  one  used), 
the  values  for  the  three  alternative  displacement 
norms. 

For  the  PI  Dynamic  Solution:  six  values  are 
printed  —  the  Iteration  number,  the  displacement 
norm,  the  maximum  acceleration  Increment  and  the 
values  for  three  alternative  displacement  norms. 

For  the  FREQ  Subanalysis;  the  values  of  the  steady 
state  solution  vector,  the  response  amplitude 
operator  data  for  displacements  and  tensions. 

For  the  MNR  Solution;  A  message  when  the  tangential 
stiffness  matrix  Is  calculated.  A  message  when 
the  1-D  search  Is  Initiated  and  terminated. 

Element  tensions. 

The  values  of  DMU,  BETA  and  the  mod1fV:d  diagonal 
terms  of  the  tangential  stiffness  matrix  when 
numerical  damping  Is  used. 

A  message  when  alternating  estimates  are  detected 
and  the  next  "guess"  is  between  them. 

The  current  positions  of  the  nodes. 

For  the  PI  Dynamic  Solutions:  A  message  whenever 
alternating  estimates  are  detected  and  the  next 
"guess"  is  between  them. 

For  the  RFB  Dynamic  Solutions:  A  message  when 
tangent  stiffness  matrix  calcul a  ted . 

Element  tensions 

Effective  loads 

Incremental  displacements 


For  the  FREQ  Subanalysis;  The  values  of  the  ship's 
roll  angle. 


For  the  MNR  Solutions:  Current  values  of  the 
residual  norm,  the  components  of  the  residual, 
the  force  vector,  the  total  displacement,  and 
the  incremental  displacement. 

IBG  is  reset  to  7  after  5  steps. 

For  the  PI  Dynamic  Solutions:  Current  values  of 
the  force  residual,  acceleration  increment,  and 
displacement  increment. 

For  the  FREQ  Subanal^ysis;  The  values  for  the 
wave  induced  ships  forces  (local  coordinates), 
the  force  vector  (global  coordinates),  steady 
state  response  vector. 

Details  of  dynamic  tension  calculation 

For  the  RFB  and  SLI  Static  Solutions:  The 


incremental  or  effective  forces  and  the  displacement 
increments. 

For  RFB  Dynamic  and  MNR  Solutions;  The  terms 
of  the  tangential  stiffness  matrix  in  compact 
storage  format. 

IBG  set  to  8  after  one  printing. 

For  FREQ  Subanalysis;  The  dimensional  1  zed  form 
of  the  ships  mass  and  restoring  matrices  (M^  and  K^) 

2 

Diagonal  terms  of  the  system  matrix,  [K-u  M  +  iwC]. 

The  parameters  used  in  interpolating  in  the 
Ship  Motion  File. 

The  ship's  added  mass  and  damping  matrices  (M.-,  C-) 
and  loading  data.  ^ 

The  values  of  fluid  velocity  and  fluid  induced 
loads  at  each  node  in  the  system. 

The  local  and  global  forms  of  the  ship's  static 
loads  (same  as  IBG  «  5) 

The  local  to  global  transformation  matrix  for  ships. 


>11 


H 


The  mass  matrix  In  compact  form  for  the  mooring 
system  in  a  FREQ  subanalysis. 

Once  The  global  stiffness  matrix  in  compact  form. 

IBG  set  to  11  after  one  printing. 

Once  Each  element  stiffness  matrix  and  local  to  global 

transformation  data. 

The  element  contributions  to  the  residual.  IBG  set 
to  11  after  one  printing. 

Once  Same  as  IBG  -  13  except  IBG  is  set  to  zero  after 

one  printing. 

A  nonzero  value  for  IBG  on  the  Master  Control  Card  will  cause  the 
printing  of  the  components  of  the  lumped  mass  matrix,  the  residual  added  mass 
(the  added  mass  to  be  removed  from  the  tangential  direction),  the  gravity  load 
vector,  plus  a  list  of  the  direction  cosines,  the  initial  length,  and  the 
unstretched  length  for  each  element  in  the  input  initial  configuration.  This 
latter  output  has  been  found  to  be  very  helpful  in  identifying  erroneous  node 
and  element  input.  If  IBG  *  3  the  same  data  as  in  the  subanalysis  list  for 
IBG  =  3  will  be  printed.  This  represents  the  initial  reference  state  and  all 
displacements,  velocities,  accelerations,  and  forces  should  be  zero.  The  nodal 
positions  should  be  the  same  as  those  printed  in  the  normal  output. 


>12 

>13 

14 


5.0 


DATA  CHECKS  AND  ERROR  RETURNS 


The  SEADYN/DSSM  program  makes  some  checks  of  the  Input  'ata  and 
attempts  to  aid  the  user  find  his  errors  by  printing  various  messages.  No 
attempt  has  been  made  to  be  comprehensive  In  this  feature  since  It  Is  very 
difficult  to  foresee  and/or  detect  many  of  the  possible  errors.  During  the 
processing  of  the  subanalyses,  checks  are  made  on  the  validity  of  the  requests 
and  the  convergence  of  the  solution  procedures.  The  messages  that  may  be 
printed  are  listed  below  with  a  brief  description  of  the  probable  cause  and/or 
cure.  The  action  taken  after  the  error  detection  Is  indicated  by  the  following 
codes : 

(F)  Fatal,  run  aborted 

(N)  Abort  analysis  case  and  seek  a  new  problem  definition 
by  searching  the  deck  for  a  NEW  SAD  card. 

(D)  Abort  present  SAD  activity  and  go  to  the  next  SAD  card. 

(S)  Skip  this  request  and  go  to  the  next  card. 

(C)  Continue  calculation  with  action  as  indicated. 


5.1  SYSTEM  DESCRIPTION  CHECKS 

(N)  NODES  OMIHED  IN  GENERATION 

Message  is  followed  by  a  list  of  ones  and  zeros  (ten  per 
line)  corresponding  to  the  nodes  In  the  system.  The 
zeros  indicate  which  nodes  were  not  defined  either  by 
input  or  generation  sequence.  All  nodes  must  be  accounted 
for.  Following  the  integers  the  nodal  coordinates  are 
printed  to  aid  In  finding  the  error.  The  problem  usually 
comes  from  improper  node  generation  input. 


(N)  CATENARY  GENERATION  CASE  2  NOT  WORKING 

Message  is  followed  by  five  numbers.  They  are  values 
for  the  following  items: 

X  =  the  distance  from  node  N1  to  the  point 
where  the  catenary  is  tangent  to  the 
bottom.  In  this  case  it  will  be 
negative  indicating  the  first  node  is 
above  the  tangent  point. 

X.p  =  horizontal  distance  from  node  N2  to  the 
tangent  point 


XSPAN  =  Horizontal  distance  between  N1  and  N2 


YSPAN  =  Vertical  distance  between  N1  and  N2 
C  =  H/w  for  catenary 


The  option  for  generating  a  non-tangent  catenary  (Case  2) 
has  not  been  developed.  If  this  was  not  intended,  check 
the  input  for  the  positions  of  nodes  N1  and  N2  and  the 
values  for  H  and  w. 


ELEMENT  GENERAlTON  ERROR 

First  element  was  not  input,  or  element  cards  are  not 
in  increasing  element  number  order,  or  last  element 
was  not  input. 


MORE  THAN  TEN  ELEMENTS  CONNECTED  TO  BUOY/ANCHOR  xxx 

The  number  given  indicates  the  buoy/anchor  number.  Array 
dimensions  presently  limited  to  ten.  Check  element  cards 
or  reformulate  description. 


MOORING  BUOY  AT  NODE  xxx  HAS  TOO  MANY  SLAVES. 
FOUND  xxx 

Present  dimension  limits  it  to  ten. 


IMPROPERLY  DEFINED  MOORING  BUOY  AT  NODE  xxx 
NO.  OF  SLAVES  FOUND  =  xxx 

Mooring  buoys  require  at  least  two  slaves. 


SHIP  DATA  INPUT  ERROR  ON  xxx  (Ship  No.) 

NO.  OF  SHIPS  ON  FILE  =  xxx 
LOAD  FUNCTION  OPTION  =  xxx 

Blank  Card  for  moored  ship  data  requested  definition 
of  ship  from  load  file  with  no  load  file  defined. 

OR 

Attempted  ship  scaling  with  no  load  file  defined. 
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(F)  ERROR  IN 

HEADING  =  xxxx 

LAST  TABLE  ENTRY  =  xxxx 

SYMMETRY  FLAG  =  xxxx 

Heading  requested  which  exceeds  the  largest  value  in 
the  table.  Check  ship  load  input  table. 


(F)  DEPTH  CORRECTION  ERROR  —  SHIPS  DRAFT  EXCEEDS  WATER  DEPTH 
Gross  input  error.  Check  for  dropped  deck. 


(C)  WARNING— UNITS  DO  NOT  APPEAR  TO  BE  CONSISTENT 
GRAVITATIONAL  ACCELERATION  FROM  TAPE  IS  xxxx 
UNIT  LABELS  FROM  TAPE  ARE  xxx  xxx  xxx 

The  ship  motions  file  conversion  factors  do  not 
properly  convert  the  GRAV  on  the  file  to  the 
GRAV  specified  in  this  run.  Calculation  still 
proceeds . 


(N)  TAPE  POSITIONING  OR  FORMAT  ERROR 
UNABLE  TO  FIND  SHIP  DATA 
ITEM  xxx  LAST  READ  IS  FOR  HEADING  xxxx 
AND  WAVELENGTH  xxxx. 

WANTED  HEADING  xxxx  WITH  WAVE  LENGTH  xxxx. 

The  ship  motion  file  is  not  formulated  properly  or 
other  input  or  equipment  malfunction  has  made 
reading  the  file  impossible. 


5.2  SUBANALYSIS  OPTION  ERRORS 

(N)  UNRECOGNIZED  ANALYSIS  OPTION  =  xxxx 

TRY  TO  GET  TO  NEXT  CASE 

Usually  indicates  improper  numbers  of  cards  or 
cards  out  of  sequence. 


INSUFFICIENT  STORAGE  TO  PROCEED 
COMMON  SIZE  =  XXX 
STORAGE  NEEDED  =  xxx 
HALF  BANDWIDTH  =  xxx 
DEGREES  OF  FREEDOM  =  xxx 
BASE  SIZE  =  xxx 

Subanalysis  request  (DEAD, LIVE, DYN)  cannot  be 
processed  due  to  storage  limitations.  Problem 
must  be  reformulated  or  NCOM  increased  (see 
Section  3.4) 


NOT  ENOUGH  STORAGE  FOR  MODE  SHAPES 
NEED  xxxx,  WITH  xxxx  AVAILABLE 

See  Section  3.4 


NOT  ENOUGH  STORAGE  FOR  FREQUENCY  SOLUTION 
NEED  xxxx,  V'tril  xxxx  AVAILABLE 

See  Section  3.4 

INCONSISTENT  ANALYSIS  REQUEST  ON  RESTART 

Attempted  a  restart  of  DEAD, LIVE  or  DYN  and  the 
file  was  not  from  that  type  of  subanalysis. 


TAPE  LABEL  xxxxxx  DOES  NOT  AGREE  WITH  xxxxxx 

The  label  check  failed  on  restart.  The  first 
six  characters  on  the  RESTART  Title  Card  did 
not  agree  with  the  first  six  characters  on 
the  Analysis  Title  Card  which  saved  the  file. 


BLANK  COMMON  ON  TAPE  LARGER  THAN  SPACE  AVAILABLE 

The  data  saved  on  the  file  requires  more 
storage  than  is  presently  available,  (see 
Section  3.4)  The  two  numbers  printed  are 
the  required  and  current  values  of  NCOM. 


NUMBER  OF  PAYOUT  ENDS  EXCEEDS  LIMIT  OF  5 
Check  problem  specifications. 


NUMBER  OF  SHIPS  =  xxxx 

DYNAMIC  SOLUTION  PRESENTLY  LIMITED  TO  ONE  SHIP 

Frequency  domain  solution  requested  with 
more  than  one  ship. 


SPECTRUM  ERROR.  NO  FREQUENCIES  FOUND  WITH  SIGNIFICANT  WAVE  HEIGHTS 

Check  spectrum  parameters  and/or  frequency 
range. 


INVALID  COMPONENT  TYPE  xxxx 
ASSUME  END  OF  INPUT 

Occurs  when  random  response  requests  are  being 
processed  and  a  card  is  encountered  which  does 
not  have  SHIP,  NODE,  TENS  or  DONE  in  columns 
1-4.  The  items  to  this  point  are  processed 
and  the  case  aborted  with  the  additional  message: 


POSSIBLE  SEQUENCE  ERROR  —  CASE  TERMINATED 

(See  previous  explanation) 

The  dynamic  response  file  is  not  written. 


SHIP  OUTPUT  REQUESTED  WITH  NO  SHIP  IN  THE  SYSTEM 

Random  response  request  for  ship  ignored 
when  NSHIPS  /  1. 


REQUESTED  DYNAMIC  EFFECTS  WITH  NO  FREQUENCY  DOMAIN  FILE  PROVIDED 
ABORT  ADEQUACY  CHECK 

Dynamic  response  file  was  not  saved  in  the 
previous  FREQ  SAO  (see  INP(6)  on  SAO  Card). 


IMPROPER  RECORD  NUMBER  ON  FREQUENCY  DOMAIN  FILE 
ABORT  ADEQUACY  CHECK 

The  requested  record  on  the  dynamic  response 
file  was  either  zero  or  greater  than  the 
number  of  records  written.  Check  INP(l) 
on  CHEK  SAO  Card. 


NO  DYNAMIC  TENSION  PROVIDED  ON  ELEMENT  xxxx 
IGNORE  DYNAMICS 

The  random  response  data  for  this  element  was  not 
requested  by  a  TENS  Card  for  this  wave  heading  in 
the  FREQ  SAO. 


CAPACITY  AND  COMPONENT  ID  BOTH  ZERO 
SKIP  REQUEST  FOR  xxxx  (CTYPE) 

Check  input  card. 


INVALID  COMPONENT  TYPE 
SKIP  REQUEST  FOR  xxxx 

The  component  number  is  not  one  recognized 
by  the  component  inventory.  See  Section 
3.2.3.11  and  Appendix  A. 


THERE  IS  NO  BUOY  AT  NODE  xxx 
SKIP  REQUEST  FOR  xxxx 

Check  node  number  given. 


MORE  THAN  TWENTY  LINES  ON  ANCHOR  AT  NODE  xxx 
SKIP  REQUEST  FOR  xxxx 

The  fixed  node  where  anchor  capacity  check  is 
requested  has  too  many  elements  connected  to 
it  to  make  the  check. 


NO  LINES  CONNECTED  TO  NODE  xxx 
SKIP  REQUEST  FOR  xxxx 

Check  node  number. 


ANCHOR  NOT  ON  BOTTOM  AT  NODE  xxx 
SKIP  REQUEST  FOR  xxxx 

Anchor  weight  is  not  sufficient  to  hold  the 
lines  at  this  node  and  it  has  been  lifted 
off  the  bottom  or  has  not  reached  there. 


TOO  MANY  LOAD  SETS  FOR  DRIFT  FORCE  ITERATION 
REFERENCE  STATE  NOT  UPDATED 


The  previous  static  analysis  used  3  load  sets.  This 
leaves  no  room  to  store  the  drift  forces. 


INVALID  NODE  NUMBER  =  xxx 
SKIP  THIS  REQUEST 

Regular  wave  response  data  was  requested  for  a  node 
number  which  is  less  than  1  or  greater  than  the 
number  of  nodes  in  the  model. 


DATA  NOT  AVAILABLE  FOR  xxxx 
SKIP  THIS  REQUEST 

Regular  wave  response  data  was  requested  for  a  frequency 
outside  of  the  range  which  was  generated  on  the  RAO  file. 


INVALID  CALCULATION  OPTION  =  xxxx 
CASE  TERMINATED 

The  frequency  domain  calculation  option  was  not 
RAND,  REGb  or  DONE.  Cards  are  out  of  sequence 
or  card  mispunched. 


SHIP  MOTION  FILE  ERROR 
WAVELENGTHS  NOT  IN  DECREASING  ORDER 

Check  format  of  Ship  Motion  File. 


SHIP  MOTION  DATA  EXCEEDS  LIMITS 


NOB 

NOH 

NOK 

NRV 

LIMITS 

5 

30 

30 

8 

VALUES  READ 

xxx 

xxx 

xxx 

xxx 

The  ship  motion  file  has  arrays  larger  than  the 
dimension  in  SEADYN/DSSM. 

NOB  =  number  of  Froude  Numbers 
NOH  =  number  of  wave  headings 
NOK  =  number  of  wavelengths 
NRV  =  number  of  roll  angles 


SOLUTION  OPTION  EXECUTION  MESSAGES 


EQUATION  DECOMPOSITION  FAILED  ON  ROW  xxx 

The  simultaneous  equations  In  the  subroutine 
SLVBAN  (DEAD.  LIVE.  DYN.  FREQ)  or  subroutine 
COMBAN  (FREQ)  are  singular  or  sufficiently 
ill-conditioned  to  appear  singular.  The  row 
of  the  matrix  is  calculated  from  (3*(I-1)+J) 
where  I  is  the  node  number  and  J  is  the  direction 
(1,  2  or  3).  Check  node  I  for  proper  constraint. 
This  also  occurs  in  poorly  tensioned  (soft)  initial 
configurations.  Check  for  zero  tensions,  etc. 

(see  Section  3.5) 

If  this  occurs  with  solution  option  1  or  2  or 
in  the  steady-state  response  calculations  in 
FREQ  SAO  the  case  is  terminated. 

If  it  occurs  during  a  MNR  solution  (options  10 
or  0)  various  attempts  are  made  to  remove  the 
singularity  and  repeat  the  step.  Falling  in 
these,  the  case  is  terminated. 


ROW  xxx  OF  STIFFNESS  MATRIX  HAS  NEGATIVE  DIAGONAL  TERM  »  xxxx 

Message  follows  previous  message  on  solution 
options  1  or  2. 


SINGULAR  EQUATIONS  WITH  NUMERICAL  DAMPING  FACTOR  OF  xxxx 

This  message  is  printed  from  the  MNR  solution. 

See  note  1  of  Section  3. 2. 3. 2  for  further 
explanation.  This  message  is  printed  after  the 
three  trials  described  in  that  note.  This  will 
be  followed  by  a  reduction  in  step  size  subject 
to  the  limits  of  INP(14)  on  the  SAO  Card. 


LAST  TWO  RESIDUAL  AND  DISPLACEMENT  NORMS  xxxx  xxxx  xxxx  xxxx 
KOUNT  =  xxxx 

The  four  values  printed  are  the  residual  and 
displacement  norms  for  iteration  i-1  and 
iteration  i  in  the  MNR  solution.  This 
message  signals  an  increasing  norm  indicative 
of  divergence  or  a  lack  of  convergence  in  the 
number  of  iterations  given  by  KOUNT.  This 
will  be  followed  by  a  reduction  of  step  size 
subject  to  the  limits  of  INP(14)  on  the  SAO 
card. 
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(C.N)  STEP  SIZE  REDUCED  TO  xxxx  ON  INCREMENT  xx 
RESIDUAL  AND  DISPLACEMENT  NORMS  xxxx  xxxx 

This  message  follows  the  reduction  of  step 
size  In  a  MNR  solution.  It  will  follow 
either  of  the  previous  two  messages  when 
the  number  step  size  reductions  is  within 
the  limits  imposed  by  INP{14)  on  the  SAO 
card. 


(N)  DIVERGENCE  ON  STEP  xxx  AT  xxxxx  (load  factor  or  time) 

Signals  abort  of  MNR  solution  after  INP(14) 
successive  step  size  reductions. 


(F.N)  BUOY/ANCHOR  NO.  xxx  MOVED  TOO  FAR  BEYOND  LIMIT  IN  ONE  STEP 
DECREASE  STEP  SIZE 

Indicates  a  check  for  surface  or  bottom 
constraints  detected  large  motion  of  the 
body  as  It  approached  the  limit  (e.g.,  buoy 
moved  out  of  water  more  than  buoy  diameter 
In  one  step).  On  the  MNR  Iterations 
(solution  option  0  or  10)  and  DI  solutions 
(option  5)  the  step  size  will  be  decreased. 

In  all  other  situations  the  run  1s  aborted. 


(C)  DIVERGING  NONLINEAR  ITERATION  AT  TIME  xxxx 
STEP  XX  WITH  A  TIME  INCREMENT  OF  xxxx 
LAST  TWO  NORMS  xxxx  xxxx  KOUNT  =  xx 

Signals  a  lack  of  convergence  in  the 
Direct  Numerical  Integration  (option  5) 
time  domain  solution.  The  various 
situations  which  lead  to  this  message 
are: 

a)  Buoy/Anchor  moved  too  far  (see  previous  message) 

b)  Acceleration  more  than  doubled  in  one  iteration 
on  the  component  with  the  largest  acceleration. 

12 

c)  The  displacement  norm  exceeds  1  x  10 

d)  The  displacement  norm  increased  in  three 
successive  iterations 

e)  The  largest  displacement  increment  exceeds 
a  magnitude  of  1  x  lOlO 

f)  Number  of  iterations  exceeds  INP(13)  on  SAO  Card. 
This  will  be  followed  by  a  reduction  in  time  step 
subject  to  the  limits  of  INP(14)  on  the  SAO  card. 


TIME  STEP  REDUCED  TO  xxxx  ON  STEP  xxxx 
NORM  =  xxxx 

This  message  follows  the  reduction  of  step 
size  in  a  DI  solution.  It  will  follow 
the  previous  message  when  the  number  of 
step  size  reductions  is  within  the  limits 
imposed  by  INP(14)  on  the  SAO  Card. 


DYNAMIC  SOLUTION  DOES  NOT  CONVERGE  AT  TIME  xxxx  WITH  A  TIME  INCREMENT  OF  xxxx 
LAST  TWO  NORMS  xxxx  xxxx 

Signals  abort  of  DI  solution  after  INP(14) 
successive  step  size  reductions. 


;• 

6.0  SAMPLE  PROBLEMS 

No  simple  set  of  demonstration  problems  can  explore  all  of  the  capabilities 
of  a  program  such  as  SEADYN/DSSM.  The  program  was  designed  to  be  very  general  In 
the  types  and  complexity  of  problems  It  can  deal  with.  Even  though  only  one 
finite  element  type  Is  presently  available  In  the  program.  It  can  deal  with  a 
wide  variety  of  problems.  Reference  2  explored  some  of  the  capabilities  and 
compared  the  accuracies  and  efficiencies  of  the  various  solution  options  prior 
to  the  Introduction  of  the  moored  ship  capability.  Reference  6  deals  specifically 
with  the  static  and  dynamic  analyses  of  moored  ships,  and  a  number  of  example 
problems  are  presented  there. 

In  this  section  five  sample  problems  will  be  presented.  Four  of  them 
are  from  the  set  discussed  In  Reference  2  and  the  fifth  evaluates  the  static 
response  of  a  moored  ship  to  changes  In  wind  and  current  direction.  The  input 
cards  for  each  case  are  presented  with  the  discussion.  It  will  be  Instructive  for 
the  new  user  to  compare  these  cards  with  the  Input  descriptions  of  Section  3.0. 

6.1  THE  SNAP-THROUGH  RESPONSE  OF  A  TWO-BAR  TRUSS 

A  popular  problem  for  demonstrating  geometrically  nonlinear  effects  [17] 

Is  a  shallow  two-bar  truss  supported  by  a  spring  and  driven  through  the  snap- 
through  range  by  a  single  point  load.  A  comparison  of  the  performance  of  the 
three  static  analysis  solution  options  are  compared  In  Figure  6-1.  The  input 
cards  for  the  3  solutionsare  shown  In  Table  6-1.  The  data  represents  the  SLI 
solution  followed  by  restarts  to  get  the  MNR  and  FBK  solutions. 

6.2  STATIC  RESPONSE  OF  THE  SEACON  II  STRUCTURE 

Kretschmer,  et  al  [18],  describe  a  cable  test  structure  known  as 
SEACON  II.  Figure  6-2  presents  the  pertinent  physical  characteristics  of  the 
structure.  The  static  response  of  the  system  to  an  ocean  current  Is  presented 
in  Figure  6-3.  The  model  was  rather  crude  (2  elements  per  arm  and  eight  elements 
per  leg)  but  the  intent  of  the  problem  was  more  In  demonstrating  the  solution 
procedure  than  in  getting  accurate  results.  The  structure  is  quite  flexible  and 
it  offers  a  good  demonstration  of  the  difficulties  encountered  in  getting  a 
stable  dead  loaded  configuration.  The  natural  starting  point  for  the  problem  is 
a  compatible  arrangement  of  the  unstretched  cables.  In  the  final  dead  loaded 
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FIGURE  6-1  TWO-BAR  TRUSS  SNAPTHROUGH  RESPONSE 
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TABLE  6-1  INPUT  DATA  FOR  TWO-BAR  TRUSS 


WATER  DEPTH .  2860  ft 

DEPTH  AT  BUOYS .  500  ft 

DISTANCE  BETWEEN  ANCHORS... 6600  ft 

LENGTH  OF  EACH  LEG . 4080  ft 

LENGTH  OF  EACH  .■\RM .  1000  ft 

CABLE  DIAM .  0.060583  ft  (0.727  in.) 

CABLE  WEIGHT  IN  SEA  WATER....0.  31  Ib/ft 

CABLE  EA . . .  1.  7735  x  10°  lbs 

NET  BUOYANT  FORCE . 1745  lbs 

(EACH  BUOY) 


BUOY  DIAM 


5,5833  ft 


FIGURE  6-2  THE  SEACON  11  STRUCTURE 


DISPLACEMENTS  EXAGGERATED  FOR  VISUALIZATION 
(XX)  VERTICAL  DISPLACEMENT  +  UP 
DISPLACEMENTS  ARE  IN  FEET 


Figure  6-3.  Deflection  of  SEACON  II  Delta 
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TABLE  6-2  INPUT  DATA  FOR  SEACON  II  STATIC  RESPONSE 


configuration  the  buoys  are  nearly  200  feet  below  their  position  in  the  unstretched 
(straight  line)  configuration.  This  means  they  must  drop  almost  10%  of  their 
initial  height. 

The  input  data  listed  in  Table  6-2  shows  a  procedure  where  the  dead- 
loaded  configuration  was  obtained  by  assuming  initial  tensions  (Method  1  of 
Section  3.5)  and  proceeding  with  a  10  step  DEAD  solution  using  the  RFB  method. 

This  was  followed  by  a  MNR  correction  to  that  configuration.  The  current  was 
then  applied  in  a  LIVE  subanalysis  using  the  MNR  method. 

6.3  THE  TOWED  BODY  PROBLEM 

Gibbons  and  Walton  [19]  report  experimental  data  on  the  steady-state 
configuration  of  a  body  towed  at  sea  by  a  single  bare  cable.  Figure  6-4  shows  the 
results  of  a  finite  element  analysis  of  the  problem.  The  data  shown  in  the 
column  marked  "FLUID  LOAD  ERROR"  represents  a  solution  with  the  fluid  drag  load 
calculated  from  an  erroneous  application  of  Morrison's  Equation  and  the  independence 
principle.  The  towed  body  was  modeled  as  an  equivalent  sphere  with  the  data  shown 
in  Figure  6-4. 

The  dynamic  response  of  the  body  to  a  sudden  turn  is  shown  in  Figure  6-5. 
The  model  used  in  this  case  had  fewer  elements  than  was  used  for  the  results  in 
Figure  6-4.  The  input  cards  for  the  dynamic  case  are  presented  in  Table  6-3. 

The  analysis  proceeded  from  a  simple  DEAD  subanalysis  with  the  line  straight  down, 
to  a  LIVE  subanalysis  to  get  the  steady-state  solution  and  then  to  a  DYN  sub¬ 
analysis  with  a  dynamic  initial  condition.  Table  6-3  also  shows  the  TVARY  subroutine 
used  to  produce  the  course  change. 

6.4  ANCHOR  SUSPENSION  DYNAMICS 

The  response  of  a  system  to  imposed  motion  is  demonstrated  on  the 
system  shown  in  Figure  6-6,  with  the  results  shown  in  Figure  6-7.  The  input 

cards  are  given  in  Table  6-4.  The  TVARY  subroutine  shown  in  the  table  produces  the 

sinusoidal  motion. 
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TABLE  6-3  INPUT  DATA  FOR  TOWED  BODY  TURN 
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BUOY:  BUOYANT  FORCE  11, 200  LBS 

DIAMETER  6. 0  FT 

ANCHOR:  WEIGHT  (TS  WATER)  25, 000  LBS 

DIAMETER  5, 0  FT 


UNSTRETCHED  WEIGHT 

LENGTH  EA  IN  WATER  DIAM 

LINE  (ft)  (lbs)  (Ib/ft)  (ft) 

1  .  4000  7  X  10® . 0  - 0.25 

2  _ 6500  7  X  10® . 1.  0 - 0.  0917 

3  . 7000  6  X  10® _ 0 - 0.167 
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Figure  6-6.  Description  of  Anchor  Suspension  System 
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TABLE  6-4  INPUT  DATA  FOR  ANCHOR  SUSPENSION  DYNAMICS 


6.5 


STATIC  EXCURSION  OF  DD-692  IN  A  FOUR  POINT  MOOR 


A  finite  element  model  of  the  DD-692  ship  in  a  four  point  slack  moor 
is  shown  in  Figure  6-8.  The  position  of  the  ship  relative  to  its  quiescent 
position  is  plotted  in  Figure  6-9.  The  loading  was  a  wind  of  30  knots  and  a 
surface  current  of  2  knots.  They  were  assumed  to  have  coincident  headings 
and  were  applied  initially  as  beam  loading  and  varied  incrementally  to  head 
loading  on  increments  of  ten  degrees.  The  quiescent  state  was  recovered  with 
a  restart  and  loading  was  varied  from  beam  to  following  seas  in  ten  degree 
increments.  The  input  cards  are  presented  in  Table  6-5. 
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TABLE  6-5  INPUT  DATA  FOR  DD-692  EXCURSION 
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APPENDIX  A 


LOGIC  OVERVIEW  OF  THE  SEADYN/DSSM  PROGRAM 

The  SEADYN/DSSM  Program  consists  of  a  simple  main  program  and  a  set 
of  more  than  50  subroutines.  The  main  program  simply  establishes  the  size  of 
unlabeled  COMMON  and  calls  the  driving  routine  (SEADYN)  which  controls  the 
subdivision  of  unlabeled  common  into  arrays  which  are  sized  to  meet  the  require¬ 
ments  of  the  specific  problem  being  analyzed.  The  actual  execution  of  the 
various  options  of  the  program  are  controlled  from  the  MANIPR  subroutine  which 
is  called  by  SEADYN.  The  subanalyses  were  developed  essentially  in  modular 
form  with  the  majority  of  the  interfacing  done  through  MANIPR. 

Specific  information  about  each  of  the  subroutines  and  macro-flow 
charts'  of  each  of  the  subanalyses  are  provided  in  the  Program  Maintenance  Manual 
[5].  A  macro-flow  chart  of  the  overall  analysis  flow  is  presented  in  Figure  A-1. 
The  first  block  corresponds  to  the  main  program.  The  next  three  levels  are 
contained  in  the  SEADYN  subroutine.  The  remainder  represents  the  MANIPR  subroutine 
and  subroutines  it  calls. 
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FIGURE  A-1  MACRO  FLOW  CHART  OF  SEADYN/DSSM 


APPENDIX  B 


SHIP'S  LOADING  FUNCTIONS 


This  Appendix  is  an  excerpt  from  the  DESMOOR  Program  user's  manual  [9]. 
It  describes  the  tabular  and  analytical  approaches  to  obtaining  static  loads 
on  a  ship  subjected  to  winds  and  surface  currents.  Identical  procedures  are 
used  in  SEADYN/DSSM  and  DESMOOR  insofar  as  the  ships  loads  are  concerned.  The 
ship  loading  file  generated  by  either  one  of  the  programs  can  be  used  by  the 
other. 

The  references  for  this  Appendix  are: 

3.  Anon.,  "Design  Manual iHarbor  and  Coastal  Facilities,"  Dept,  of  Navy, 
Naval  Facilities  Engineering  Comnand,  Washington,  DC,  NAVFAC  DM-26, 

July  1968. 

4.  Hughes,  G,,  "The  Air  Resistance  of  Ship's  Hulls  with  Various  Types 
and  Distributions  of  Superstructures,"  lESS,  1932. 

5.  Rossell,  H.E.  PRINCIPLES  OF  NAVAL  ARCHITECTURE,  SNAME,  1939. 

6.  Altman,  R. ,  "Forces  on  Ships  Moored  in  Protected  Waters,"  HYDRONAUTICS, 
Inc.,  Tech  Rept.  7096-1,  July  1971. 
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The  external  loadsapplied  to  the  mooring  system  are  all  assumed  to  be 
applied  through  the  ship.  Three  sources  of  loading  are  considered:  wind  loads, 
surface  currents,  and  point  loads  representative  of  working  loads.  The  point 
loads  are  specified  by  giving  the  coordinates  of  the  point  where  the  loading  is 
applied  and  then  giving  the  components  of  the  lateral  load  and  the  yaw  moment. 

It  is  assumed  that  the  point  loads  do  not  change  their  magnitude  or  global 
directions  as  the  ship  moves  to  a  new  position.  Specification  of  the  wind  and 
current  loadings  is  somewhat  more  complicated.  The  DESMOOR  program  provides  two 
approaches  to  defining  these  loads.  The  first  is  in  the  form  of  loading  tables 
which  give  load  coefficients  versus  ship's  heading  relative  to  the  flow.  The 
second  approach  uses  approximate  analytical  expressions. 

The  tabular  approach  is  based  on  the  procedures  given  in  NAVFAC  Design 
Manual  26  (DM-26)  [3].  The  DM-26  approach  utilizes  experimental  curves  for  the 
forces  and  moments  for  various  headings  of  wind  and  current  for  a  set  of 
"representative"  vessels.  Similarity  scaling  is  then  applied  to  get  loading  values 
for  ships  other  than  the  test  models. 

The  DM-26  procedure  begins  with  a  set  of  load  measurements  obtained  from 
subscale  tests  on  a  representative  ship's  model.  The  measurements  give  values 
for  the  lateral  and  longitudinal  forces  and  yaw  moment  versus  flow  heading.  These 
measurements  represent  the  combined  effects  of  such  phenomena  as  profile  and 
friction  drag,  lift  induced  side  forces,  and  shifts  in  the  center  of  pressure. 
Tables  of  these  measurements  can  be  specified  as  input  to  the  DESMOOR  program  or 
a  special  ship  loading  file  can  be  generated  and  saved  for  use  with  both  the 
DESMOOR  and  SEADYN  programs. 

Given  the  headings  of  wind  and  surface  current  relative  to  the  ship,  the 
loading  coefficients  are  obtained  by  linear  interpolation  in  the  tables.  In  the 
event  that  there  are  tables  provided  for  more  than  one  velocity,  the  table  for 
the  velocity  nearest  the  one  specified  in  the  analysis  will  be  used.  This  is 
determined  by  comparing  the  squares  of  the  velocity  ratios. 
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After  the  load  coefficients  are  obtained  from  the  tables  they  must  be 
scaled  to  account  for  differences  in  the  conditions  modeled  in  the  test  and  those 
being  analyzed.  The  scaling  accounts  for  differences  in  flow  velocity,  water 
depth  and  ship  geometry. 


The  formulas  for  adjusting  for  the  effects  noted  above  are  given  in  DM-26, 
and  they  are  restated  here  for  completeness. 


WIND 


Fc  =  A  /A. 

s  f  ms  s  ts 


e  f  me  e  te 


=  C  V‘M  A  L/A.  L. 
w  m  m  s  ts  t 


(2-20) 

(2-21) 

(2-22) 


where 

F^  =  lateral  force  on  ship 

Pg  =  longitudinal  force  on  ship 

=  yawing  moment  on  ship 

F„„  =  lateral  force  on  model 

ms 

F  =  longitudinal  force  on  model 
me 

=  yawing  moment  on  model 
V  =  wind  velocity 

A^  =  side-projected  areas  above  the  water  line  of  ship  being  analyzed 
Afj  =  side-projected  area  above  the  water  line  of  modeled  ship 
Ag  =  end-projected  area  above  the  water  line  of  ship  being  analyzed 

Afg  =  end-projected  area  above  the  water  line  of  modeled  ship 

L  ■  Length  of  ship  being  analyzed 


=  Length  of  modeled  ship 


“  ~~z 

v/ 


''t 


S  =  linear  scale  of  the  model 

Vy  =  wind  velocity  used  in  model  test 


CURRENT 

^2  ^  ^1  ‘-W2^^W1 


'^1  ^  ''2  r-Wl^'-W2 


^2  '  *^1  ^2'-W2^^l'-Wl 


(2.23) 


(2.24) 


(2.25) 


(2.26) 


(2.27) 


(2.28) 


V>'i. 

f  A  i.*» 


where 

h  =  depth  of  water 
V  =  velocity  of  current 

=  water  line  length  of  vessel 
F  =  lateral  or  longitudinal  resisting  force 
M  =  yaw  resisting  moment 
A  =  displacement 

Subscript  1  denotes  the  vessel  for  which  the  test  was  made,  and  subscript  2 
denotes  the  vessel  being  analyzed. 

When  the  velocity  from  Equation  (2.26)  does  not  correspond  to  one  of  the 
tables  given  for  the  model  test  then  the  forces  and  moments  will  be  selected  from 
the  tables  corresponding  to  the  velocity  nearest  the  value  of  in  Equation  (2.26), 

It  will  then  be  necessary  to  adjust  the  values  by  the  square  of  the  ratio  of  the 
V,|  velocity  and  the  velocity  represented  in  the  tables, 


It  is  quite  likely  that  the  depth  at  the  proposed  mooring  site  will  not  be 
the  same  as  that  obtained  for  hg  in  Equation  (2.25).  In  that  event,  a  correction 
for  depth  is  required.  DM- 26  suggests  that  the  correction  be  made  assuming  an  in¬ 
verse  relationship  with  the  side  resistances  at  the  two  depths  in  question.  The 
curves  given  in  Graph  124  (EC-2)  of  DM-26  are  used  along  with  Equation  (2.25)  for 
this  purpose.  The  data  is  given  in  tabular  form  and  the  side  resistances  are  ob¬ 
tained  by  logarithmic  interpolation.  The  resistance  for  a  depth  greater  than  that 
in  the  table  will  be  the  last  value  in  the  table. 


The  adjustments  for  current  velocity  and  depth  are  summarized  by  the  following 
equations: 

Vi^ 

''s2  *  ^h  7^  ^2  (2.29) 

''tl 


';2 


M- 


'h  A 


=  f.  M 


el 


1 

T 


'tl 


^1  Si  ^  '^ti 


2  ‘^p2  S  h 


'7^ 

^ti 


(2.30) 

(2.31) 


where 


fjj  =  the  depth  scaling  factor 

V^1  =  the  velocity  at  which  the  test  data  was  obtained 
A  =  the  propeller  projected  area 
Cp  =  the  propeller  drag  coefficient 


The  primes  indicate  the  value  adjusted  to  the  desired  conditions  for  the  mooring 

site.  Equation  (2.30)  reflects  the  adjustment  in  the  longitudinal  force  recommended 

by  DM-26  with  the  assumption  that  _L  =  2.88  with  A  in  ft^  and  V  in  knots. 

Assuming  the  specific  weight  of  sea^wateP  is  64  Ib/ft^  and  the  acceleration  due  to 

gravity  is  32.2  ft/sec^  then  C_  =  1.00.  The  form  using  pC  rather  than  2.88  is 

P  2  ^ 

required  to  make  the  procedure  dimensionally  independent. 
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The  approximate  analytical  expressions  for  ship's  loading  are  based 
primarily  on  the  work  of  Hughes  [4]  and  standard  naval  architectural  formulas  [5]. 
The  wind  loading  is  given  by 

F  =  K  pg  (A^  sin^e  +  cos^e)  cos{a-6)  (2-32) 

where 

K  =  constant,  0.6 

F  =  resultant  wind  force 

p,  =  density  of  air 
0 

V  =  wind  velocity 

6  =  wind  heading  relative  to  the  bow 

a  =  heading  of  the  resultant  wind  force  relative  to  the  bow 

The  heading  of  the  resultant  force,  o,  is  approximated  as  a  function  of  e 
in  a  7^^  order  polynomial  as  follows: 

a  =  0.0715608  +  7. 954381 6  -  0.32545616^ 

+  0.00731316^  -  9.3966  x  10‘^6^ 

+  6.85008  X  10"^^  -  2.6323  x  10‘^6® 

+  4.1453  X  10’^^6^  (2.33) 

In  Equation  (2.33)  both  6  and  a  are  measured  in  degrees. 

The  distance  between  the  ship  forward  perpendicular  and  the  center  of 
wind  pressure,  X^p,  can  be  approximated  as  a  polynomial  function  of  the  wind 
direction,  e.  This  relationship  is 

-P  =  0.2004112  +  0.00486416  -  4.52442  x  lO'^e^ 

+  5.45736  x  10‘^e^  -  3.78789  x  10"®e^ 

+  1.02881  X  lO’^^e^ 


(2.34) 


Here,  as  above,  e  is  measured  in  degrees.  The  yawing  moment  due  to  wind  is  then 
approximated  by 


M„  =  FL 


sin  a 


(2.35) 


Analytical  expressions  for  the  resistances  from  current  effects  utilize 
the  approach  presented  by  Altman  [6].  These  expressions  are  summarized  below: 


^  (h/H)^-l. 

(2.36) 

'  0.215  p,y  L  H  sin  9 
w  w 

(2.37) 

1  „/  (S„C„  +  ACp)  cos  e 

(2.38) 

'^s 

(2.39) 

where 

7  =  displaced  volume 

Cj  =  wetted  surface  coefficient,  input  with  ship  description 

Cn  =  hull  resistance  coefficient,  input  with  ship  description  or 

calculated  as  +  0.0005  (2.40) 

=  residuary  resistance  coefficient 

Cr  =  frictional  resistance  coefficient,  =  — — (2.41) 
Rg  =  Reynolds  number  for  the  hull 

M  =•  Fg  L^p  (2.42) 

Lj.p  =  distance  from  midships  to  hull  center  of  pressure 

Lcp  =  L  [LgQ  +  0.00226  (  9  -  90»}]  for  0*’  <  9  <  180” 

irn  =  L  n:„„  +  0.00226  (9  -  270”)]  for  180”  <  9  <  360” 


(2.43) 


Lqq  =  ratio  of  distance  to  center  of  pressure  at  e  =  90®  to  the  distance 
to  the  center  of  hull  side  area 

Various  of  these  terms  require  further  discussion.  The  hull  resistance 
coefficient ,  Cj^,  represents  the  sum  of  various  coefficients  for  different 
sources  of  hull  resistance.  This  coefficient  may  be  input  or  calculated  in  the 
computer  program.  When  no  input  is  given  for  it  will  be  calculated  as  the 
sum  of  a  residuary  resistance  coefficient,  a  frictional  resistance  coefficient, 
and  a  fouling/surface  effect  coefficient.  The  fouling/surface  effect  coefficient 
will  be  given  an  arbitrary  value  of  0.0005.  T>»e  frictional  coefficient  will  be 
calculated  from  Equation  (2.41)  and  the  residuary  coefficient  will  be  obtained 
from  linear  interpolation  in  a  digitized  form  of  Figure  38  of  the  Altman  report  [6]. 
This  method  of  obtaining  is  limited  to  low  flow  velocities  since  wave-making 
resistances  are  ignored. 

The  longitudinal  location  of  the  center  of  pressure  for  a  hull  skewed  with 
respect  to  the  flow  is  estimated  by  Equations  (2.43).  This  requires  an  estimate 
of  the  ratio  of  the  distance  to  the  center  of  pressure  and  the  center  of  area  for 
beam  flow,  Ugg.  This  factor  will  be  estimated  by  linear  interpolation  between 
the  values  for  the  ships  DD-692  and  EC-2  using  the  block  coefficient  as  a  reference. 
Altman  gives  the  values  for  Ugg  for  the  DD-692  and  EC-2  as  D.056  and  -0.138, 
respectively.  (Negative  means  aft  of  midships.) 

It  should  be  emphasized  that  these  analytical  expressions  are  to  be 
viewed  as  a  convenient  alternative  to  the  DM-26  experimental  curve  procedure. 

They  are  not  represented  as  highly  accurate,  and  it  remains  to  be  demonstrated 
that  they  are  capable  of  giving  reliable  approximations  of  the  ships  loading. 
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APPENDIX  C 

DESCRIPTIONS  OF  THE  SHIP  LOAD  FILE 
AND  COMPONENT  INVENTORY  DATA 

This  Appendix  deals  with  the  special  storage  requirements  of  the 
Ship  Load  File  (logical  unit  10)  and  the  arrays  for  the  mooring  component 
Inventories. 

The  Ship  Load  File  can  be  setup  and/or  used  by  either  the  DESMOOR  or 
the  SEADYN/DSSH  Program.  The  file  contains  one  logical  record  for  each  ship 
cataloged  on  the  file  and  It  Is  written  In  a  binary  form  with  the  following 
FORTRAN  statement: 

WRITE(10)NWIND.NTHETW,WNDVEL,WNDHED.WN0COE,SCALE,NCRNT,NTHETC,CURVEL,CURHED, 

CURCOE , TOEPTH , TBLOCK , TSLT, TSAE . TSAS . TSWL , TSB , TSD , TSDSP . TSAP , SHCAP , WLBL , CLBL , LLBL , VLBL 

A  number  of  the  Items  In  the  list  are  arrays  and  they  are  written  In 
their  entirety  using  the  Implied  OO-LOOP  feature  of  FORTRAN  I-O  statements.  A 
description  of  each  Item  Including  the  dimensions  of  the  arrays  Is  given  below. 


VARIABLE 

DESCRIPTION 

NWINO 

No.  of  wind  velocity  tables 

NTHETW 

No.  of  headings  1n  each  wind  table 

WNDVEL{5) 

Array  of  wind  velocities 

WNDHE0(20) 

Array  of  wind  headings 

WN0C0E(20.3,5) 

Array  of  wind  load  coefficients  giving 
values  for  up  to  20  headings  for  end 
force,  side  force  and  yaw  moment  for 
up  to  5  wind  velocities. 

SCALE 

Scale  for  wind  load  tests  (A  means  1/A^^  ; 

NCRNT 

No.  of  current  velocity  tables 

NTHETC 

No.  of  headings  in  each  current  table 

CURVEL(5) 

Array  of  current  velocities 

CURHED(20) 

Array  of  current  headings 

CURCOE (20. 3, 5) 

Array  of  current  load  coefficients  giving 
values  for  up  to  20  headings  for  end 
force,  side  force  and  yaw  moment  for  up 
to  5  current  velocities. 

-162- 


VARIABLE 


DESCRIPTION 


TDEPTH  Water  depth  for  test 

TBLOCK  The  ship's  block  coefficient 

TSLT  Total  ship  length 

TSAE  End  projected  wind  area 

TSAS  Side  projected  wind  area 

TSWL  Water  line  length 

TSB  Beam  at  mid-ships 

TSD  Draft  at  mid-ships 

TSDSP  Volume  displacement 

TSAP  Propeller  projected  area 

SHPCAP(12)  _  Title  of  12. six-character  Hollerith  words 

WLBL  Wind  force  label,  six-character  Hollerith  word 

CLBL  Current  force  label,  six-character  Hollerith  word 

LLBL  Length  label,  six-character  Hollerith  word 

VLBL  Velocity  label,  six-character  Hollerith  word 

The  mooring  component  Inventories  are  contained  In  arrays  which  are 
generated  by  DATA  statements  In  the  COMPOAT  subroutine.  The  data  is  stored  in 
labeled  COMMON  as  follows: 

COMMON/COMPNT/  ATYPE(3,6),ANCTAB(5,16,6),NANCR(6),BTYPE(3,2),BU0TM(7,6,2), 
NBU0Y{2).HTYPE{3,4),HAWTAB(2.23.4),HAWSIZ(23),NHAWS(4),CHAIN(31,3),NCKN, 

HLDFAC(6) 


VARIABLE 

ATYPE 

ANCTAB 


DESCRIPTION 


Label  for  each  of  6  ,  ichor  types.  Each 
label  Is  3  six-character  Hollerith  words. 

Storage  for  anchor  data  for  up  to  6 
anchor  types  (3rd  dimension).  The  array 
allows  up  to  16  anchors  of  each  type, 
listed  In  order  of  increasing  weight. 

Five  Items  are  given  for  each  anchor. 

WORD  CONTENTS 

1  Anchor  weight 

2-4  Hollerith  labels  for  federal 

stock  number 

5  Holding  power 

The  holding  power  is  calculated  from  the 
weight  and  the  factors  given  in  HLDFAC. 


VARIABLE 


DESCRIPTION 


NANCR  An  array  listing  the  number  of  anchors 

of  each  type  in  the  inventory. 

BTYPE  Label  for  each  of  two  buoy  types.  Each 

label  is  3  six-character  Hollerith  words. 

BUOTAB  Storage  for  buoy  data  for  up  to  two  buoy 

types.  The  array  allows  up  to  six  buoys 
of  each  type,  listed  in  order  of  increasing 
buoyancy.  Seven  items  are  given  for  each 
buoy: 

WORD  CONTENTS 

1  Outside  diameter 

2  Height 

3  Weight 

4  Nominal  buoyancy 

5  Maximum  Buoyancy 

6-7  Hollerith  label  for  federal 

stock  number 

NBUOY  An  array  listing  the  number  of  buoys  of 

each  type  in  the  inventory. 

HTYPE  Label  for  each  of  4  hawser  types.  Each 

label  is  3  six-character  Hollerith  words. 

HAWTAB  Storage  for  hawser  data  for  up  to  4  hawser 

types.  The  array  allows  up  to  23  sizes 
for  each  type.  The  two  items  given  for 
each  size  are  the  tensile  strength  and  the 
weight  per  100  units  of  length. 

HAWSIZ  A  list  of  the  hawser  sizes  given  in 

increasing  order.  It  is  assumed  that  all 
four  hawser  types  have  the  same  list  of 
sizes.  If  no  entry  is  available  for  a 
given  size  then  set  the  strength  to  zero 
and  the  weight  to  some  small,  non-zero 
number. 

NHAWS  Set  to  23  for  each  hawser  type. 

CHAIN  Storage  array  for  chain  data.  Presumes 

only  one  chain  type  (stud-link  chain). 

The  array  contains  31  entries  for  size, 
strength  and  weight  per  unit  length. 


VARIABLE 


DESCRIPTION 


NCHN  The  number  of  chain  entries  (31). 

HLDFAC  A  list  of  the  holding  power  factors  for 

the  6  anchor  types.  This  Is  the  number 
which  multiplies  the  anchor  weight  to 
get  the  holding  power  In  firm  sand. 

The  component  data  presently  listed  In  the  Inventories  assumes  weights, 
buoyancies  and  strengths  are  In  units  of  pounds.  Lengths  and  buoy  dimensions  are 
In  feet.  Hawser  and  chain  sizes  are  In  Inches. 

It  would  be  possible  to  alter  the  contents  of  the  Inventories  by 
rewriting  the  COMPDAT  subroutine.  The  subroutine  Is  called  once  each  time  a 
component  check  or  design  selection  Is  made.  It  Is  presently  set  up  to  avoid 
re-calculating  entries  after  the  first  call.  The  program  uses  the  Items  In  the 
arrays  to  search  for  the  appropriate  entries.  The  appropriate  values  are  then 
selected,  scaled  without  altering  the  arrays,  and  used.  Minor  changes  could  be 
made  without  changing  the  calling  program.  More  extensive  changes  or  generalizations 
may  require  some  modifications  (not  major)  In  the  calling  program. 

A  listing  of  the  present  version  of  COMPDAT  and  the  Inventories  It 
contains  are  presented  on  the  following  pages. 
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COMPCNENT  INVENTORY 


r,- v" j 


ANCHOR  TYPE  =  NAVY  STO  STOCKLESS 


HEIGHT 

FED.  STOCK  NO, 

HOLD.  POWER 

.300006*04 

C204fl-5 16-7758 

IfIRH  SAND) 
’,210006*09 

.500006*04 

C204fl-516-775;f 

,350006*09 

.600006*04 

02040-516-7796 

,420006*05 

.700006*04 

C2040-516-7759 

,490006*05 

.900006*04 

C2fl40-516-77S<» 

,630006*05 

w 

.100006*05 

C20  4a-272-2244» 

,700006*05 

.130006*05 

02040-272-2249 

,910006*1)9 

•  145006*05 

02040-272-2246 

•  101506*  06 

.180006*05 

C2040-5 16-7799 

,12^006*06 

.200006*05 

C2040-272-2247 

,140006*06 

.250006*05 

C2Q40-272-2242 

,175006*06 

.300006*05 

C20 40-272-2241 

,210006*06 

.400006*05 

02040-277-242^ 

,280006*06 

w 

ANCHOR  TYPE  = 

NAVSHIP(LHT) 

HEIGHT 

FED.  STOCK  NO. 

HOLD.  POWER 

w 

.  100006*03 

H20 40-377-8608 

fflRH  SAND) 
,283746*04 

.150006*0  3 

H2040-377-86ai 

,395656*04 

w 

.200006*03 

H2040-377-6602 

,500916*04 

.300006*03 

H2040-377-8603 

,698486*08 

.900006*03 

H2040-377-8604 

,106196*05 

Si^ 

.750006*03 

H2040-377-86flS 

,148076*05 

.100006*04 

H2040-377-86a6 

,187466*05 

.  150006*0  4 

H2040-377-8607 

,261406*09 

'w 

.200006*04 

H204fl-377-8b08 

,330956*05 

.290006*04 

H2040-377-8609 

,397406*05 

.300006*04 

H2040-377-8610 

,461486*05 

.400006*04 

H2040-377-8611 

,584266*05 

.900006*04 

H2040-378-5633 

,701576*05 

.600006*04 

H2040-378-5634 

,814706*09 

.  100006*05 

H2040-377-8612 

,123856*06 

.130006*05 

H2040-377-8613 

,153986*06 

ANCHOR  TYPE  =  NAVFAC  STATO 


HEIGHT 

.200006*03 
.30000E*04 
.600006*04 
.900006*04 
. 1ZQ00E*0S 
.  190006*09 


F60.  STOCK  NO, 

2CF2040-aC0-9699 

2CF2a40-702-7864 

2CF2040-702-678S 

2CF2040-702-6788 

2CF2040-702-6787 

2CF2040-80t-7938 


NOLO.  POWER 
(FIRH  SAND) 
,400006*04 
,603  006*05 
,120006*06 
,180006*06 
,240006*06 
,300006*06 


I  % 

s"* 

•  •  K  *  .  N  • 

*  •  4  .  4  •  'v 

.'S-vS 


.  * . ' 


BUOY  DATA 


STEEL  STUO-LIKK  CHAIN 


SIZE 

STRENGTH 

•  7500 

48550. 

.8750 

65280. 

1.0000 

84500. 

1.1250 

106080. 

1.2500 

130070. 

1.3750 

156330. 

1.5000 

185060. 

1.6250 

216030. 

1.7500 

249210. 

1.8750 

284540. 

2.0000 

322000. 

2.1250 

3  61530. 

2.2500 

403100. 

2.3750 

446660. 

2.5000 

492190. 

2.6250 

539620. 

2.7500 

588930. 

2.8750 

640070. 

3.0000 

693000. 

3.1250 

747680. 

3.2500 

804070. 

3.3750 

862130. 

3.5000 

921810. 

3.6250 

983080. 

3. 7500 

10 45900. 

3.8750 

1110210. 

4.0000 

1176000. 

4.1250 

1234200. 

4.2500 

1311790. 

4.3750 

1381330. 

4.5000 

1452  930. 

HEIGHT /LENGTH 
5.5556 
7»7T7i 
9.6444 
12.2222 
15.0000 

17.7778 
21.1111 
24.4444 

28.3333 

32.7778 

36.6667 

41.1111 

46.6667 

51.6667 

57.7778 

63.3333 
70.0000 

76.6667 

83.3333 

91.1111 

98.3333 

106.1111 
114.4444 

122.7778 

131.1111 
140.0000 

148.8889 

158.8889 
170.0000 

183.8889 

198.3333 


Code  1  Code  2  Code  3  Code  4 

2-IN-t  NYLON  2-IN-l  POWER  9RAI0  2-IN-t  STRBL  BRAID  12  ST.  OlUC  STREIR 

sue  STRCNCIH  NEICHI/ISCL  strength  4EICHT/10CL  STRENGTH  HEIGHT/lOtL  STRENGTH  HEICNT/IIBL 

.SGOd  SJQO.  6. so  TSOJ.  S.90  TEDC.  T.90  6300.  6. NO 
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APPENDIX  D 


DESCRIPTION  OF  THE  SHIP  MOTION  FILE 

The  data  for  the  motion  equations  for  a  ship  driven  by  harmonic  waves  is 
provided  to  the  SEADYN/DSSM  program  through  the  Ship  Motion  File.  This  set  of 
data  is  assumed  to  be  on  a  sequential  binary  file  on  logical  unit  08.  This 
Appendix  describes  the  format  of  that  file.  The  notations  of  Reference  4  are 
used. 


The  equations  of  motion  for  a  ship  moving  in  waves  on  a  free  surface 
are  assumed  to  have  the  following  form: 

6  .  iWpt 

("jk * ''jkK  *  Vk *  Vk  ■''j*  -J’’ . ® 

The  terms  of  this  equation  are  assumed  to  be  provided  on  the  Ship  Motion  File  in  a  non 
dimensional  form  reflecting  the  effects  of  the  wave  length  of  the  surface  wave 
and  the  relative  heading  between  the  wave  and  the  ship. 

The  relationships  between  the  terms  of  equation  (D-l)  and  the  non- 
dimensional  terms  on  the  file  are  given  below: 

Mj[j  »  M  .  .  GMU(J,K)  (0-2) 

Units:  FT^ 

Ajk  ”  .  DA(J,K)  (D-3) 

Units:  FT^ 

=  M  .  .  DB(J,K)  (D-4) 

Units:  FTl""^ 

B*4  =  M  .  Y57l  .  .  B44S  {l^) 

Units:  FTL 

r =  M  n  I  n( jfk)-l  nrf.i.Kl 


(D-5) 


fD-6) 


A 


(B0D{J,1)  +  iB0D(J+3,l))  J-1,3,5 

P  E  M  n  I  j  )"^ 

*^3  ^  '  9  ■  I  •  (BEV(J,1)  +  iBEV(J+3,l))  J-2,4,6 


(D- 


Units: 

where: 

M  ■  Ship's  mass  (TMAS) 

Units:  FT^  L"'' 

L  =  Ship's  length  (ELL) 

Units:  L 

g  =  gravitational  acceleration  (6RAV) 

Units;  LT"^ 

=  the  Ith  roll  angle  index 
n(j,k)  =  m(j)  +  m(k) 

m( j )  =0  for  j  ^  3 
1  j  >  3 

i  =  /T 

The  coefficients  are  assumed  to  be  linearized  for  a  unit  motion  anplitude  and  wave 
height.  The  usual  units  for  the  Ship  Motion  File  are: 

F  -  long  tons  (2240  lbs) 

L  -  feet 

T  -  seconds 

Angles  and  angular  responses  are  assumed  to  be  in  radians. 

In  addition  to  the  ship  motion  coefficients  the  Ship  Motion  File  provides 
coefficients  for  estimating  the  steady  state  approximations  of  the  second  order, 
wave  induced  drift  forces.  These  dimensionless  coefficients  are  used  to  estimate 
the  drift  forces  once  the  magntiudes  of  the  ship  responses,  ,  are  obtained. 


1 

This  presumes  the  wave  motion  is  given  by  e  where  n.  is  real  and  the  ship 

iwt  ^  ° 

response  is  e  . 


The  coordinate  system  presumed  for  ship's  motions  and  forces  is  a  right- 
handed  cartesian  system  with  its  origin  at  the  ship's  center  of  gravity,  its  x  axis 
positive  aft,  its  y  axis  positive  starboard,  and  its  z  axis  positive  upward.  The 
angular  convention  for  the  relative  heading  between  the  ship  and  the  waves  assumes 
the  following: 


Wave  Heading  Description 

0°  Following  seas 

90°  Beam  seas  with  waves  traveling  from 

port  to  starboard 

1 80°  Head  seas 

270°  Beam  seas  with  waves  traveling 

from  starboard  to  port 

The  differences  between  the  wave  heading  convention  and  the  ship's  coordinates 
should  be  noted. 

The  Ship  Motion  File  is  organized  in  logical  records.  The  specific  con¬ 
tents  of  each  record  will  be  described  below.  There  are  seven  distinct  record 
types.  The  first  two  records  contain  data  which  is  independent  of  wave  heading 
or  wavelength.  Record  types  3  through  7  are  dependent  on  heading  and  wavelength 
and  are  repeated  in  a  nested  loop  fashion.  The  overall  form  is: 

Record  1 
Record  2 


Length 


Heading 


Record  3 
Record  4 
Record  5 
Record  6 
Record  7 


The  wave  headings  are  assumed  to  be  listed  in  decreasing  order  with  +180°  being 
the  largest  allowed.  The  interpolation  routines  assume  the  values  given  for  +180° 
will  be  used  for  -180°,  therefore  data  for  -180°  need  not  be  given. 
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The  wave  lengths  are  assumed  to  be  listed  In  decreasing  order  (i.e.,  increasing 
frequency  order). 

The  individual  records  of  the  file  are  described  in  terms  of  the  Fortran 
read/write  lists  associated  with  each  record. 

RECORD  1  NAME  1,  NAME  2,  NAME  3 

Three  Hollerith  variables  providing  identifying  data. 

RECORD  2  (TITO(I),  1-1.  12),  WORD,  WORD  2,  WORD  3,  ELL,  BEAM,  DRAFT,  TVOL, 

TMAS,  TPST,  Z6,  CBV,  NOB,  (FN(I),  I»l,  NOB),  NOH,  (HDGl(I),  1*1,  NOH), 

NOK,  (BAM(I),  1=1,  NOK),  VNY,  GRAY,  NRV,  (RANG(I),  1*1,  NRV), 

((GMU(I,J),  J-1,  6),  1=1,6),  ((DC(I,J),  J«l,6),  1=1,6) 

TITO  =  Hollerith  title  consisting  of  12  6-character  words 
WORD  =  length  unit  label  (6-character  Hollerith) 

WORD  2  ■  force  unit  label  (6-character  Hollerith) 

WORD  3  =  moment  unit  label  (6-character  Hollerith) 

ELL  ■  Ship's  length  (L) 

BEAM  =  Beam  (L) 

DRAFT  «  Draft  (L) 

TVOL  =  Ship's  volume  is  obtained  from  (ELL/2)^  .  TVOL 
TMAS  *  ship's  mass  (FS^L’b 

TPST  =  Longitudinal  distance  from  c.g.  to  forward  most 
station  is  obtained  from  (ELL/2)  .  TPST 
ZG  =  vertical  distance  from  water  line  to  c.g.,  (+  up)  (L) 

CBV  =  vertical  distance  from  water  line  to  center  of  buoyancy  (+  up) 
is  obtained  from  ELL  .  CBV 

NOB  =  number  of  speeds  (SEADYN/DSSM  expects  only  one) 

FN(I)  *  the  Froude  numbers  for  each  speed 
NOH  *  number  of  wave  headings 

HDGI(I)  =  the  wave  headings  listed  in  decreasing  order  starting  with 
180°  and  proceeding  no  further  than  -180°  +  40 
NOK  *  number  of  wave  lengths 

DAM(I)  =  non-dimensional  wave  length  in  decreasing  order, 

X  =  ELL  .  BAM(I) 


VNY  s  fluid  viscosity  (L^T) 

_2 

GRAV  =  gravitational  acceleration  (LT  ) 

NRV  *  number  of  roll  angles 

RANGE(I)  the  values  of  roll  angles  (radians)  listed  in 
increasing  order 

GMU(I,J)  =  the  non-dimensional  mass  matrix 

DC(I,J)  s  the  non-dimensional  hydrostatic  restoring  matrix 

RECORD  3  W.  HDGI(MM),  JJ,  FN(JJ).  LL,  BAM(LL) 

=  heading  number 
JJ  speed  number 
LL  =  wave  length  number 

RECORD  4  ((DAd.J),  J=1.6).  1=1,6),  ((DB(I,J),  J=l,6),  1=1,6) 

DA(I,J)  =  the. non-dimensional  added  mass  matrix  for  the  current 
combination  of  heading,  speed,  and  wave  length 
0B(I,J)  =  the  non-dimensional  wave  damping  matrix 

RECORD  5  (BOD(I),  B0D(l+3),  BEV(I),  BEV(I+3),  1=1,3) 

BOD,  BEV  =  the  non-dimensional  wave  force  coefficients. 

RECORD  6  (B44S(I),  1=1,  NRV) 

B44S(I)  =  the  non-linear  roll  damping  terms  which  are  added  to 
the  damping  matrix  depending  on  the  size  of  the  roll 
angle  (non-dimensional). 

RECORD  7  ((TX(I,J),  1=1,7),  J=l,7),  ((TY(I,J),  1=1,7),  J=l,7),  ((TM(I,J).  1=1,7), 
0*1.7),  (TP(I).  1=1,7) 

The  non-dimensional  drift  force  coefficients. 


APPENDIX  E 


LINEARIZED  FLUID  DAMPING  ON  CABLE  ELEMENTS 

An  approximate  form  of  the  fluid  damping  effects  are  developed  in  this 
Appendix.  To  obtain  this  approximation  it  is  assumed  that  flow  field  induced 
drag  is  accounted  for  separately  in  a  static  analysis  and  only  the  dynamic  re¬ 
sponse  of  the  cable  element  contributes  to  the  drag  damping.  This  assumption 
neglects  a  cross  term  involving  the  local  flow  field  velocity  and  the  element 
velocity  which  arises  from  squaring  the  relative  velocity  [2]  .  A  one  dimen¬ 

sional  example  will  clarify  the  matter.  Given  a  body  moving  parallel  to  the 
flow  in  a  uniform  stream.  If  the  stream  velocity  is  and  the  velocity  of  the 
body  is  Vg,  then  the  drag  loading  on  the  body  is  proportional  to 

(V,  -  Vb)2  .  .  ZV^Vj  *  (E-l) 

The  approximation  introduced  here  amounts  to  saying  that  is  applied  to  the 
system  to  obtain  a  static  reference  position  of  the  body  and  Vg  represents  the 
velocity  of  oscillation  about  that  state.  The  drag  contributions  from  Vg  are 
linearized  and  the  cross  term  -2VgVg  is  neglected. 

Assume  an  element  in  its  local  coordinate  system.  Further  assume  that 
the  drag  force  per  unit  length  has  the  following  linear  form: 

where  is  a  equivalent  linearized  damping  coefficient 

V,  W  components  of  the  element  velocity  in  two  orthogonal  directions 
normal  to  the  element. 

Using  the  finite  element  approximation  with  a  linear  shape  function  leads  to 
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0  0  0  10  0 

0  1  0  I  0  % 

0  0  1  I  0  0 

0  0  0^0  0 

0  %  0  I  0  1 

0  0  %  1  0  0 


(E-5) 


The  definition  of  the  element  damping  matrix  will  be  complete  once  an 
expression  for  Cg  is  obtained.  This  expression  can  be  obtained  by  considering  a 
segment  of  cable  oscillating  harmonically  in  a  transverse  mode  in  a  fluid.  Assum¬ 
ing  the  deflected  shape  is  closely  approximated  by  a  sine  curve  facilitiates  the 
development. 


Y  -  Yq  Sin  wt  sin 

Y  »  Yq  w  Cos  wt  Sin  ^ 


The  normal  drag  per  unit  length  is 


(E-6) 

(E-7) 


:-d 


'yM 


f  =  >.  Pf  d 

The  work  done  in  one  cycle  over  the  length  1  is 

TT 

1  2w 


■7/ 


fdxdy  *  2 


'’f  Sd 


TT 

2w 


/  / 


TT 

2^7 


/  I 


1 


fdxydt 


TT 

2^7 


.  3 

(y)  dxdt 


(E-8) 


(E-9) 


An  equivalent  linear  form  for  the  normal  drag  force  per  unit  length  can  be  written: 


=  C^y 
e  e 


(E-10) 
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The  equivalent  work  term  is; 


r  r  .  ^ 

“e  “  i  j  “  2^e  I  I 


.  2 

(y)  dxdt 


(E-ll) 


Requiring  that  the  linearized  term  must  dissipate  the  same  energy  per  cycle  as  the 
non-linear  term  yields 


(  ^  f' 

h _ lo 


(y)  dxdt 


TT  1 

i:  I 


{E-12) 


(y)  dxdt 


Substituting  from  equation  (E-7)  and  carrying  out  the  integration  yields; 


Ce  =  Pf  C^dwyQ  (E-13) 

These  results  are  applied  in  the  computer  program  by  assuming  y^  =  0.01 
for  calculating  when  starting  from  a  static  state  and  selecting  yQ  to  be  the 
largest  of  the  normal  components  of  the  nodal  velocities  on  each  element  on  sub¬ 
sequent  iterations.  A  minimum  of  two  iterations  are  required. 
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The  equations  for  the  added  mass  and  damping  coefficients  for  mooring 
buoys  which  are  given  in  equations  (2-46)  through  (2-49)  on  page  22  proved 
to  be  cumbersome  for  programming.  A  polynomial  curve-fitting  of  the  original 
curves  given  by  Kim[12]  was  used  in  the  program.  The  polynomial  coefficients 
are  given  below. 


FUNCTION  a 
M..  1 


^0 

"1 

"2 

"3 

®4 

^5 

^6 

«7 

1.0620 

-0.4090 

5.3299 

-10.0143 

7.6387 

-2.9089 

0.5509 

-0.0414 

1 . 7945 

1 . 3362 

-10.9227 

18.0521 

-13.8920 

5.5725 

-1.1249 

0.0902 

0. 

-1.5252 

7.2144 

-8.6447 

5.0535 

-1.6218 

0.2750 

-0.0192 

0. 

4.3747 

-9.8378 

10.7232 

-6.6657 

2.3676 

-0.4439 

0.0339 

The  equation  form  is  F=  i  a.(a‘) 

i=0  ^ 


This  page  added 
JULY  1978 


